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Vorwort 


Dieser Band ist Bestandteil einer Reihe über die Anwendung programmierbarer Taschenrechner in 
Naturwissenschaft und Technik. Er versteht sich nicht als Lehrbuch, sondern als Grundlage und An- 
regung zur Erstellung eigener Programme für den jeweils vorhandenen Rechnertyp (jeglicher Nota- 
tion). Aus dieser Sicht ist auch das breite Spektrum der Anwendungsbeispiele zu sehen. Es ging mir 
bei den Programmen in erster Linie um eine klare, übersichtliche Form und nicht um die Ausnutzung 
bestimmter Typ-Eigenheiten., 

Dieses Buch wendet sich an Ingenieure und Techniker, sowie auch an Studenten der Universitäten 
und Fachhochschulen. Weiterhin soll es als Anregung zum Einsatz des programmierbaren Taschen- 
rechners beim praktischen Physikunterricht in allgemeinbildenden Schulen dienen und gewisse, noch 
herrschende Vorurteile abbauen. Es lassen sich Zusammenhänge demonstrieren, die auch experimen- 
tell nicht darstellbar sind (z.B. die Planetengesetze). Darüberhinaus zeigen sich auf natürliche Weise 
die Beziehungen zwischen Bewegung und mathematischem Gesetz. Aber auch dem interessierten 
Laien wird durch die kurze Einführung zum Themengebiet eine Einarbeitung ermöglicht. Es wird 
außerdem eine gewisse Kenntnis in der Taschenrechnerprogrammierung, insbesondere zu den Typen 
TI 58/59, vorausgesetzt. Eventuell ist das in der Literatur angegebene Einführungsbuch [1] zu 

lesen. 

Der Inhalt dieses Bandes umfaßt die Technische Mechanik mit ihren Teilgebieten Kinematik, der 
Lehre von den allgemeinen Bewegungsvorgängen, und der Dynamik, der Lehre von den Kräften. 
Letzere unterteilt sich wiederum in die Statik, der Lehre vom Gleichgewicht der Körper, und der 
Kinetik, der Lehre von den Körperbewegungen durch Kräfte. Da die Kinematik Grundlagen der 
Kinetik behandelt, die Statik jedoch allgemeine Grundlagen betrachtet, ist in diesem Buch die 
Reihenfolge Statik/Kinematik/Kinetik gewählt worden. Die Programme sind so allgernein gehalten, 
daß sie ein möglichst umfassendes Teilgebiet dieser Gliederung erfassen. 

Durch die Anregung von Herrn H. J. Niclas, Lektor im Vieweg Verlag, entstand dieser Band. Ihm 
und dem Verlag Vieweg möchte ich an dieser Stelle für die freundliche Aufnahme danken. Weiterer 
Dank gebührt Herrn K. Nielsen, Produkt Marketing Manager von Texas Instruments und Herrn K.H. 
Burkart, für ihre hilfreiche Unterstützung. Zuletzt gilt mein besonderer Dank all denjenigen, die 
direkt oder indirekt, zur Entstehung dieses Buches beigetragen haben und damit insbesondere meiner 
Frau (Ulrike) für ihre tatkräftige Unterstützung bei der Manuskriptbearbeitung. 


Hamm, September 1979 Harald Nahrstedt 
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1 Einführung 


Bei der Lösung naturwissenschaftlicher und technischer Probleme, gewinnen neben den analyti- 
schen Methoden die numerischen immer mehr an Bedeutung. Die Entwicklung einer exakten Lö- 
sung (sofern sie überhaupt existiert), ist in der Regel mit erheblich größerem Aufwand verbunden, 
als der Einsatz eines Näherungsverfahrens. Zumal die so erhaltene Näherungslösung der exakten 
Lösung beliebig angenähert werden kann. Die dabei auftretenden umfangreichen Rechnungen er- 
ledigt üblicherweise eine EDV-Anlage. Die zunehmenden Speicher- und Strukturerweiterungen 
programmierbarer Taschenrechner lassen in dieser Hinsicht ihren Einsatz immer interessanter er- 
scheinen. Durch ihre Ortsgebundenheit und den damit verbundenen direkten Einsatz nicht nur 
am Arbeitsplatz, bilden sie eine sinnvolle Ergänzung vorhandener größerer Anlagen. 


1.1 Algorithmen und Fiußdiagramme 


Jedem automatisierten Prozeß liegt ein Algorithmus zugrunde. Umgekehrt ist auch bisher kein an- 
derer Weg bekannt, einen Prozeß zu automatisieren, als ihn zu algorithmisieren. Das heißt endliche, 
linear folgende Regeln festzulegen, nach denen ein vorhandener oder zu konstruierender Automat 
durch Eingabewerte und sinnvolle Umformungen, Ausgabewerte (Ergebnisse) erzeugt. Da ein pro- 
grammierbarer Taschenrechner ein ebensolcher Automat ist, bedarf es zu seiner Nutzung solcher 
Algorithmen. Die exakte Formulierung eines Algorithmus in der Weise, daß sie vom Rechner „ver- 
standen’’ und nachvollzogen werden kann, nennt man Programm. Damit haben wir alle Stufen der 
Programmentwicklung angedeutet. 

Sie beginnt bei der Problemanalyse und der Feststellung durchzuführender Regeln. Den so, mitunter 
schriftlich fixierten Algorithmus, kann man mit Hilfe eines Flußdiagramms graphisch anschaulich 
wiedergeben. Wie immer, wenn die Umgangssprache unzureichend ist, bedient man sich einer spe- 
ziellen Sprachform. So ist es in der Technik die Technische Zeichnung und in der Informatik das 
Fiußdiagramm. Es besteht vorwiegend aus den in Bild 1.1 gezeigten einfachen Sprachsymbolen. 


Marken (Labels) 


Ein- und Ausgabe I | Unterprogramme 
— 


allgemeine Operationen 


NE 


Start/Stop Abarbeitungsrichtung 
. 5 ätzli 
bedingte Springe  ----- il zusätz iche Bemerkungen 
zu einer Operation 


Bild 1.1 


Das nachfolgende, sehr einfache Beispiel, zeigt die Schritte zur Programmentwicklung exemplarisch. 


Gesucht ist ein Programm zum Wurzelziehen aus beliebig reellen Zahlen. Kurz 


fix) av, xEeir 


Die Problemanalyse liefert einen Gültigkeitsbereich der Funktion für alle positiv reellen Zahlen. Ein 
möglicher Algorithmus ist damit: 


. Lies x ein 

. Ist x <0, dann weiter bei 6. 
. Bilde y=v/x 

. Gib y aus 

Stop 

. ‚Fehlermeldung‘ 

. Stop 


sıaon2un. 


Es gäbe auch z.B. noch die Möglichkeit, nach der Fehlermeldung Ix| zu bilden und nach 3. zu 
gehen, etc. 

Das Fiußdiagramm in Bild 1.2 macht den Algorithmus noch anschaulicher. Der für den T158/69 
modifizierte Algorithmus, also das Programm, hat dann die in Tabelle 1.1 wiedergegebene Form. 
Bild 1.3 zeigt die Programmanwendung an zwei Beispielen. 


WURZELZIEHEN 
17. x 
4, 122105826 Tx 
Ausgabe 
Fehler! 
WURZELZTEHEN 
-17. x 
FEHLER!" 
Bild 1.3 
Bild 1.2 Anwendung des Wurzelprogramms 
Fiußdiagramm zum Problem Wurzelziehen an zwei Beispielen 


Mit Hilfe der in Bild 1.1 beschriebenen Symbole, läßt sich auch die Wirkung des Dsz-Befehls 
(decrement and skip on zero) in seiner hauptsächlichen Anwendung als Zähler erklären. Bild 1.4 
zeigt diese Anwendung im Flußdiagramm. 


Tabelle 1.1 Programm-Wurzelziehen 


Bedingte Abfrage, 


Start und Ausdruck Ausdruck — Berechnung und 

Wurzelziehen gestrichelte Linie — Ausdruck 

D00 76 LEL ü43 63 DP 075 

001 ii A [I Tun TE 

002 25 CLR 45 0202 Orr 

003 69 OF age 00 0 urs 

004 00 00 ar le = 73 

005 04 4 ads 00 cd os0 

006 03 3 43 0258 081 

007 04 4 cso 000 Ügs 

008 01 1 üsi 0208 083 

009 03 3 use 000 084 

010 05 5 (s2 02058 185 

ii 04 4 os4 0000 ÜSE . 

012 068 & 05553 OF us? E) 

013 01 1 ose 01 01 LEE 62 

014 07 7 DS? 52 OF 0893  (E 

015 6° IF SE 0202 au 95 

Die 01 159 89 OP 031 28 

ir 02 28 nel 02 03 a2 98 

018 07 7 el 83 DF a3 98 

019 04 4 gez 04 04 09431 

020 06 & DES 6a DPF 

021 02 2 üe&4 05 05 

022 04 4 DES 95 AV : 

023 01 1 6 53 OF Fehlermeldung bei 

024 07 7 Ger nd dc negativern Eingabe- 

025 02 2 Er - wert 

026 03 3 TE u Te Br) 195 TE LE 

027 69 OF Orc 63 OF De ae 

028 02 08 rl 04 U 097 63 OF 

023 0 1 DIE 00 

030 07 7 033 02 2 

031 03 3 100 01 1 

032 01 1 101 01 1 

03 00 u Eingabe und Aus- i0o2 07 >? 

034 00 0 druck 103 89 OF 

eu 0 072 91 RrS ee 

037° 00 ö Dr2 6@ OF 105 br e 

038 00 0 74 06 De won 

039 6? DF 108 u7 7 

040 03 03 103 Di 1 

041 6°? OF 116 0? 7 

042 05 05 mit a3 3 
112 05 5 
113 077 
1il4 02 2 
115 689 OP 
t1l6 02 08 
il’ 62 OF 
115 05 05 
112 28 AIW 
120 98 AIW 
‚el 28 AM 
122 95 HIV 
1232 91 Res 


Ein Zähler i wird auf die Anzahl n gesetzt, die eine Befehlsfolge 
A abgearbeitet werden soll. In der Befehlsfolge A läßt sich der 
Zähler auch direkt oder indirekt als Adresse verwenden. 

Das in allen Flußdiagrammen auftretende Zuweisungszeichen (:=) 
ist nicht mit dem mathematischen Gleichheitszeichen zu verwech- 
seln. Die in Bild 1.4 verwendete Zuweisung 


i:=171 


ist mathematisch wenig sinnvoll, bedeutet aber informatisch die 
Verminderung des Wertes i um 1. Damit ist die Funktion des 
Zuweisungszeichens beschrieben. Der rechts des Zeichens stehen- 
de Ausdruck wird als Ergebnis der links stehenden Variablen zu- 
gewiesen. 


Bild 1.4 Die Anwendung des 
Dsz-Befehls als Zähler für 
Programmschleifen 


1.2 AOS-Technik 


Auf dem derzeitigen Markt für programmierbare Taschenrechner unterscheiden wir zwei Systeme, 
die mit algebraischer und die mit umgekehrter polnischer Notation. Mit den Wirkungsweisen und 
Zusammenhängen habe ich mich in [10] umfassend auseinandergesetzt. Da diesem Buch der T158/59 
zugrunde liegt, will ich auf das algebraische Organisationssystem (AOS) noch etwas eingehen. Eine 
umfangreiche Beschreibung finden Sie in [1]. 

Jedem Computer, auch einem programmierbaren Taschenrechner, liegt die in Bild 1.5 dargestellte 
Grundstruktur zugrunde. 


JL 


Bild 1.5 
Grundstruktur eines Computers 
(> Daten, — -» Steuerbefehle) 


Die darin enthaltene Speichereinheit läßt sich, wie in Bild 1.6 dargestellt, in 3 logische Einheiten 
aufteilen. Die Grenze zwischen DSE und PSE ist beim T1/58/59 variabel und man spricht von einer 
sogenannten dynamischen Speicherverwaltung. Die ASE hat bei der AOS-Technik einen Haupt- 
speicher, den sogenannten Akkumulator. Außerdem bedient sich die ASE einer Vielzahl von Hilfs- 
speichern, zur Speicherung von Zwischenergebnissen. Diese Speicherung ist, im Gegensatz zur um- 
gekehrten polnischen Notation, ohne Einfluß vom Benutzer und geschied nach den Gesetzen der 
Algebra. Den graphischen Zusammenhang gibt Bild 1.7 wieder. 


4 


Arbeits- 
Speicher- 
Einheit 


Daten- 
Speicher- 
Einheit 


DSE 


Programm- 
Speicher- 
Einheit 


Bild 1.6 


PSE Gliederung der Speichereinheit 


1. Hilfsspeicher 
2. Hilfsspeicher 
= 3. Hilfsspeicher 


Akkumulator 


Bild 1.7 


Arbeitsspeichereinheit (ASE) 
bei algebraischer Notation 


Interner Datenfluß ohne 
Einfluß vom Benutzer 


Wird nun ein geladenes Programm in der PSE aufgerufen, so veranlaßt dieses, über die Steuerein- 
heit, den Transport von Werten, aus der DSE, der PSE oder über die Eingabe, in den Akkumulator. 
Das Rechenwerk, veranlaßt durch die Steuereinheit, verknüpft die Werte aus Akkumulator und 
Hilfsspeichern zu sinnvollen Ergebnissen. Das Endergebnis steht nach Abschluß aller Operationen 
im Akkumulator und kann dann an die DSE oder zur Ausgabe weitergegeben werden. 


1.3 Allgemeine Programmiergrundlagen 


Da dieses Buch als Anregung zum eigenen Programmieren gedacht ist, sollen an dieser Stelle helfende 
Grundlagen und Tips vermittelt werden. 

Um den Algorithmus eines anstehenden Problems zu entwickeln, empfiehlt sich die Methode der 
strukturierten Programmierung. Dabei wird durch eine schrittweise Untergliederung des Problems 
(top-down-design) und deren Einzellösungen, das Gesamtproblem gelöst. Bild 1.8 zeigt ein Problem 
P und dessen Untergliederung in 3 Teilprobleme. Diese unterteilen sich je nach Möglichkeit wieder 
in Teilproblerne, usw., bis sich die Lösung eines Problems anschaulich aus der Summe der Einzel- 
lösungen ergibt. 

Liegt ein umfassender Algorithmus in Flußdiagrammform vor, kann die Umsetzung zum Programm 
erfolgen. Dieser Vorgang dürfte bei Kenntnis des Rechners keine großen Schwierigkeiten bereiten. 
Nach der Programmerstellung empfiehlt sich jedoch eine kritische Betrachtung zwecks Optimierung 
des Programmschrittbedarfs. Auf folgende Punkte ist dabei zu achten: 


1. Wiederholt gleiche Programmschritte lassen sich zu einem Unterprogramm zusammenfassen. 
Dies ist jedoch nur bei der Ersparnis von mehreren Programmschritten sinnvoll, da das Programm 
an Übersichtlichkeit verliert. 

2. Wiederholt ähnliche Programmschritte mit wechselnden Daten lassen sich durch indirekte Pro- 
grammierung im Programmschrittumfang vereinfachen. Diese komplexe Programmiermethode 
erfordert allerdings ein gewisses Maß an Verständnis und Übersicht. 


Bild 1.8 Aufteilung eines Problems in Teilprobleme 


Man kann auch schon beim Aufstellen des Flußdiagramms auf gewisse Rechnereigenheiten Rück- 
sicht nehmen. Weiterhin sollten bestehende Programmrestriktionen und daraus resultierende Fehler 
notiert werden, um bei einer späteren Benutzung unnötige Fehler zu vermeiden. Dies gehört aber 
schon eigentlich ins nächste Kapitel. 


1.4 Dokumentation 


Um eine sinnvolle Programmpflege betreiben zu können, d.h. ständig neu gewonnene Erkenntnisse 
in alten Programmen zu verwerten, bedarf es einer gründlichen Programmdokumentation. Es 
empfiehlt sich also die Benutzung von Formblättern. Zu einer umfangreichen Programmdokumen- 
tation gehören, in der Reihenfolge ihres Entstehens: 


oo PONn— 


. Berechnungsgrundlagen 

. Flußdiagramm 

. Speicherplatzbelegung 

. Programmplatzbelegung 

. Programmbeschreibung (Eingaben, Ausgaben, etc.) 
. Testbeispiel 

7. 


Angabe der Restriktionen und möglicher Fehler 


In einer solchen Programmbibliothek lassen sich dann auch bei späteren Problemen Teillösungen 
finden, Schreiben Sie dazu Ihre Programme abschnittweise, wie ich es in diesem Buch ebenfalls 
getan habe, auf und sparen Sie nicht mit Kommentaren. 


2 Statik starrer Körper 


Eine Kraft kann nicht unmittelbar, sondern nur anhand ihrer Wirkung beobachtet werden. Ihr 
Wirken zeigt sich in der Verformung eines Körpers oder der Änderung seines Bewegungszustandes. 
Aufgrund der fundamentalen Bedeutung der Kraft für die Mechanik, sollen im ersten Teil dieses 
Kapitels ein paar grundlegende Anschauungen und Programme behandelt werden. 


2.1 Kraft, Moment und Gleichgewicht 


Die Kraft hat, für ihre mathematische Behandlung am starren Körper, den Charakter eines linien- 
flüchtigen Vektors. Sie unterliegt damit den Gesetzen der Vektoralgebra und ist durch Größe, 
Wirkrichtung und Angriffspunkt eindeutig bestimmt. 


2.1.1 Reduktion einer räumlichen Kräftegruppe 


Dieser allgemeinste Fall einer Kräftereduktion beinhaltet alle möglichen Besonderheiten, wie später 
noch nachfolgend dargestellt wird. 

Wir betrachten die Anordnung von n Kräften F;; j= 1... n; bezüglich eines beliebig rechtwinkligen 
Koordinatensystems (e} , ea, e3), nach Bild 2.1, mit ihren Angriffspunkten a, an einem imaginären 
starren Körper. Die Reduktion dieser Kräfte bezüglich des frei gewählten Ursprungs u ergibt die 
aus den Vektoren 


F=yF, (2.1.1) 


und 


n 
M= %,Mu= Dax, (2.1.2) e3 
i=1 i=1 


resultierende Kraft und Moment. 


Bild 2.1 
Kräfte im Raum 


Die für ein Programm wichtige Komponentendarstellung ergibt sich unter Einführung eines freien 
Zählers i und zweier davon abhängiger Zähler p(i) und q (i} mit der Zuordnung 


i pP q 
1 2 3 
2 3 1 
3 1 2 


aus den Gleichungen 


A 
Fra N A. (e1,23 (2.1.3) 
i=t 
und 
n n 
Mu= \ My = N (&jp Fia a Fiol, i=1,2,3 (2.1.4) 


1 i-1 


Diese Schreibweise ist gezielt auf eine indirekte Programmierung abgestimmt. Auf eine Anwendung 
des Matrix-Programms des Standard Software Moduls habe ich hier verzichtet, um das Prinzip der 
indirekten Programmierung anschaulich zu demonstrieren. Das Matrix-Programm findet unter 2.1.2 
Verwendung. 

Der funktionale Zusammenhang der Zähler i, p,q kann in einer Ringanordnung von 3 Speichern 
und dem darin Herumschieben der Werte 1, 2, 3, wie in Bild 2.2 dargestellt, erreicht werden. 


i 
KERSER 
» ENREEREBEE NT 


Ringspeicher und Zähler 


Die betragsmäßigen Größen von resultierender Kraft und resultierendem Moment ergeben sich aus 
den berechneten Komponenten nach dem pythagoräischen Ansatz 


= 
a; 
i=1 


IF, y 
u 
IM‚I= 1/7 S Ma; (2.1.6) 


de 


(2.1.5) 


und 


Bild 2.3 
Darstellung der Dyname 


Wird die resultierende Kraft F, senkrecht aus der Ebene, die F, und M, aufspannen, so um den 
Vektor & verschoben, daß resultierende Kraft und Moment die gleiche Richtung haben, so bezeich- 
net man dieses Vektorpaar (F,, M,) als Dyname (Kraftschraube). Der Ortsvektor a5 des Fußpunktes 
f bezüglich des Ursprunges u (Bild 2.3) ergibt sich in seiner Komponentendarstellung aus 

1 


ale IF,12 (Fp Mug — Fra Mup), i=1,2,3 m) 


Mit Hilfe des Parameters p 


F,M, 1 % 
5 = \F,M 1 
Pi LE) 
läßt sich der auf den Fußpunkt f (Bild 2.3) bezogene Momentenvektor M; berechnen 
3 
T 
M=pF,=p a (2.1.9) 
ie 1 


Außer der Möglichkeit der Komponentenangabe, gibt es zur Richtungsangabe eines Vektors zum 
Ursprung und den Koordinaten, die Angabe der Richtungswinkel. Unter Betrachtung von Bild 2.4 
ergeben sie sich aus den Gleichungen 


a; = arccos ın|. i=1,2,3 (2.1.10) 
e3 
az a 
a, Bild 2.4 
s Richtungswinkel eines Vektors 
1 
% e2 


Einl i 
LBLA /Eingabeprogr. Dre " ae En 
teilweiser Auswertung 


Fr Eingabe 7 
/ 
eu 1=1,23,3 f 
Fi. ie1,2,3 Ä 
EF)sSr,+F, + --- -{ Summenbildung der Kraftkomponenten 
- : ge Summenbildung der 
EM.=EM, - ---- 
E Mui = Mi + (ap Fa "3a zu: 4 Momentkomponenten 


Ben Auswertungsprogramm zur 
e LBLB /Auswertung F---- Bestimmung der resultierenden Kraft 


Er und des Moments 
| 2. | 
| IF,| = | N“ Fu Ei mn 4 Betrag der resultierenden Kraft 
[EEBER RER. Ber 
ui 


u... 1 Betrag des Moments 


Keane { Ortsvektor zur Dyname 


2 
I 
n 
} 
ce 
: 
1 
A) 
3 
= 
[7 
5 
u 


3 
M = p S F ...- 4 Moment der Dyname 


/ Ausgabe 7 


& = zus r Hilfsprogramm zur Berechnung 
7 Eingabe 7 "der Richtungswinkel von 


/ 
K3 j n=1: M n=2:% 


n=3: M, n=4: F, 


& 
a, = arccos —, 
lal 


Ausgabe 
a. i=1,2,3 


( Stop/Eingabeende ) Bild 2.5 Flußdiagramm zur Kräftereduktion 
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Die Programmanalyse ist damit abgeschlossen und wir kommen zur Darstellung des Algorithmus 
in Flußdiagrammform. Er ergibt sich zwangsläufig aus der Abfolge der aufgestellten Gleichungen. 
Bild 2.5 ist eine mögliche Form. 

Zu Anschauungszwecken ist dieses erste Flußdiagramm etwas aufwendiger kommentiert als die 
nachfolgenden Diagramme. Das auf den T158/59 zugeschnittene Programm lautet, bei einer Spei- 
cherplatzbelegung nach Tabelle 2.1, wie unter Tabelle 2.2 wiedergegeben. Eine Kommentierung 
der Ein- und Ausgabedaten durch das Programm, habe ich aus Übersichtlichkeit unterlassen. Ta- 
belle 2.2 gibt also das ‚nackte‘ Programm wieder. Dies gilt gleichfalls für alle nachfolgenden Pro- 
gramme. 


Tabelle 2.1 Speicherplatzbelegung zum Programm Kräftereduktion 


00 Zähler und Ind. Adr. 


01 Fa pFra = M4g 10% Fa =, Fra 
02 Fa} F, pFa=Ma !M, 11 ZFa=Fo gF, 
03 Fu pFri = My, 12 & Ri = Fra | 
04 eı3 a3 13 IM;| 
05 ei2 e | a2) % 14 la,| 
06 e1 Er27 15 IM, | 
07° EMys= Ma 16 IF,| 
08 ZMy2= Ma |\M, 17 p 
09 2} My = Mu1 18 
19 } Zähler + Ind. Adr. 
20 
21 24 
22 } 1. Ringzähler 25 ?}2. Ringzähler 
23 26 
” 
113 2 416 5 
1.Fal | 2) 1 3 | | 514 6 
3) 2 1 6,5 4 
[ib 
‚19 8 10) 12 11 
2. Fall 81 ı 9 | | 11 10 12 | 
9 8 7 12 11 10 


Eine ebene Kräftegruppe läßt sich mit diesem Programm ebenfalls behandeln. Dabei wird lediglich 
die 3. Komponente der Vektoren Null. 

Haben alle Kräfte einen gemeinsamen Angriffspunkt, so ist es sinnvoll, den Koordinatenursprung 
in diesen zu legen. Die Koordinaten des Angriffspunktes werden dann mit Null eingegeben. Ein 
Moment tritt für dieses System nicht auf, erscheint also mit Null als Ausgabewert. Sind resultieren- 
de Kraft und Moment Null, gilt also 


F},=0 und M,;=0, füralle i=1,2,3 (2.1.11) 


so befindet sich der, diesen äußeren Kräften ausgesetzte, starre Körper im Gleichgewicht. 
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Tabelle 2.2 Programm Kräftereduktion zum T158/59 


Eingabeprogramm: 


Berechnung zur 


Start resultierenden Kraft Ringshiften 
ann TE LEL ülE 21 Re 31 10 E° 
ol 11H 0l2 23 FET 
au2 47 CME alt TE 5T# Zählerverminderung 
l2 08 & 015 co 00 + Rücksprung 
coas4 18 C*" oile Ar DEE 32 17 B' 
üir ac oc 033 97 IS2 
Eingabe le 79 > 054 00 08 
cos TE LEL 43 92 Aw 033 65 & 
DCE 85 + ü3E 61 GTO 
nn IE & Berechnungsvorbe- 027 ss + 
u: u = reitung 
og 42 5T0O or A 
one ca dr Fe 7 
010 75 LEL ee; 
oil 79 - 
Auswertungsprogramm: 
Bestimmung des 
Start Ortsvektors/2. Teil 
028 TE LEL us2 10 E' 080 423 RCL 
039 12 E 0S2 178" 081 16 16 
ü4dü 935 hly üs4 47 DEE 082 35 178 
Üdt 95 AD“ 055 000 00 083 49 PRD 
ü4E 01 1 ÜS6 89 + 084 04 04 
D43 025 (85 49 FRD 
tt 13° Bestimmung der 086 05 05 
Vektorbeträge/1. Teil 08? 49 PED 
Bestimmung des Orts- sr 88 06 06 
vektors/1. Teil uSE 083 42 STO 
45 0c 6 059 090 17 17 
o4e IE A" DEU 091 33 88 
47 TE LEL 061 09249 PRD 
4 55 + Ü&E 093 14 14 
4a 193 I° 063 
SO TE S5T# ed 
051 18 15 0865 
Bestimmung des 
Parameters 
(94 093 11lü 49 FED 12 IuB 1328 75 LBL 
0939 18 111 1° 17 1 Er 139 42 STD 
ÜBE TE . 1 44 N 140 73 RC# 
097 SE Bestimmung des i 13 13 141 18 18 
08 7 Momentes’M 1: Ir e' 14234 TH 
099 18 118 12 “7 132 143 72 5ST* 
100 65 = 113 1 ca 0 144 18 183 
101 TR Rus 114 1: 44 ZUM 145 178° 
102 197 197 115 146 97 DSZ 
1U2 85 + 116 Bestimmung der 14’ 00 00 
104 ir E*° 117 Vektorbeträge/2. Teil 145 42 STO 
105 #7 ISZ 118 ci 1 
10e 00 00 1193 a6 Ausgabe + Ende 
107 52 EE 120 2 15 A! 1932 00 0 
102 00 0 121 r I | 150 22 INW 
109 95 = 122539 44 5LUM 151 90 LST 
IE. 2 on on 152 91 RrS 
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Unterprogramme: 


Setzen der Zähler 


153 76 LEL 
154 16 RA' 
155 42 STD 
156 18 18 
157 85 + 
158 03 3 
159 85 + 
160 42 STO 
161 19 19 
162 03 3 
163 95 = 
164 42 STD 
165 20 20 
166 03 3 
167 42 570 
168 00 00 
169 92 RTN 
Verminderung der 
Zähler um 1 

170 7& LEL 
171 17 B" 
172 01 1 
173 94 +7- 
174 44 SUM 
Hilfsprogramm: 
Start 

237 76 LEL 
238 13 C 
239 9% PET 
240 98 AIW 
241 65 # 
242 03 53 
243 42 STD 
244 00 (Wü 
245 95 = 
246 42 STD 
247 18 :8 
248 02 2 
249 09 9 
250 42 5”0D 
251 19 19 
252 00 0 


253 42 570 
254 20 7020 
255 76 LEL 
256 34 TE 


175 18 158 
176 44 SLM 
ir7 19 19 
ira 44 ZUM 
1793 z0 Zu 
150 92 RTH 


Setzen der Ring- 
zähler 
181 
1e2 
183 
154 
185 


- 
[nt 
| 
Oo& 

rare 

ie} 
Dem} 


r£oo 


pr 


+ (0 r 


zär 72 
ESE 15 
2S9 TE 
a u 5 4 
zei 35 
2e2 44 SUM 
geE2 oO 2O Eu 


Zählerverminderung 
um 1 + ENERSREUNG 


127 19% 19 
138 178" 
139 97 DEZ 
zoo 00 6 
zul 38 


Berechnung mit 
Ringzähler 


Tinkicfı 


Ffotsrofofsfatarfatofafotarafafe 


mm. 
7 


Berechnungsvorbe- 


reitung 


u 
A 
ri 
ni 
m 


rafıı BAR af 


Ausgabe 


Ringzählervertau- 
hung 
1 pe} 


[= 
ziS 
- 


ats 


250 43 RüL 
281 20 20 
EEE 
293 22 IMY 
4 39 008 
235 93 PET 


Zählerverminderung 
um 1 + Rücksprung 
ol. 

34 +. 
44 SUM 
15 
ar 
win] 
45 
ol 
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2.1.2 Zerlegung einer Kraft 


Der Vorgang der Kräftereduktion ist ein umkehrbarer Prozeß. Aus diesem Grunde läßt sich eine 
Kraft wieder in Komponenten zerlegen. Eine Zerlegung ist eindeutig möglich, wenn 3 unabhängige 
Richtungen für die Komponenten gegeben sind. Es soll somit die Kraft F in die Komponenten Fj, 
j= 1,2, 3, zerlegt werden. Diese Komponenten bilden bezüglich des Koordinatensystems die 
Summen 


3 e3 
F= 2, Fi: (2.1.12) 
i=1 
Bild 2.6 
Kraftkomponenten 02 


e > .. Fa Es 


Unter Benutzung der Richtungswinkel (Bild 2.4) a;; ergibt 


3 
Fr cosa;=F, i=1,2,3, (2.1.13) 
i=1 


alle Komponenten längs der Koordinatenachsen. In abgekürzter Schreibweise für 


Cosa, = 2; (2.1.14) 
wird damit 
3 
ba Fai=F, i=1,23. (2.1.15) 


ie1 


Wie oben angedeutet, ist dieses Gleichungssystem unter bestimmten Bedingungen eindeutig lösbar, 
denn wir erhalten 3 Gleichungen mit 3 Unbekannten. 
Zur Lösung dieses linearen Gleichungssysterns 


AF=F (2.1.16) 


bietet sich das Solid-State-Softwareprogramm ML-O2 an. Dabei entsteht eine Schwierigkeit, denn 
durch die Verwendung als Unterprogramm ist die Eingabe über R/S nicht möglich, da das Unter- 
programm durch INV SBR ins Hauptprogramm zurückspringt. Listet man sich diese Programmteile 
aus, so ergibt sich für unser Problem folgende, durch ML-O2 vorbestimmte und in Tabelle 2.3 wieder- 
gegebene, Speicherplatzbelegung. Unter Festlegung des Berechnungsalgorithmus in Flußdiagramm- 
form (Bild 2.7), ergibt sich das entsprechende Programm in Tabelle 2.4. 
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LBLA/Start 


Programm darf nur dann mit RIS 

_-—--—-—- - weiterlaufen, wenn IAl #0 
Nur dann sind 3 unabhängige 
Richtungen gegeben! 


Ausgabe 


Fi, j = 1, 2, 3 . N 
Bild 2.7 Fiußdiagramm zur Kraftzerlegung 


R/S/Ende 


Tabelle 2.3 Speicherplatzbelegung 


vorbestimmt: 

08 211 11 951 14 431 21 Fı 

09 412 44j 12 822 47; 15 932 a3i 20 Fa F; 
10 413 13 923 16 933 22 F3 

frei gewählt: 

00 05 F, | 

01 7 Zähler 06 Fz, | Fi 

02 07 Fa 
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Tabelle 2.4 Programm Kraftzerlegung zum TI 58/59 
Determinantenbe- 


Start/Einlesen rechnung und Ausgabe 
(ST 22 INY 083 03 3 
055 &4 FI* 084 42 STD 
059 £ ü85 07 07 
[u T=20] 086 36 PGM 
sl üsr 02 02 
[uT-0e 088 13 C 
Te) 089 91 RS 
084 = 
u 5 = E ? Lösungsberechnung 
1 Sa Di (nur möglich wenn 
a en Determinante # 0) 
Ges 00 090 25 CLR 

092 02 02 
093 15 E 

Umspeicherung 

Berechnung - 4a Bm 

Richtungskosinusse Ei BL Ag ar 

15 TE cEL In ab a 

die 12 6 ur Bas 20 ‘56 
a 096 99 PRT 
a 097 43 RCL 
On 098 21 21 
3 099 99 PRT 
den 100 43 RCL 
0 101 22 22 
8: 102 99 PRT 
oe 103 91 RS 


Das Programm kann nach Ausgabe der Determinante nur dann mit R/S weiterlaufen, wenn diese 
ungleich Null ist. 


2.1.3 Stützkräfte in Tragwerken 


Wir kommen zur Anwendung der Gleichgewichtsbedingung. Das Ziel der statischen Auslegung von 
Tragwerken ist es, daß diese bei Einwirkung äußerer Kräfte ihre vorgesehene Ruhelage beibehalten. 
Die Belastungs- und Stützkräfte müssen sich also im Gleichgewicht befinden. Die Untersuchung 
darüber kann natürlich mit dem zuvor bestimmten Programm geschehen. 

Da aber ein Großteil von Berechnungen sich auf ebene Fachwerke beschränkt, soll nachfolgend ein 
Programm zur Berechnung ebener Fachwerke nach dem Knotenpunktverfahren entwickelt werden. 
In diesem Programm kann dann auf spezielle Belange des Berechnungsverfahrens eingegangen wer- 
den. 

Unter einem ebenen Fachwerk versteht man idealisiert ein Gebilde aus geraden Stäben, die in ihren 
Endpunkten (Knoten) durch reibungsfreie Gelenke miteinander verbunden sind. Die äußeren Be- 
lastungskräfte greifen dabei nur in den Knoten an. Durch diese Idealisierung können in den Stäben 
nur Zug- oder Druckkräfte übertragen werden. Insbesondere gilt damit für jeden Knoten des Fach- 
werks die Gleichgewichtsbedingung. Betrachten wir einen Knoten bezüglich eines Koordinaten- 
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Bild 2.8 
Ebenes Fachwerk 


Bild 2.9 
Kräfte im Knoten 


systems zur Festlegung der Kraftangriffsrichtung, denn die Kräfte können nur längs der Stabrich- 
tung wirken, so muß die Gleichgewichtsbedingung 


n 
Nr =0 (2.1.17) 
i=1 
erfüllt sein. In Komponentenschreibweise heißt dies 
n n 
D,Fieosa;=0 und IF; sina; = 0. (2.1.18) 
i=1 i=1 


Damit liegt ein lineares Gleichungssystem von zwei Gleichungen vor. Folglich lassen sich damit 
zwei unbekannte Stabkräfte bestimmen. 
Seien also F,, i= 1,2 die unbekannten und F,, j = 1,..., n bekannte Stabkräfte, so ergibt sich 


n 
F, cosa; + Fa cosaa = - > F, cosa; 
1 
(2.1.19) 
Fi sina, + Fz sinaz =- ) Fjsing. 
In Matrixschreibweise 
AF=F (2.1.20) 


n T 
mit F=(F,, F,)" und F’= (- 5 F,cosa;, $ Fr ‚sina,) ; sowie der Matrix 
ji=1 i=1 


# ( cosa, cos =) 
sina, sindo 


Die Lösung läßt sich wiederum sehr gut mit dem Solid-State-Softwareprogramm ML-02 bestimmen. 
Bild 2.10 zeigt das Flußdiagramm und nach Tabelle 2.5, der Speicherplatzbelegung, folgt in Ta- 
belle 2.6 das Programm. 

Die Eingabe der Wirkrichtung der unbekannten Kräfte sowie die Eingabe der bekannten Kräfte und 
ihrer Wirkrichtungen muß nicht nach einem bestimmten Umfahrsinn, ähnlich dem zeichnerischen 
Cremonaplan-Verfahren, erfolgen. Die Reihenfolge der Eingabe der Wirkrichtungen stimmt mit der 
Reihenfolge der Ausgabe berechneter Stabkräfte überein. Die richtige Behandlung wird im Anwen- 
dungsbeispiel gezeigt. 
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LBLB /Berechnung 


LBLA /Start 


Eingabe 


a, i=1,2 


cosa;, sing; 
(P>R) 


Ausgabe 
Fj 
i=1,2 


R/S /Ende 


Dieser Programmteil wird solange durchlaufen, bis alle 
--- bekannten Kräfte an einem Knoten mit ihren Richtungen 
eingegeben sind. Danach weiter mit B 


Bild 2.10 Fiußdiagramm zu Stützkräfte in Tragwerken 


Tabelle 2.5 Speicherplatzbelegung 


00 Zähler 08 cosaı 14 —- F,sina, | F, 

01 F; 09 sina, C 

07 2 10 c0saz 15- F,cosg |Fz 
11 sinag 
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Tabelle 2.6 Programm Stützkräfte in Tragwerken 


Start 
000 TE LEL »T Eingabe bekannter 
001 1198 Ep Kräfte + Winkel "4 
002 47 CME To 42 13 
003 02 2 11 42 te n° 
004 42 5TO T ud 
005 07 rn D 45 Lösungsbestimmung 
io de 5 

Winkeleingabe und i E 
Koeffizientenbe- Berechnung der N 45 
rechnung Determinanten lid a 
00& ci 1 MET E12 osc 
007° 32 41T IS] 
008 00 0 
009 91 RS 
010 9@ PRT re 
011 37 PR a 
012 42 STO 14 14 
013 03 09 nt 35 FET 
D14 32 X1T n +2 RÜL 
015 42 STD 7 15035 
016 08 08 & 34 FRT 
017? 0 1 44! il Re2 
018 32 81T LE Li 
019 S1 R“S 44 SM [ 

12 12 in 


2.1.4 Biegeträger 


Die idealisierte Vorstellung, daß äußere Kräfte nur in den Knotenpunkten eines Fachwerkes an- 
greifen, trifft in der Realität nicht zu. Durch auftretende Flächen- und Punktlasten außerhalb der 
Auflager werden Stäbe auch auf Biegung beansprucht. Natürlich gibt es auch einfache Biegeträger. 


Das nachfolgende Programm soll die Möglichkeit bieten, neben einem Rumpfprogramm z.B. die 
Gleichung der elastischen Linie eines beliebigen Belastungsfalls einlesen zu können. 
Dies setzt die Aufteilung 


Magnetkartenseite 1: Spezialprogramme 
Magnetkartenseite 2: Rumpfprogramme 


voraus. Gehen wir von der Annahme aus, daß die Spezialgleichungen als Hauptprogramme abge- 
faßt werden, so wird das Rumpfprogramm als Unterprogramm ausgebildet. Eine mögIche Zuordnung 
wäre: 


LBL E: Eingabe 

LBL A: Gleichung der Balkendurchbiegung 
LBLB: Gleichung der Balkenneigung 
LBLC: Gleichung des Querkraftverlaufs 
LBLD: Gleichung des Momentenverlaufs 
usw. 
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Betrachten wir als konkretes Beispiel einen nach Bild 2.11 einseitig eingespannten Träger, durch 
Punkt- und Streckenlast belastet. Das Moment an einer beliebigen Stelle x ergibt sich nach Bild 2.11 
aus dem Ansatz i 


Mi = FI-n tal %. (2.1.21) 


Bild 2.11 Einseitig eingespannter Träger mit Punkt- und Streckenlast 


Durch Einsetzen in die allgemeine Gleichung der elastischen Linie 


„__Mi 
vn (2.1.22) 
ergibt sich 
„__ıL - | _.u2 
y. gr + x)°). (2.1.23) 
Eine Integration führt auf die Gleichung der Balkenneigung 
u ee Fa I a Ta 
Va-c (Fix ee +). (2.1.24) 
Mit der Randbedingung y’(x = 0) = 0 folgt cı=0. 
Eine nochmalige Integration führt auf die Gleichung der Durchbiegung 
1 (F\2_F a,922_9, 3,9 a ) 
=-— [n2-xn+ 1222 -< = £ 
y er 51% 6“ +71 gr er +Ca (2.1.25) 


Mit der Randbedingung y(x =0) = 0 folgt ca = 0. 


Damit erhalten wir die 3 Spezialgleichungen: 


= . (ar (3 (2) - (ZY)+a (6 (#)-4(2) + ))). (2.1.26) 
ve) a1 


< 


und 


mo 3 (1-*)(2r+a (1-%)) (2.1.28) 
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in denen die Laufvariable x nur im Quotienten x/l auftritt. Bezeichnen wir diese Gleichungen allge- 
mein mit f (x), so zeigt Bild 2.12 das Flußdiagramm zur Berechnung von f(x) über den ganzen 
Träger mit der Schrittweite Ax. Tabelle 2.7 zeigt die Speicherplatzbelegung mit Einteilung fester 
und freier Belegung. Die feste Belegung ist für ein einwandfreies Funktionieren des Rumpfpro- 
gramms notwendig. Auf die freie Belegung greifen nur die Spezialprogramme zurück. Das entspre- 
chende Programm zu diesem allgemeinen Anwendungsfall zeigt Tabelle 2.8. 


LBL__/Berechnung 


Ausgabe: 


x“ 
fix;) 


LBL EE / Unterprogramm 


ja 
X INVSBR /Rücksprung 


Tabelle 2.7 Speicherplatzbelegung 


LBL E/Eingabe 


Eingabe 
(feste Belegung) 
L 
Ax 


Eingabe 
(freie Belegung) 
F,O,EI 


R/S / Eingabeende 


R/S /Ende 


Bild 2.12 Flußdiagramm zum Biegeträger 


feste Belegung: freie Belegung: 
00 Zähler 05 F 
0 x/ 06 Q 
02 x 07 E 
03 1 08 I 


04 Ax (Schrittweite) 
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Tabelle 2.8 Programm Biegeträger zum T158/59 


Eingabeprogramm 
Gca TE STE Ei 65 x 
un 4 43 RCL 
Ss 


RTH 


na Te 


III EI TI TI IR 
[SL Ed ARTEN EN 
[0 
nn 
x 


a 


ne v 65 x 

B 43 RÜL 
al ne 03 03 
se en 1793 55 - 
15 Ice 1e0 02 2 
4 [122 1e1 95 = 
ge an 1E2 99 PRT 
51 ar 183 681 G”D 

- 184 12 6 

73 

ur Unterprogramm 

075 Rumpf 

5 240 TE LBL 

078 241 52 EE 

“17 SS 242 43 RÜL 

a 243 02 02 

Balkenneigung 24 Tr 244 5. 

en: 245 43 RÜL 

[ra TE LEL i: EEK 245 04 04 

cEü 12 &8 133 01 4 Sr 5: 

det a A 1go na 58 248 42 870 

ugz Fl SER 14ı 000 

ü&3 52 EE 142 55 = 

DE4 DENT 14389 0 

JES 93 PRT 144 95 = 

NEE 55 7 145 93 PRT 

Der SENT 146 61 570 

DEE 935 2 147 12 8 

E93 0 4E SO 

sc 01 ci 

o21 02.5 Momentenverlauf 

En BR 148  LEL 

1! 42 FÜL 143 Ü 

Wi 


nd 
z 


Ü 1 
iu ES = 1 af 
90 067 02 2 18 i 
dl Ü BSH 18 = 
od & 053 42 FÜL 1 r = 
54 s 1c0 ui [1 1E8 = Unterprogramm 
as n 5 1 : on 1 Sr ? = Potenzieren 
D46 163 10 E' Er SE oe 
| Ie4 54: 1e0 01. 268 32 MT 
ICH 85 7 Let 793, 28 269 53 
en 1e2 42 5°0 270 43 RCL 
ae 1 a m 271 01 0 
Me IEd 65 = ar2 45 yx 
I sn 185 43 RÜL 273 32 X-T 
ill 03.23 1ER ODE U zur 54 
111 6908 ler 890 - 275 92 RYN 
112 42 RÜL les 025 


2.1.5 Anwendungsbeispiele 


Die jedem Abschnitt folgenden Anwendungsbeispiele sollen Ihnen die Verwendung der zuvor ent- 
wickelten Programme demonstrieren. Mit einem Beispiel lassen sich jedoch nicht alle notwendigen 
Informationen vermitteln. Es liegt an Ihrem Verständnis der Programme, inwieweit Sie diese ein- 
setzen. Anwendungsbeispiele demonstrieren auf anschauliche Weise, welche Programmabläufe 
sinnvoll oder umständlich sind. Änderungen gehören zur allgemeinen Programmpflege. Ergänzungen 
und Kommentierung der Programme können, je nach Rechnertyp, hinzukommen. 


BER 
An einem starren Körper greifen an den angegebenen Punkten folgende Kräfte an: 


a; = (3,0,0); Fı = (0, 3,0) 
a2 = (0, 1,0); F> = (0,0,4) 
a3 = (0,0, 2); F3 = (5,0,0). 


Gesucht ist ihre Wirkung auf den Körper. 


Eingabe (Aufruf mit A): 


Sinn 


t |» - 
5 2 le 


Die Eingabe Angriffspunkt/Kraft muß in dieser Reihenfolge geschehen. Die Reihenfolge der Index- 
werte ist belanglos. Nach jeder Eingabe beginnt das Eingabeprogramm erneut. Sind alle Daten ein- 
gegeben, wird mit B das Berechnungsprogramm aufgerufen. 


Ausgabe: 

Nach Berechnungsende wird durch 1 INV 2nd List der Inhalt der Datenspeicher, beginnend mit 01, 
ausgelistet. Nach Ausgabe von 17 stoppen Sie mit R/S die Ausgabe ab. Die Zuordnung Datenwert / 
Datenbezeichnung ergibt sich aus Tabelle 2.1. 

Die Bilder 2.13 und 2.14 geben den Sachverhalt dieses Beispiels bildlich wieder. 


Giüc al) 
zueccde 


Lew} 


6.88 01 Ma 
5.16 DE My2 M+ 
8.6 03 Myı 
0. r& 04 a3 
-0,58 08 a2 EI; 
-0.26 Üb a1 
9 7 Musa 
10 [jal= Mu2 Mu 
4. 09 Muı 
4. 10 Fra 
3. 11 Fr2 Fr 
5. 12 Fri 
12. 16223664 13 IMEI 
. J907Sr7286l 14 lag! 
14. 03566885 15 IMU! 
7. 071067812 18 IF, 
1.72 17 p 
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1 
Maßstäbe: | Weg, Kraft 


2 
— Moment 
äbe: 1/4 
Bild 2.13 Maßstäbe: J Weg 
Resulitierende Kraft und Moment 2 PIREBER 
1 Kraft 
— 


Bild 2.14 Darstellung der Dyname 


Richtungswinkel (Aufruf mit C): 


1. M Zu &% 2. Mu 4. Fr 
F2.4418 ay . 45. 
44,5 a2 64. 29590975 ag 


50. 11885 


ag 79. 550069801 a3 


Hier haben Kraft und Moment die gleiche Richtung. 
Das ihre Richtungswinkel gleich sind, beweist Hilfs- 
programm C. Durch Eingabe einer Kennziffer, siehe 
Bild 2.5 und Aufruf C wird die Berechnung der 
Richtungswinkel veranlaßt. 


AR PER 


Eine Kraft F = (2.25, - 7.5, - 2.25) soll in drei Komponenten zerlegt werden. Deren Richtungen 
sind durch die Vektoren 


a4 = (4,-5, -0.5) 
a, = (- 3.5, -5, 2.5) 
a3 = (0.5, -5, -5.5) 


vorgegeben. 
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Eingabe: 


nınen 


[ 
De 


' 
in“ 


Die Reihenfolge der Vektoreingabe ist mit der Reihenfolge 
der Komponentenausgabe identisch. 


natur m 
nen inne ® 


ı Zt 


_Nn 
in® 
ne 
— 
[] 
o 


ae 


Nach Beendigung der Eingabe wird durch den Aufruf des Solid-State-Softwareprogramms ML-02 
die Determinante des, für dieses Problem vorhandenen, Gleichungssystems berechnet und ausge- 
geben. 


Nur wenn dieser Wert # Null ist, darf das Programm mit R/S 


Fı wieder gestartet werden. Es erfolgt dann die Ausgabe der 
Fa Komponenten. 
Fa 


Bild 2.15 
Kraftzerlegung 


a3 
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Hinweis: 
Wird das Hilfsprogramm LBLC der Kräftereduktion nach der Berechnung eingelesen, so lassen sich 
mit der Zuordnung: 


n= 7=73-F, 

n = 10+*3-a,l[F}), 
n = 1373-3 (F,)}, 
n= 16+3-a3z(F,), 


die Richtungswinkel der einzelnen Vektoren bestimmen. 


a 
Das in Bild 2.16 dargestellte Tragwerk wird mit der eingezeichneten Kraft belastet. Gesucht sind die 
Stabkräfte. Die Auflagerkräfte ergeben sich aus dem äußeren Gleichgewichtszustand zu 


A+B+F=0 | A=- 5X 36=-09 
12XA+3XF=0 | B = -3.6+09=-2.7 


Bild 2.16 
Tragwerk 


ko 


-2,7 Bild 2.17 Größen und Bezeichnungen am Fachwerk 
3,6 


Für die Beschriftung der jeweiligen Knotenberechnung habe ich das in Tabelle 2.9 wiedergegebene 
Programm benutzt. Es berechnet für die Zahlen n = 1 bis 99 zu jeder Ziffer z ihre Kennziffer k für 
den Drucker nach der Gleichung 


k=z+1+Intlz+ 7)X 2. (2.1.29) 


Die Beschriftung der n-ten Knotenberechnung erfolgt durch den Programmaufruf nE. 
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1. KNOTEN? 


zu be- 
rechnen- 
de Kräfte 


bekannte 
Größen 


2. KNOTEN?! 


15%. 4349488 
233. 1301024 


 BEZSEFILLI 


-1,.125 
306, 590976 


1. OTTFÜEDFE 
‚Silzessaed 


3. ENDTEH? 


112. 1955905 
150, 


0. 21914503 


0,675 

Ü. 
‚SI13E2E3E4 
53. 12010235 


.des0rre4e! 
1. 157142557 


1. Winkel 
2. Winkel 


Determinante 
# Null 


1. Kraft A 
1. Winkel 


1. Kraft Sı 
2. Kraft Sa 


4. ENDTEHE 


24, 2014055 


D.F12314902 


Bild 2.18 Der erste Knoten 
mit seinen Daten 


5. EHOTEHS 


= bHUOTEH? 
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Tabelle 2,9 Beschriftungsprogramm zur Knotenberechnung 


ve LBL Sie 03 03 
iS E si? 8 + 
47” CMS 215 04 4 
53 DF 21° 00 0 
oc Du a20 02 2 
42 5STD 221 06 6 
O1 01 222 023 
r& LEL 322 01 1 
43 REL 224 9% = 
Si m [! s25 69 DPF 
u D Eu: rs ze Cl 01 
[ui] =: S327 033 
35 4E 228 02 2 
a2 {1 228 03 3 
[ zz IN“ 330 v7 
HT 31 01 i 
232 077 
322033 
ol a334 01 1 
= 335 06 6 
ı 226 022 
2] 237 8% OP 
= 2235 02 02 
A 339 6% DP 
EET 40 05 05 
41 91 RS 


-4- 
Für den in Bild 2.19 dargestellten, einseitig eingespannten Träger mit Strecken- und Punktlast, sind 
Durchbiegung, Balkenneigung und Momentverteilung gesucht. 


Bild 2.19 


Einseitig eingespannter Träger 
mit Strecken- und Punktlast 

E = 2,1: 10° N/cm? 

I =6: 109 cm® 


Die Eingabe erfolgt durch den Aufruf E: 


Trägerlänge I! (cm) 
Schrittweite x (cm) 


Pan] Einzellast F (N) 
Solndl, Streckenlast Q (N) 
2100. E-Modul (N/cm?) 
EINTETNIHFR Axiales Widerstandsmoment (cm?) 
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Durchbiegung (Aufruf A): 


2. Stelle x (cm) 
-0. 01852351 Durchbiegung y (cm) 


ER 18, 34, 
-, corzüu4e328 -, [E25 555 -, 23T L0SE 
Er ie 
-. (157571423 - UAEHEZOIDE a 
2. 
-, VEFZZSTIER 2 


-, 0413359788 


-, [OYSZEST14 


Neigung (Grad) (Aufruf B): 


& Stelle x (cm) 
-. 1046451645 Neigung y’ (grd) 


4. 
-. 2005659435 -, 58 


Ein 1E. Ze 
. 2ZE811I9175 I ELZOSTLTE TEHISERÜHLL Best 


E 
. 3ETERZOTEE 


10, en 
-, 495707952 Berne a en 


-.S5ldael23737 


Momentverlauf (Aufruf C): 


2 Stelle x (cm) 24, 
1102000. Moment M (Ncm) recdo, 


4. 10. 18, 22, 


1008000. FSsooco, SES000, Z4H2ZOCL, 1 FEN, 


3 3E 
150000, 2000, 


8, 
18000. ) 


su, 
Ü, 


8. 14. zu, 2 
zzand 


832000. dl, ELLE, 


hi 
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Bild 2.20 
Durchbiegungs-, Neigungs- 
und Momentenverlauf 


- 1° 7-0,5cm 
viy 


2.2 Seiltheorie 


Anders wie im vorausgegangenen Abschnitt, geht man bei der mathematischen Formulierung der 
Theorie gespannter Seile davon aus, daß diese keine Biegemomente übertragen können, mit anderen 
Worten also biegeschlaff sind. Daraus resultiert, daß Querkräfte nicht und Zugkräfte nur tangential 
zur Seilkurve auftreten. Dies trifft in gleicher Weise auch auf Ketten zu. 


2.2.1 Seil unter Eigenlast 


Wir setzen ein Seil mit konstantem Querschnitt und homogener Massenverteilung voraus. Wir ideali- 
sieren das Seil weiterhin zu (unendlich) vielen und (unendlich) kleinen starren Stücken, die unter- 
einander durch reibungslose Gelenke verbunden sind. Diese können nur Zugkräfte tangential zur 
Seilachse übertragen. Die Betrachtung eines solchen infinitesimalen Seilelements liefert nach Bild 2.21 
folgenden Kraftansatz in vektorieller Schreibweise 


(+ ar x (8 + a8) Hat) + (r+S )xadl = 0. (2.2.1) 


Bild 2.21 
Seilelement 
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Da wir nachfolgend ausschließlich ebene Belastungszustände betrachten wollen, vereinfacht sich 
unser Ansatz zu der Betrachtung nach Bild 2.22 mit den Gleichungen 


H=H+dH (2.2.2) 


S+ds 


V+dV 
und 


Vradi=V+dV. (2.2.3) 


Bild 2.22 
Ebenes Seilelement 


Die Seilkraft S zerlegt sich in die Komponenten H und V. Aus der Gleichgewichtsbedingung folgt 
sofort, daß 


dH =0. (2.2.4) 
Das heißt, der Horizontalzug H ist längs des Seiles konstant. Weiterhin ist 

dV = qdl. (2.2.5) 
Das Neigungsverhalten des Seilelements an der Stelle (x, y) ergibt sich aus der Gleichung 

vr (2.2.6) 
Das Verhalten längs des Seiles bestimmt damit die Differentialgleichung 

2 

Ma (2.2.7) 
mit 

di=]/1+ (Dan (2.2.8) 

dx 

folgt weiterhin 

d?y q_/ (av 

Gay (5) ui, (2.2.9) 


2.2.2 Die Seilkurve als Variationsproblem 
Die Differentialgleichung (2.2.9) läßt sich näherungsweise durch eine Differenzengleichung der 
Form 


q 
Atana = vr (tan«)? : Ax (2.2.10) 


ersetzen. Sie ist die Grundgleichung unseres Berechnungsalgorithmus. Unter Vorgabe der Strecken- 
last q, des Horizontalzuges H und eines Wegelements Ax (je kleiner Ax, umso besser die Approxi- 
mation an die tatsächliche Kurve), läßt sich die Neigungsänderung eines Seilelements gegenüber dem 
Nachbarelement bestimmen. Auf diese Weise erhalten wir iterativ die Bestimmung aller Seilelemente. 
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Den Berechnungsalgorithmus in Flußdiagrammform zeigt Bild 2.23. Unter Festlegung der Datenre- 
gister nach Tabelle 2.10 folgt unter Tabelle 2.11 das Programm. 

Der Vorteil dieses Programms ist bei kleiner Schrittweite eine sehr genaue Annäherung an den tat- 
sächlichen Seilverlauf und die Möglichkeit, auch komplexere Vorgänge wie unter 2.2.5—3— gezeigt, 


berechnen zu können. 


LBL A /Start 


Eingabe: k 
@1,X1,Yı, 1, %2 
Ax 


Ausgabe 


1x, Yili 
ai 
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Druckersteuerung durch Faktor k 


k 


k 


RIS /Programmende 


=0: Ausgabe nur für den Punkt (xz, y2) 
mit Anzeige von x; beim Durchlauf 
= 1: Protokoll aller Daten 


Bild 2.23 Flußdiagramm zur Seilkurve 


Tabelle 2.10 Speicherplatzbelegung 


O0lgq 
02 H 
03 Ax 


10 tana, 


07x 
08 a, 
09 k 


11 i 


Tabelle 2.11 Programm Seilkurve als Variationsproblem 


Eingabe 

000 r& LBL 
001 11 A 
002 47 CMS 
003 09% 9 
004 42 STO 
005 00 00 
006 41 RS 
00799 PRT 
008 72 ST* 
009% 00 00 
010 97 DSE 
011 00 0 
012 00 00 
013 06 06 
014 96 ADY 
015 43 RÜL 
0l& 08 08 
017 30 TAN 
018 42 STO 
019 10 WU 
020 01 1 
021 42 STO 
022 11 1i 
023 32 KIT 
024 43 RÜL 
025 09 09 
026 67 EU 
027 16 A" 
028 86 S5TF 
0293 01 01 
Start/Berechnung 
030 7& LBL 
031 1& RA° 
Neigungsänderung 
032 01 ı 
033 85 + 
034 43 RÜL 
035 10 10 
036 33 8 
03795 = 


038 34 TK 


039 65 & 
c4Uü 42 RCL 
ai 02 02 
042 693 8 
üÜdE 45 FÜL 
a4d4 Cl 
045 55 = 
de 42 RuL 


047 02 08 
04E 85 = 
049 44 sum 


sc 10 10 
Höhenänderung 

051 42 RÜL 
052 10 10 


067 34 
EEE 44 SUM 


069 040g 
Aufruf/Ausdruck 
corü 43 REL 
Ür 11 11 
ürs 17 BB! 


Zählererhöhung 


ar 
(rd 


ir 
IS 


Abfrage + Rücksprung 


pi 
44 
11 


1 


SuM 


11 


SBR/Ausdruck 
089 Te LEL 
ÜSE Tr E' 
er 8 IFF 
se 01 01 
05% £1 GTO 
cal 93 FRT 
031 42 RÜL 
az 0707 
a2 @@ PRT 
c34 42 FÜL 
[U U u 
096 92 FET 
(a7 42 REEL 
n9E 10 un] 
099 2EINMY 
100 ZC TAN 
101 ?9 FR 
102 43 RÜL 
102 04 04 
104 92 FET 
105 42 FÜL 
106 02 DE 
iur 820 + 
102 43 FÜL 


104 jo u 
110 ZEIMY 
111 ZU TAH 
112 39 C05 
1l3 95 = 
114 99 FRT 
115 92 AM 
1iE 9Z RTH 
x-Anzeige beik = 0 
il’ TE 

118 816 
119 42 RC 
120 n Ü 
121 32% 
122 43 RC 
122 0707 
124 FT GE 
129 99 PRT 
12& 43 FÜL 
127 07 v 
120 6 PAL 
129 92 RTH 


[PEN TENERIGEIGEIGFEFGHERN 
DR ı 5 SD re 


GL LLC 


m 
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Der Nachteil dieses Programms ist, daß es sich hier um ein Variationsproblem handelt, d.h. man 
muß mit mehreren Berechnungen, unter Variation von Anfangsneigungswinkel «| und Horizontal- 
zug H, sich auf den zweiten Aufhängungspunkt ‚einschießen’. Dies wird in 2.2.5—1- anschaulich 
vorgeführt. Zu diesem Zweck ist auch die Kennziffer k im Programm vorhanden. Für k = O wird 
nur der letzte Programmdurchlauf, also die Daten für den 2. Aufhängungspunkt, ausgedruckt. Hat 
man sich dann ‚eingeschossen’, wird mit k = 1 der gesamte Seilverlauf mit der Schrittweite Ax 


protokolliert. 
Es liegt nun nahe, die Lösung der Differentialgleichung oder eine Näherung dafür zu suchen, um die 
Startwerte des Iterationsvorganges besser bestimmen zu können. 


2.2.3 Die exakte Lösung 

Bringt man die Differentialgleichung (2.2.9) auf die Form 
EEE. 

v1l+ y2 H 


dann liefert eine erste Integration 


dx, (2.2.11) 


‚a 
arsinhy = yxte- (2.2.12) 


Sei der tiefste Punkt der Seilkurve die Stelle (xo, Yo), nach Bild 2.24, dann bestimmt sich die Inte- 
grationskonstante c; aus der Randbedingung y'ixo) =0 zu 


q 
cı = pe (2.2.13) y 


(x2. Ya) 


(x1,yı) 


Bild 2.24 
Seilkurve 


(xo, Vo) 


Wir erhalten somit 
v A q 
y- sinn (2x xol]. (2.2.14) 
Durch nochmalige Integration wird daraus 


H q 
y- qeosh (3 (x=x0)) +0. (2.2.15) 


Mit dem Koordinaten-Nullpunkt in (x, yı), Bild 2.24, bestimmt sich c, aus der Randbedingung 
y(x1}=0, zu 


cz = = cosh (3 xo). (2.2.16) 
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Die allgemeine Lösung lautet damit 


v=& (cosn (2-01) - cosh (2%)). (2.2.17) 


Die Unbekannten x, und yo ergeben sich mittels Randbedingungen aus den Gleichungen 


3 _Harcann (2) 
x garanhı (2.2.18) 
und 
H q 
va (oosh (H x) = 1). (2.2.19) 
LBLB /Start Tabelle 2.12 
Speicherplatzbelegung 
011 
02 q 
03 ya 
04 x, 
05 H 
06 ya/I 
07 gxo/H 
artanh (x) 
„Amftes 
2 1-x 
Ixl<1 
cosh (x) 
er +e’* 


Ausgabe * 
H 
xo 
Yo 


Bild 2.25 


Flußdiagramm zur exakten Lösung 


R/S / Programmende 
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Sind also die Aufhängungspunkte, die Seillänge und das spezifische Gewicht des Seiles bekannt, 
können alle übrigen Werte ermittelt werden. Der Horizontalzug H ergibt sich aus der Gleichung 


a2 /P-w im (22) (2.2.20) 

2H BE = sin >H). 2. 
Durch Kürzung auf die Form 

ucg = sinh (u) (2.2.21) 


2 2 
qx “-y ee E N R 
mit u= ze und co = V .—: erhalten wir die goniometrische Gleichung 
x 
2 


uco - sinh (u) = 0, (2.2.22) 


deren Lösung nur graphisch oder durch eines der üblichen Näherungsverfahren, wie Newton-Cotes 
oder Regula falsi, bestimmt werden. Mit der Größe u ist dann auch der Horizontalzug H gegeben. 
Ein ausführliches Beispiel dazu ist 2.2.5-1-—. 


Tabelle 2.13 Exakte Lösung 


u: 

ve LBL Ta Tu 473 085 55 + 
7 - gs + 08E 43 RCL 
47 CHE 1 Oi 1 ce? 05 05 
04 4 = “= 188 95 = 
42 5STD , B 25 + 089 42 STO 
oc 00 2i DE & DEN 090 07 07 
1 R=S Üge 9350 = oil 1 091 22 IN 
45 PRT 023 24 TH ran 092 23 LNX 
TE ST 42 REL 093 85 + 
cn cn ce 06 (94 43 RCL 
ar 1532 33 = 95 07 07 
TWIN] e3 LN# 096 94 +7- 
con 00 = £ E 097 22 INY 
ce 08 .L 43 RÜL 098 23 LNK& 
SE AI ge Je os 09 199 95 = 
5 ss = 100 55 + 

42 RCL 43 RÜL 101 02 2 

04 04 02 ma 102 95 = 

23 7 35 = 103 75 - 

02 2 a4 +8 - 104 01 A 

Er 85 + 105 95 = 

a 3 PRT 43 REL 106 65 x 
re 42 STO 4 04 107 43 RCL 
SEC EEE) us = 108 03 05 
oe ss + 109 55 + 
a Be [) u 3 110 43 RCL 
404 42 RCL [ 93 = 111 02 02 
Se m cd: 0% “3 FET lie 95 = 
ee Sg = 113 94 +7- 
42 FÜL Yo 114 993 PRT 

od 0 = Fo Te 115 98 ADV 

a = 08: 42 RCL 1l& 98 ADV 

4? STO od de 02 117 91 RS 
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Die wichtigsten Gleichungen wollen wir wieder in einem Programm zusammenfassen. Die Benutzung 
eines Näherungsverfahrens wollen wir dadurch umgehen, daß wir für sinh {u} die ersten beiden Ter- 


me der Reihenentwicklung setzen. Wir erhalten so 
3 
u 
- {u+-)=0Q 
UC9 (u 6 ) 0 


Durch Ausklammern von u folgt 


2 
1-5 =0. 


Diese quadratische Gleichung ist lösbar nach der Formel 


u=-y6lco-1), 


dabei ist nur die positive Lösung real. 


2.2.4 Straff gespanntes Seil unter vertikaler Einzellast 


Wenn ein Seil straff gespannt verläuft, läßt sich nach Bild 2.26 
annähernd 


dx 
U = Cosa (2.2.26) 
in die Gleichung (2.2.7) einsetzen 
y_a cd 
dx NH ecosa' (2.2.27) 


Eine zweifache Integration liefert 


a 2 + + 
= —— x? tc4X tn. 
YT IH cosa 1 2 


Aus der Randbedingung y (x = 0) = 0 folgt unmittelbar 
[7 Zi 0. 


Mit der Randbedingung y (xa) = y2 folgt 


2 
= 2 — 2 + 01Xa, 
y2° IH cosa 2 Cı%2 
also 
an Dez 
17% 2Hcosa 7 


Wir erhalten anschließend 


q 


a 2H cosa 


Y2 
KK? -xXaxX) + x. 
2 x 


(2.2.23) 


(2.2.24) 


(2.2.25) 


di 


dx 


Bild 2.26 Element eines 
straffen Seiles 


(2.2.28) 


(2.2.29) 


(2.2.30) 
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(x2, Ya) 


Xo. Yo) = (un, vo) 


Bild 2.27 Straff gespanntes Seil unter Einzellast 


Durch den Trick, in die beiden Aufhängungspunkte des Seiles, siehe Bild 2.27, jeweils einen Koordi- 
naten-Ursprung zu legen, erhalten wir folgende Randbedingungen: 

yv(x=0)=0, 

v(u=0)=0, 

yIx=xo)=h+tvlusuo). 


LBL C /Start Tabelle 2.14 

Speicherplatzbelegung 

Eingabe 01 rF 

X2, Y2 02 q 

qrF 03 ya 

04 %2 

05 yo 

. 06 x%o 

Eingabe 07. 

x. Vo 08 7} 


& (2 
a4 = arctan %o 


c0804 
vo M Yo-Y2 
X27%o 


Ausgabe 
H 


Bild 2.28 
Flußdiagramm zum gespannten Seil 
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Die Seilkurve ergibt sich für das jeweilige Koordinatensystem nach Gleichung (2.2.30) zu 


y 

VON = za Kor) BR (2.2.31) 
q v 

v(u) " Zsa, ul+ ou. (2.2.32) 


Aus einer Gleichgewichtsbedingung für annähernd ‚horizontalen’ Seilverlauf folgt 
F = Hy’(xo) + Hv’{uo). (2.2.33) 


Durch Einsetzen der Ableitung von (2.2.31) und (2.2.32) für xo = uo, erhalten wir daraus die Ge- 
brauchsformel 


2 cosaı cosa2 


HImS— sag nu (2.2.34) 


Mit Hilfe dieser Formel läßt sich der Horizontalzug näherungsweise vorbestimmen. Sie bildet den 
Abschluß der Berechnungsprogramme für das Themengebiet Seiltheorie. 


Tabelle 2.15 Programm gespanntes Seil 


Eingabe/Start a1 H: 
000 7r&E LBL DE 55 = 42 FEL 
001 13 C 027 43 RÜL DE DE 
002 47 cHS ST T5 sa - 
003 04 4 029 99 = 42 FRÜL 
004 42 5TO 030 22 IH or ur 
005 00 00 021 z0 TAH 32 205 
006 91 Rr5 Ü32 42 570 Ba + 
007 _ 99 PRT 033 F 07 Sa ı 
008 72 ST* 42 RL 
009 00 00 “2 ug 
010 97 DSZ 034 423 FÜL va 
611 00 DÜ gas 05 05 42 
012 00 00 DE 75 - DE 
013 0& 06 U37 "CL eg 54 
014 98 AD“ nbet) u: ii 55 

133 = a 5 
Eingabe/xo, vo ügc - Ei 
015 76 LEL c41 3 a 25 ; 
016 = RCL o4E I a 95 = jr u u 
01795 ADY IND 3C 1 oe. icon 75 - 
018 91 R#S c44 242 RÜL 101 43 REL 
019 99 PRT ü4S 2 02 DE 102 DE 0% 
020 42 STD CE 4 55 - 1032 54 5 
021 0& 08. gr u 1 104 54 
022 931 RS 45 nr BE Zu 105 95 = 
023 939 PRT 043 7 85 + 10& 393 FRT 
024 42 STD 0S0 E42 RÜL 107 £1 670 
025 05 05 üsl a3 0 0 108 43 RÜL 

Es — 
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2.2.5 Anwendungsbeispiele 
Io 


Zwischen zwei im Abstand von 50 m stehenden Masten der Höhen 30 m und 20 m, soll eine Leitung 
mit dem spezifischen Gewicht von 50 N/m eine lichte Höhe von 10 m nicht unterschreiten. Siehe 
Bild 2.29. 


Bild 2.29 
Hochleitung 


Die Wahl der richtigen Seillänge ergibt alle weiteren Startdaten. Bei einer Seillänge von |= 60m 
ergibt Programm B die gewünschten Werte: 


Eingabe: x2 


v2 


_A 
van 


Ausgabe: 


Mit diesen Startwerten berechnet Programm A iterativ die Seilkurve. Ein Problem bleibt die Be- 
stimmung des Startwinkels &,. Seine Variation liefert einen Winkel von a1 = - 57.8°. 
Nachfolgend lediglich der Ausdruck für die richtige Kurve. Sie ergab sich nach 9 Ansätzen. Dabei 
wurden die ersten Berechnungen mit einer größeren Schrittweite durchgeführt. 


Eingabe: 14 = 
-57,8 =0ı 
Ü, =x 
=, =yı 
Sl, =x2 
U, =h 
l. = x 
1200. = H 
Sl, =q 
Ausgabe: =si= 
=x= 
=yi® 
=aj= 
= |; = 
=5= 
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6. 

&. 

22. 02307994 
-49, 16327807 
9. 391006159 
1835. 31247 


? 
r 
20. 92959082 
-47. 55696358 
11. 47280184 
1778. 154822 


19. 89784319 
-45. 89521812 
12. 90963975 
1724. 205491 


Er 

ch 

18. 2259638 
-44, 18297067 
14. 30411081 
1673. 365273 


10. 
10. 
18. 01218738 
-42. 42033175 
15. 65872756 
1625. S40de1 


11. 
il. 
17. 15485333 
-40. 60760835 
16. 97593005 
1580. 642655 


12. 

12. 

16. 35240271 

-38. 74532201 
18. 25808724 

1538. 585601 


12. 

13. 

15. 60337523 
-36. 83421669 
19. 50750393 
1499. 300024 


14, 
14. 
14. 90640692 
-34. 8752775 
20. 7264235 
1462. 703485 


15. 

15. 
14, 26022686 
-32. 8697356 
21. 217032203 
1428. 730237 


126€ 16365548 
-30, 31309255 
23. 08146211 
1397, 316093 


i 
Tofu 
NE 5 DE 


11. GrOSa7rE 
-17. 


RER 
2 
10. PaSTr89r 
-19. 38257974 
30. 66904852 
1244. 533922 


6 8500024. 


“oO 


- Achr 
no 


BEI 


Für das ‚Einschießen‘ auf den zweiten Aufhängungspunkt wird die Kennziffer k = 0 eingesetzt. Es 
erfolgt nur der Ausdruck der letzten Berechnung. Bild 2.30 gibt die ersten 4 Versuche wieder. Die 
ersten 3 wurden mit der Schrittweite 5 und der letzte mit 1 durchgeführt. Die richtige Kurve mit 

ausreichender Genauigkeit zeigt Bild 2.31. 


Bild 2.30 ‚Einschießen’ auf den zweiten Aufhängungspunkt 
u} Ü, 
4 EU 
1. Il 
Kirn Kun 
Sc, Sl, 
LE, Ü, 
c 1200 
So, 


a, =-57,8°/Ax=1 


50. 
SD, 
15. 293379478 
40. 54601101 
51. 41850058 
1579. 188017 


Bild 2.31 
Seilkurve der Hochleitung 


-2_ 


Auf einer Länge von 1000 m soll eine Seilbahn einen Höhenunterschied von 300 m überwinden. Die 
Last beträgt 20kKN und das spezifische Tragseilgewicht ist 100 N/m. Wir wollen für diese Seilbahn 
einen maximalen Durchhang von 50 m zulassen. Damit können wir iterativ für jeden Punkt einer 
Parallelen im Abstand 50 m zur Steigung 


300 


x - 50 = 0.3x - 50 


y” 1000 


1000 m 


Bild 2.32 


Seilbahn 
(x4,Yı) 
20 kN 


Seilkräfte und Seillänge bestimmen. Doch zunächst errechnen wir den notwendigen Horizontalzug 
für äquivalente Stellen auf der Parallelen mittels Programm C. 


Eingabe: 1000. =x% 
200, =ya 
100. =q 
eocöoo, =) 
ES, 
Ausgebe: 350. = xg 145, 
55. = yo lea 147552, 44 
14rE32. 7732 =H 


Wie zu erwarten war, erhalten wir für die Koordinaten (500, 100) den größten Horizontalzug. Da 
hier für das Seil auch der größte Durchhang zu erwarten ist, wollen wir für mittigen Lastangriff die 
Seilkurve ermitteln. Die direkten Winkel! betragen &; = 11.3° und a = - 21.8°. Wenn wir diese 
etwas neigen, bekommen wir gute Startwerte. 10 Versuche ergaben dann folgende Daten: 


Linke Kurvenhälfte: 


Eingabe: 1, =k 
2. i =-04 

U, =x 

Mo =yı 
Ol, = %2 

I 
zu. = Ax 

je2llod,. =H 

100, =a 

Ausgabe: 


20. 0E40096% 
162295. 60 


43 


5 

160. 

14, Pri1rr96l 
7, T29215436 
160. 7444 


163520. 2983 


21. 
420. 
72, 31598446 
16. 68947164 
427. 4359057 
1659232. 8876 


22. 

440. 

79. 16987314 
17. 3641732 
445. 3909842 
169544. 2917 


23, 
460. 
35. 68229929 
18. 03635315 
469, 4244691 
170451. 4124 


24. 
480. 
92, 45423272 
18. 70592929 
490, 5398446 
1rildd. 3415 


25. 
Ss00. 
22, 48668256 
19. 37282115 
Ss11. 7402066 
171933. 1744 


Rechte Kurvenhälfte: 


Eingabe: 


=k 
= 
x 
=yı 4. & 
=x2 320. E80, 
=| 2596. 6440656 236. 6420322 
= Ax 27,.51403013 26. 24949603 
=H 1?0. 4264744 235. 2094279 
=q 1S277E, 4275 180742, 807 
Ausgabe: en = ” 
=“ 860, 
=y ; 227. 0544696 
=a; 2.8 23. 61204032 
sh 212. 2095169 257. 3887281 
=$; 181745, 6756 179767, 6639 
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2. 

340, 

217. 7405498 
24, 97120658 
279. 4511192 
178820. 0924 


3 

820. 

208. 6338305 
24. 32705813 
301. 3999845 
177899. 943 


10. 

800. 

199, 927911 
23. 67966028 
323. 2396893 
177007. 0704 


11. 

r80, 

191. 4264322 
23. 02908032 
344. 9TUIB1E 
176141. 3334 


wıpın 


in Tı 


14387589 


172215. 0245 


In Ten 


iz 

1 4 B 
Se. 2II4TTE 
169549, 9158 


2l. 
>50. 
120, Jei11d6r7 
16. 26444255 
537, 1299097 
1883942, 2617 


[un Ar ur ur RER 


Ina ® 


efgermın 
[um = DEE u 3 zer SS 2 


nie nıns 
STERNE 


- Mminm® 


‚2902854 
239. 2294221 
18£300. 8275 


Bild 2.33 
Seilkurve der Seilbahn 
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Die ermittelte Seilkurve zeigt Bild 2.33. Eine Betrachtung der Stelle, an der die Kraft angreift zeigt 
uns jedoch, daß diese Berechnung nur eine grobe Näherung ist. An dieser Stelle ergibt sich das in 
Bild 2.34 dargestellte Kräfteverhältnis. Nach der Gleichgewichtsbedingung können die an der Last- 
angriffsstelle vorhandenen Seilkräfte S und S’, mit ihren vertikalen Komponenten V und V', näm- 
lich nur eine etwas geringere Vertikallast V’ tragen. Es ist 


V=S sin(19.4) = 56999.4 N 
V’=$’ sin (13.6) = 39304.0 N 
und damit 
BO. 
L’ = 17695.4 N <L = 20000 N vi! 
U 
Bild 2.34 Kräfte am Lastangriffspunkt 
Beer. 


Dieses Beispiel soll Ihnen einen Einblick geben, welche komplexen Probleme sich mit dem aufge- 
stellten Programm bewältigen lassen. 

Dazu betrachten wir eine Seilanordnung mit schief angreifender Einzellast und unterschiedlichen 
Streckenlasten. Ein rechnerischer Ansatz ist nicht möglich. Unter Annahme des Horizontalzuges 
und des Anfangswinkels ergeben sich die ersten 10 m linearen Abstand, wie in Bild 2.34 dargestellt. 
Über diese Entfernung soll lediglich das spezifische Seilgewicht wirken. 


(x1,yı) 


Bild 2.34 1. Seilstück 


Eingabe: 1 ö & 


Ausgabe: 1. 
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RET 


2 


8. 

39. 41299024 
46, (9616287 
13. 20279154 
1442. 064866 


9 
9 
38. 84607901 
-44,. 03622086 
14. 6938046 
1391. 013058 


10. 

10. 

37. 54871843 
-41. 30355151 
16. 03740275 
1343. 598157 


Soll dort nun eine Kraft mit einem horizontalen Anteil von F},; = 300 N und einem vertikalen 
Anteil F,=400N angreifen, so ergibt sich aus der Gleichgewichtsbedingung der Kräfte 


V = Ssina, = 897.4 N 
H’=H+Fy = 1300 N 
V’= V-Fy = 197.4N 


’ 


a, = arctan — = 8.63° 
H 


Damit haben wir den Ausgangswinkel für den weiteren 
Seilverlauf. Nach weiteren 10 m linearen Abstand 
wollen wir die Streckenlast ändern. 


Bild 2.35 Kräfteverhältnisse 


(,yı) 
(x2, v2) 


Bild 2.36 2. Seilstück 


Ausgabe: 


Betrachten wir nun einen Ansatz für eine Streckenlaständerung. Für ein Seilelement mit veränderter 
Streckenlast ist nach Bild 2.37 


H=H' 
V+qAl=V, 


Das heißt, wir können mit den zuletzt berechneten 
Daten weiterrechnen und müssen nur die Streckenlast 
qin q ändern. Dies soll für weitere 10 m linearen Ab- Bild 2.37 Streckenlaständerung am 
stand gelten. Seilelement 
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gl (x2. v2) 
u) Bild 2.38 3. Seilstück 
a’ 


x,Yı) ee 
192 10m 


Eingabe: Ne =kK ae 
u 46. 52452124 
38 . 62. 92374956 
n. 38. 22682475 
26, u 2817. 628555 
= Ax 
= 9 
=q 29, 
49. 1808962 
#i 56. 97617 
IE N 40. 69486607 
=y; 3208, 453713 
=, 
a 10. 
"5 30. 


51. 71696982 
53, 22565508 
43. 51424262 
3665. 189516 


Abschließend sollen noch weitere 10 m linearer Abstand mit dem spezifischen Seilgewicht 
q = 50 N/m folgen. 
Den so ermittelten Seilverlauf zeigt Bild 2.39. 


8 

38, 

77, 09860749 
74, 56736274 
70, 09768586 
4535. 285055 


Eingabe: 


1200, 
So. 


Er 

39. 

0. 825335879 
75. 12279236 
73. 39511093 
S06ß, 652584 


Ausgabe: 


10. 

40. 

54. 74241106 
75. 67793079 
78. 0376017 
3255. 238013 
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Bild 2.39 
Seilkurve eines kombiniert belasteten Seiles 


2.3 Reibung 


Eine bei jedem technischen Prozeß beteiligte und leider auch ebenso notwendige Kraft, ist die Rei- 
bungskraft. Sie entsteht an der gemeinsamen Berührfläche zweier Körper, wenn diese sich gegenein- 
ander bewegen oder bewegt werden sollen. Ihre Wirkrichtung ist immer der Bewegung entgegenge- 
setzt. Werden die Körper aus der Ruhe heraus bewegt, so muß die Haftreibung, ansonsten die Gleit- 
reibung überwunden werden. 

Eine allgemeingültige Gleichung für alle Reibungsfälle ist leider nicht gegeben, so daß ich mich nach- 
folgend auf spezielle Anwendungsbeispiele beschränke. 


2.3.1 Keil und schiefe Ebene 


Die in der Natur aufs vielfältigste verwendete schiefe Ebene, läßt sich auf folgendes Grundprinzip 
zurückführen. Nach Bild 2.40 wird ein Körper | auf einem Körper Il aus der Ruhe oder gleichförmig 
mittels der Kraft F, bewegt. Eine Gewichtskomponente erzeugt die für die Reibung grundlegende 
Normalkraft. Aus der allgemeinen Gleichgewichtsbedingung leitet sich die analytische Lösung 
sin{a@tp) 


1 cos(ßFp) (2.3.1) 


Fa = 


Bild 2.40 
Grundprinzip schiefe Ebene 
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ab. Darin gilt das obere Rechenzeichen für eine Bewegung nach oben, daß untere für die Bewegung 


nach unten, für Körper I. Ist es eine Bewegung aus der Ruhe heraus, so ist 


p = arctan ko 
und anderenfalls 


p = arctanı. 


(2.3.2) 


(2.3.3) 


Für dieses einfache Programm wollen wir uns ein Flußdiagramm sparen und nur in Tabelle 2.16 und 


2.17 die Speicherplatzbelegung und das Rechnerprogramm wiedergeben. 


Tabelle 2.16 
Speicherplatzbelegung 


01 F, 05 ı + 1 für eine Bewegung 
02 ß nach oben 
03 « - 1 für eine Bewegung 
04 u, Ho/p nach unten 


Tabelle 2.17 
Programm schiefe Ebene für den T158/59 


Start 
cooc FE LEL oil am oc os 0 
oo 11 RA 012 97 IS2 us = 
002 47 CHE 2 00 00 42 57 
a2 05 5 ol4 42 REL a4 0 
co4 42 ST (iS 25 AM +42 RC 
Ds Cu 0 u 0 
Berechnung gu + 
Eingabe 016 : RL 43 RC 
Oce Te LEL ur cs 4 0 
our 42 FÜL Di IN" 2a = 
u [17 THAH za SI 
009 LER : = = 
OLc ci 42 RÜEL BR L 


2.3.2 Gewindereibung 


CTte$O in 


+nr 


035 
36 
(37 
1138 
u39 
40 
c41 
ü42 
043 
c44 
c45 
üde 
047 
ü48 


Eine Gewindefläche entsteht durch gleichzeitige Rotation und Translation einer ebenen Kurve, nach 


Bild 2.41 
| ) 
gr Bahn der Schraubenkurve 
377 Bild 2.41 
Gewinde 


Schraubenkurve | 
(Gewindeform) | 

) 
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Nach der Form der Kurve unterscheidet man Spitz-, Flach-, Trapez-, usw. Gewinde, also die Gewin- 
deform. Bei der Gewindebewegung tritt die Normalkraft senkrecht zur Schraubenfläche nach Bild 

2.42 auf. Sie erzeugt eine Reibkraft und ein daraus resultierendes Reibmoment. Die Reibkraft wird 
dabei mittig zur Schraubenfläche abgenommen. Auch hier ergeben sich aus der allgemeinen Gleich- 


gewichtsbedingung die Endgleichungen 
F = (cosa cosß tsina@tanp) N (2.3.4) 
und 


M = (sinacosßFcosatanp) N. (2.3.5) 


Beim Flachgewinde ist ß= O0 und damit 
cosß = 1, womit als Spezialfall folgt 

F = (oosa +sinatanp) N (2.3.6) 
und 


M = (sina # cosatanp) N. (2.3.7) 


Gewinde- 
steigung 


( 


Flankenneigung 


Bild 2.42 Kräfte am Gewinde 


Da die Last allgemein bekannt ist, soll unser Programm die Normalkraft bestimmen und das daraus 
resultierende Reibmoment. Das Vorzeichen richtet sich nach einer Bewegung gegen die Last (oberes 


Zeichen) oder mit dieser (unteres Zeichen). 


Tabelle 2.18 

Speicherplatzbelegung 

00 Zähler 06 | + 1 bei Bewegung gegen die Last 

Ola - 1 bei Bewegung mit der Last 

02 F 

03 u 

04 ß 

05 a 

Tabelle 2.19 

Start Berechnung 

000 76 LBL le 43 RLL 

001 II A Hr 05 08 

002 7? CMS 18 

003 06 6 117 

004 42 STO Li 

005 00 ou 
Eingabe x 
006 7E& LBL ne 
007 43 RCL Br Des 
008 91 Res h En 
00% 93 PRT ir ee 
010 72 ST# ; EBEN 
ii 00 u u . 5 
012 97 152 ee en En 
013 00 00 ee 3a LOS 
014 43 RCL et ne 
015 98 ADY = 7 Aa 
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2.3.3 Seilreibung 


Durch die Bewegung eines Seiles über eine zylindrische Rolle, wird durch den Andruck des Seiles 
eine Reibkraft erzeugt. Die Betrachtung eines infinitesimalen Seilelements nach Bild 2.43 liefert 
durch die Gleichgewichtsbedingung 


dFg = dF, = udFn- (2.3.8) 


F 
e = Bild 2.43 
Seilreibung 
I 
Fu 
Die Normalkraftkomponente ergibt sich angenähert durch 
dfn = Fs da. (2.3.9) 
Die Betrachtung über den ganzen Umschlingungswinkel liefert das Integral 
Fs2 & 
ds [  _ 
[ — = |naa, (2.3.10) 
: Fs 
Fsı o 
womit folgt (die Integrationskonstante entfällt durch Randbedingungen): 
F 
in(2*)- us (2.3.11) 
Fs2 
Aufgelöst nach den Kräften erhalten wir die Eytelweinsche Gleichung 
Fo = Fsa .grä (2.3.12) 


Da ua > 0, ist die Seilkraft Fg] > Fga; die Kraft Fs} wirkt also immer der Bewegungsrichtung des 
Seiles nach. Der Wert von « ergibt sich aus der einfachen Umrechnung 


Rn o 


za (2.3.13) 
180 


ü= 


Tabelle 2.20 
Speicherplatzbelegung 


01 ua 
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Tabelle 2.21 
Programm Seilreibung zum T158/59 


Start 

000 TE LBL COS 42 (12 E85 dei 

001 ii A ccge 0 Cd 59 [M) 

002 47 CMS car 31 JS 95 = ) 
cüS 97 Ole 43 FRO ie) 

N co 8% ol Gl 0 u2! 
EINE TERSERRUNg io 01 DIE a1 R.S 7 
003 91 RS oil 08 01a 23 FR? ] 
004 9% PRT olz oc u zo 32 AD% in] 


2.3.4 Anwendungsbeispiele 
-1- 


Auf einer schiefen Ebene von 20° liegt eine Last von 20 kN. Die Haftreibung beträgt 0.3. Die Last 
wird durch eine Kraft um 5° zur Ebene gehalten. Gesucht ist die Haltekraft bzw. die Kraft die be- 
nötigt wird, um die Last nach oben zu ziehen. 


l. nach oben nach unten 
0. = ko Ko 
20. a a 
Er ß 
20, Fj Fı 
12. 20586228 Fa 1.2392 Fa 


Um die Last zu halten wird eine Kraft von 1.24 kN benötigt. Um sie nach oben zu ziehen müssen 
12.21 kN aufgebracht werden. 


-2- 
Auf einen Keil nach Bild 2.44 wirken Normalkräfte von I kN. Bei einer Haftreibung von 0.3 ergibt 
sich die Frage, welche Kraft zum Eintreiben des Keiles und welche zum Herausziehen benötigt wird. 


Eintreiben: Herausziehen: 
-1. 
ko NE yo 
a 12.0 
ß MB 
2-N=1KN=F, Fı t. FR 
1, S4EE3EFa5 Fz -. SSL 3SSTÜES Fa 


Bild 2.44 Keil 


Zum Eintreiben werden 0.55 kN und zum Herausziehen 0.031 kN benötigt. 
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eg 


Auf einem spitzen Bewegungsgewinde von 50° Flankenneigung und einer Steigung von 5°, sowie 
einem Flankenradius von 3cm bewegt sich eine Last von 10KN. Bei einem Reibungskoeffizienten 
von 0.2 ist nach der Normalkraft und dem Reibmoment gesucht. 


Last nach oben Last nach unten 
l. 
0 a 
22.08 ß 
[u Far u 
10. FikN F (kN) 
=. acm acm 
N (kN) N (kN) 
M (kNcm) M (kNcm) 


Die Bewegung nach unten bewirkt also eine größere Belastung des Gewindes. 


ide 
Eine Last wird über eine feststehende Umlenkrolle nach oben gezogen. Das Seil, an dem die Last 


hängt, umschlingt die Rolle mit 120°. Der Reibungskoeffizient beträgt 0.2. Die Last wiegt 500KN. 
Es ist nach der Zugkraft gefragt. 


una 
ner 


FEL, LESZ1LE 


u 
n 
o 
- 
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Während wir zuvor in der Statik den Spezialfall behandelt haben, daß äußere Kräfte am starren 
Körper keine Bewegung hervorrufen, wird in den nachfolgenden Kapiteln die durch Kräfte hervor- 
gerufene Bewegungsänderung untersucht. Dabei wurde keine streng thematische, sondern vielmehr 
eine für den Rechner interessante Einteilung vorgenommen. 

Bei der Kinematik sieht man von der Masse und den am Körper angreifenden Kräften ab und unter- 
sucht nur deren Geometrie der Bewegung. Wir betrachten damit die Lageänderung eines Punktes 
über der Zeit. Dies geschied bezüglich eines passend gewählten Koordinatensystem. 


3.1 Numerische Behandlung von Differentialgleichungen der Bewegung 


Soweit in diesem Buch die Problemanalyse auf eine Differentialgleichung führt, siehe auch Seil- 
theorie, habe ich deren numerische Integration, um Sie nur mit einer Methode zu konfrontieren, 
nach dem Euler-Cauchy-Verfahren durchgeführt. Daher will ich diese Methode kurz erläutern. 

Sei nach Bild 3.1 y= f(x} die analytische Lösung der allgemeinen Differentialgleichung y’ = f (x, y). 
Aus der Differentialgleichung folgt die Anfangsbedingung yo = f (xo, Yo) als bekannter Wert. Geht 
man nun auf der Abszisse um die Schrittweite Ax weiter, so läßt sich die auftretende Ordinatenzu- 
nahme Ay annähernd durch den Wert Ay, beschreiben. Es gilt 


Ayn 
"= 
vn (3.1.1) 


Aus der Grenzwertbetrachtung folgt 


= Jim Ein, lim -— = — (3.1.2) 
Y Ax>0 AxX Ax>0 Ax dx 2 


und damit die exakte Lösung. Bei diesem Verfahren wird also der Differentialquotient 
dy 
dx 


exakte Lösung 
v=fo) 


durch den Differenzenquotient 


Ay numerische Lösung 
Ax 
Ay 
| Ayn 
Y Ax 
x 
Bild 3.1 


Das Euler-Cauchy- Verfahren 
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angenähert. Wählt man Ax also genügend klein, kommt man der analytischen Lösung beliebig nahe. 
Den Nachteil dieses Verfahrens zeigt Bild 3.1 ebenfalls deutlich. Mit zunehmenden Schritten ent- 
fernt sich die numerische immer mehr von der analytischen Lösung. Wer genauere Verfahren be- 
nutzen möchte, kann diese unter [2] und [12] nachlesen. Als Stichwort möchte ich das oft benutzte 
Runge-Kutta-Verfahren nennen. Betrachtungen zur Fehlergröße der einzelnen Verfahren entneh- 
men Sie bitte ebenfalls der Literatur. 

Speziell in den nachfolgenden Fällen wird das Euler-Cauchy-Verfahren zweimal Anwendung finden, 
da Massenpunkt oder starrer Körper zur Änderung ihres Bewegungszustandes durch Kräfte F;, 
j=1,..., n gezwungen, diesen ihre träge Masse entgegensetzen. Als Gleichung 


n 
DFi=F,=-ma, (3.1.3) 
jal 
Die Momentanbeschleunigung ergibt sich als Differentialquotient 
dv_. 
Be (3.1.4) 


so daß die 1. Anwendung des Euler-Cauchy-Verfahrens, 
Av=- At, (3.1.5) 


die während der Zeiteinheit At auftretende Geschwindigkeitsänderung liefert. Für die Momentan- 
geschwindigkeit gilt weiterhin der Differentialquotient 


ds _. 
Far (3.1.6) 
so daß eine 2. Anwendung des Euler-Cauchy-Verfahrens, 
As=vät, (3.1.7) 


den während der Zeiteinheit At zurückgelegten Weg annähernd beschreibt. Die Momentangeschwin- 
digkeit wird dabei aus Einfachheitsgründen für die Berechnung durch die Ausgangsgeschwindigkeit 
für das Zeitintervall ersetzt. 


3.2 Bewegung des materiellen Punktes 


Die Lage eines Massenpunktes zum Zeitpunkt t wird durch den Ortsvektor r (t) bezüglich eines 
raumfesten Koordinatensystems bestimmt. Siehe Bild 3.2. In einem Zeitintervall At verändert er 
seine Lage um Ar (t, At) nach r (t + At). Der Übergang auf ein infinitesimales Zeitintervall dt liefert 
zum Zeitpunkt t die Momentangeschwindigkeit 


3 
EM rtt+An-rtt dritt 
nz PN At dt 
Ar (t, At) 
(3.1.8) ng 
r(t + At) 
Bild 3.2 
Raumbewegung 1 
2 
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Für das Bogendifferential wird mithin 
dr = v(t)dt. (3.1.9) 
Analog definiert man die zum Zeitpunkt t vorhandene Momentanbeschleunigung 


vir+Ar)- vi) dvd 


Ki er (3.1.10) 


alt) = lim 
At>0 
Die Momentangeschwindigkeit und -beschleunigung ergeben sich also als erste bzw. zweite Ablei- 
tung der Ortskoordinaten nach der Zeit. 
Die Bewegung eines Massenpunktes auf gekrümmter Bahn, läßt sich für ein infinitesimales Bogen- 
element, das eine Drehung um den imaginären Bahnmittelpunkt M mit dem Krümmungsradius p 
auffassen. Dieser ergibt sich, unter Betrachtung von Bild 3.3, aus der Gleichung 


Ivo rl. (3.1.11) 


Bild 3.3 
Bahnkrümmung 


Die bei der Bewegung aufgespannte Kreissektorfläche liegt in der sogenannten Schmiegungsebene. 
Während der Geschwindigkeitsvektor tangential zur Bahnkurve verläuft, besitzt der Beschleunigungs- 
vektor eine Tangential- (3) und Normalkomponente (a„). Das Anträgen der zu den Ortsvektoren 
r(t) einer Bahnkurve gehörenden Geschwindigkeitsvektoren v(t) bezüglich eines beliebigen Ge- 
schwindigkeitspoles (0), siehe Bild 3.4, liefert die in Bild 3.5 wiedergegebene Konstruktion. Die so 
entstehende Raumkurve wird als polarer Hodograph bezeichnet. Danach zerfällt die Geschwindig- 
keit ebenfalls in einen Normal- und Tangentialanteil 


Av = Av, + Av.. (3.1.12) 


Hodograph 


Bild 3.4 Raumgeschwindigkeiten Bild 3.5 Polarer Hodograph 
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Ihre infinitesimale Betrachtung liefert dann durch Differentiation nach der Zeit die vorangegangene 
Behauptung 


a=-anta. (3.1.13) 
Ihre Größen ergeben sich aus 
22 
v S 
nn" (3.1.14) 
und 
a =v4=8. (3.1.15) 


3.2.1 Bewegungsdiagramme 


Die Bewegungen eines Massenpunktes werden in der Regel durch Zeit-Weg-, bzw. Zeit-Geschwindig- 
keits- oder Zeit-Beschleunigungsdiagramme dargestellt. Auch eine Kombination der Bewegungs- 
größen ist üblich. 
Da in den weitaus meisten Fällen die Raumkurve durch einparametrige Gleichungen der Form 

ft; = fi(u), i=1,2,3 (3.1.16) 


oder durch Meßwerte gegeben ist, womit sich auch eine Interpolations- oder Approximationsglei- 
chung aufstellen läßt, soll dies die Ausgangsbasis für ein Programm sein. 
Der Ortsvektor eines Bahnpunktes ergibt sich in seiner Komponentendarstellung aus 


r= Verf; (3.1.17) 
i=1 
wie Bild 3.6 es wiedergibt. Die Momentangeschwindigkeit in diesem Punkt bestimmt sich angenähert 
durch die Ortsveränderung des Massenpunktes in der Zeitdifferenz At. 


f: 
& 7, (3.1.18) ea 
u 


Bild 3.6 
Räumliche Bahn 


Mit der Größe von At läßt sich die exakte Lösung beliebig genau approximieren. Eine zweite An- 
näherung des Differentialquotienten durch den Differenzenquotienten liefert die Momentanbe- 


schleunigung 


dvtu)_ Avtu) LS, Av; 
m E wese za a (3.1.19) 
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Die Beträge von Wegzunahme, Geschwindigkeitsänderung und Beschleunigung ergeben sich nach 
dem pythagoräischen Ansatz über die Komponenten k; 


1...|= DK. (3.1.20) 


Außer der Bestimmung der Raumkurve mit der Zeit als Parameter durch die Funktionen 
= ft), i=1,2,3, (3.1.21) 


gibt es noch die Bestimmung durch einen geometrischen Parameter, z.B. des Winkels p. Dessen zeit- 
liche Veränderung wird dann durch die Gleichung 
y=yplt (3.1.22) 


bestimmt. Die Ableitungen eines Ortsvektors nach der Zeit ergeben sich in diesem Fall nach den 
Regeln der Differentiation angenähert aus 


v(t) = Ay Ar (3.1.23) 
und 
_ Au &p 
alt) = Ay At‘ (3.1.24) 


Wir erhalten damit einen in Bild 3.7 dargestellten Berechnungsrumpf, in dem die Unterprogramme 
A’, ..., D’ die funktionalen Verhältnisse des jeweiligen Problems wiederspiegeln. Auf diese Weise 
lassen sich z.B. Schubkurbel-, Kolben- oder Nockenbewegungen analysieren. 


Tabelle 3.1 

Speicherplatzbelegung 

07 2 08 2 09 u, 
10 fa. 11 ff. 12 fı,.ı 
13 fa} 14 fa 15 fı; 
16 v3 -1 17 v2 -1 18 V1j-1 
19 v3 20 V3j 21 Vij 
22 a3, 23 a2j 24 47j 
25 At 26 t; 27 5; 

28 v 29 4 


In diesem Programm wurde bewußt aus Platzgründen auf eine Berechnung des Krümmungsradius 
und auf den Normal- und Tangentialanteil der Beschleunigung verzichtet. Wird die Berechnung ge- 
wünscht, so läßt sich dies mit den zuvor erläuterten Formeln leicht programmieren. Auch läßt sich 
das Programm durch indirekte Programmierung auf weit weniger Programmschritte (für TI58) 
bringen. 

Dieses Programm soll das einzige in diesem Kapitel bleiben, da sich die praktische Nutzanwendung 
erst im nächsten Kapitel ergibt. 
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LBL A /Start 


Eingabe 


to 


At Hier gilt folgende Zuordnung 


A: fılu,) 
en ur B': f2 (u) 
c: f3 (u) 


a EHRE 1 Aufruf von SBR D’ 


A;ehich,nn 


S=S.-1 +. AS, 


Als SBR E kann Umfang der Ausgabe 
beliebig festgelegt werden. Auch 
Umrechnungen sind möglich. 


Bild 3.7 Fiußdiagramm zur allgemeinen Bewegung 
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Tabelle 3.2 
Programm allgemeine Bewegung 


Start + Eingabe 


000 7& LBL 
001 11 A 
002 47 CMS 
003 91 RS 
004 2% PRT 
005 42 5TD 
006 26 26 
007° 931 RS 
008 99 PRT 
009 42 STO 
010 25 25 
o0ll 98 ADI®Y 
012 61 GTO 
013 12 & 
Ausgabe + Zeitzähler 
014 7& LEL 
015 43 RCL 
0iß& 15 E 
01? 00 0 
018 42 STD 
019 08 03 
020 42 STD 
021 07 ? 
022 43 RÜL 
023 25 25 
024 44 SUM 
025 26 26 
Start-Berechnung 
026 7& LEL 
027 12 68 
u: 

028 43 REL 
029 28 25 
030 19 D' 
031 42 STD 
032 09 0% 
fy/äfi, H 

033 16 A" 
034 42 STO 
035 15 15 
036 75 - 
037 43 RCL 
038 12 12 
039 95 = 
040 42 STD 
041 21 21 
042 33 %2 
043 44 SUM 
044 08 03 
045 43 RÜL 
046 09 09 
047 17 B" 


ü4s 
49 
050 
051 
052 
033 
054 
055 
Use 
057 
US$ 
059 
060 
061 
ja pe 
(63 
064 
065 
Ü6E 
067 
065 
069 
url 
rl 
OrE 


Asılsj: 


073 
ur74 
075 
üre 
Or? 


Yjj: 

u78 
079 
0S0 
Ds1 


Ar} 
al 
'D 


2 


= 


| 
12 
L 
E 


u 
DeKee' 


N Rai + 
on 


Alt/neu Verschiebung 


1393 24 [5 SO 43 RÜL 
100 42 STD si 15 15 
101 2825 S2 42 ST0 

s3 12 12 
3jj 54 43 RÜL 
102 42 RÜL 39 14 
aa 21 2 Ei To 
104 75 Sr 11 
105 42 RÜL ss REL 
106 E Fe} 3 13 
cr 4] zn 
ı0E 5 el 


Pfipe 

Gaga 

a 
'T£eO 


DE re EEE EEG EEE GR EIGEPETEF GEF GEF GEF REFGREFG 


& 
An} 
7 a m \ 7 7 En 7 
> hGkH4HtieD ie m m4daNlen ı Bin Haller ı Baımı 


110 

111 

112 t: 

113 43 RCL 4: 

114 AT 17 el 

11s ss 42 

1le& 42 STD Ir 

ir 8 3 42 

118 43 PEL 13 

119 19 19 42 

120 78 t6 

' ai daR L Durchlaufabfrage für 
122 u Korrektur 

123 140 
be} [1 

123 2 175 44 sum 
E> - I76 05 05 
155 s 177 43 RÜL 
125 a Ir5 05 05 
139 Ira BE MIT 
30 120 01 1 
131 181 87 EU 
132 IB, 33 
133 123 02 2 
33 184 Er Eu 
= 185 75 - 
ri 126 #1 GT0 
J 157 43 RÜL 


Korrektur zu den 
beiden ersten 
Durchläufen 


[ONE ACH ARE A Eee 
oben 


ee 


[2 u] 

42 5T0 
32 sTo 
78 LEL 
75 - 

oo Ü 

42 STO 
23 23 
&1 GT0 
43 RÜL 


3.2.2 Anwendungsbeispiele 
BE 


Ein Massenpunkt bewegt sich auf einer Schraubenlinie mit 


x = 2cosp 
y = 2sinp 
z=2% 


Und zwar mit der Winkelgeschwindigkeit 
yit) = 2t+yo 


mit p ingrd und t in Sekunden. Zur Anfangszeit to = O befindet sich der Massenpunkt an der 
Stelle yo = 45°. 


Bild 3.8 
Schraubenlinie 


x“ 


Bild 3.8 zeigt den Verlauf der Bewegung. Die Unterprogramme für dieses Problem zeigt Tabelle 3.3. 


Tabelle 3.3 
Unterprogramme zu — 1 — 


[PEIME 

wi wu 248 43 RCL 

32 RTH 249 14 14 

250 99 PRT 

251 43 RCL 

”E LEL 252 13 13 

180: 253 99 PRT 

be] 0 254 43 RCEL 

24 CE 255 09 09 

EIN 256 9% PRT 

I 857 43 RCL 

= 258 27 27 

u va a 2 259 99 PRT 
et 000 2 260 43 RCL 
En Fr 4 24: 261 28 28 
34 cE = ETH 243 28 262 93 PRT 
38 31 244 935 PRT nn BE as 
65 x 245 42 RCL Ses 39 PRT 
02 2 TE LEL 246 15 15 SEE 98 ADV 
Im D ea 29 PRT 26° 92 RIN 
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Das Programm lieferte damit folgende Werte: 


Eingabe: 0. _ to 
l. -a 
Ausgabe: o. -1 


1.414213562 eh 
1.41g213562 Via 
1. Sr0r96327 "23 "f3j 


45. 9 
0. 5 
0. 4% 
0, 3 
1 


1. 36399672 
1. 462707403 
1. 640609497 

> 
. 0987282258 
0987282258 

0. 


Da die Bewegungsverhältnisse leicht überschaubar sind, soll uns deren graphische Darstellung auch 
nicht weiter interessieren. 


-2- 

Über die Bewegung eines Massenpunktes sind folgende Werte gemessen worden: (in m) 
t x y z 

0.0 0.0 0.0 0.0 

0.1 0.125 0.062 0.020 

0.2 0.215 0.180 0.111 

0.3 0.430 0.192 0.204 

0.4 0.522 0.201 0.312 

0.5 0.533 0.240 0.421 

0.6 0.514 0.251 0,532 

0.7 0.481 0.273 0.640 

0.8 0.372 0.292 0.753 

0.9 0.370 0.304 0.802 

1.0 0.368 0.321 0.882 
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Zu diesem Problem sehen unsere Unterprogramme wie folgt aus: 


Tabelle 3.4 
Unterprogramme zu — 2 — 


fıtu) Ausgabe: 
2 FE LEL 7% LEL 240 Te 254 43 RCL 
2 IE A! 13 0° 241 1: 255 09 09 
2 ol 41 033 E42 43 : 256 99 PRT 
2 a5 = 35 = 26 268 25? 43 RCL 
si Res 1 RrS 33 PRT 258 27 327 
94 PRT 33 PRT 43 RCL 259 99 PRT 
32 RTH 35 AIY 15 15 260 43 RCL 
32 RTH 39 PRT 261 28 28 
43 RLL 262 99 PRT 
ar al 14 14 263 43 RCL 
= LEL REN. 39 PRT 264 29 29 
IB r& LEL 43 RCL 265 99 PRT 
eo 8 13 I 13 12 266 98 ADV 
ER = 32 RTH 92 PRT 26? 92 RTN 

a 


Tritt, wie in diesem Fall, die Zeit in den Funktionen als Parameter auf, ist SBR D’ ein ‚blindes’ 
Unterprogramm. Es muß aber auf jeden Fall existieren. Die Daten werden in dem jeweiligen Unter- 
programm aufgerufen (A’, B’ und C'). Die Werte ergaben: 


Grundeingaben: O, 19 Dateneingabe im LU. -% Ausgabe: 0. U 
DO. 1 -At Unterprogramm: U. -Vj 0. x 
oO 02 D. -Y 
0. -z 
0. yet 
0. -5 
D. -% 
0. 3 
0.37 
0.304 
0.802 
0.9 
0,37 
0. 304 
0.802 


1. 244885503 
. 5048762225 
12. 48759384 


. 1409574404 
1.495744 
ü, 


0,2159 DATE 0.368 
0.18 0.292 0.321 
a lli 1,753 0.882 
[u RE) lı 

ul. 372 0. 368 

I. 292 0.321 

(753 0.882 

0.5 1. 

1. 194439755 1. 32669626 


1.214064 5181075724 
r.6223336 3. 140083694 


Die graphische Auswertung ergab das in Bild 3.9 dargestellte s-t-Diagramm, das in Bild 3.10 darge- 


stellte v-t-Diagramm und das in Bild 3.11 dargestellte a-t-Diagramm. 


s{m) 


Bild 3.9 
s-t-Diagramm 


t(s) 


vim/s) 


Bild 3.10 
v-t-Diagramm 


alm/s?) 


Bild 3.11 
a-t-Diagramm 


t(s) 
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4 Kinetik 


Die Kinetik befaßt sich, wie bereits einleitend erwähnt, mit den Bewegungsänderungen, die ein 
System unter Einwirkung von Kräften erfährt. Dies wird zunächst am Massenpunkt und dann am 
starren Körper betrachtet. 


4.1 Kinetik des Massenpunktes 


Unter einem Massenpunkt versteht man einen mathematischen Punkt, in dem man sich idealisiert 
die Masse eines Körpers vereinigt denkt. Er besitzt mithin eine endliche Masse aber kein Volumen. 
Dies läßt sich bei den Ansätzen jedoch nicht immer ganz verwirklichen. 


4.1.1 Freie Bewegungen eines Massenpunktes im widerstehenden Mittel 


Die Bewegungskurve eines freien Massenpunktes unter Schwerkrafteinfluß ist als schiefer Wurf be- 
kannt. Betrachtet man einen schiefen Wurf im lufterfüllten Raum, wirkt der Bewegung des Massen- 
punktes, zu jedem Zeitpunkt seiner Bahn, eine Widerstandskraft entgegen. Nimmt man in erster Nä- 
herung an, daß die Widerstandskraft dem Quadrat der Geschwindigkeit direkt proportional ist, ge- 
hen Ansichtsfläche A und Luftdichte (Mediumdichte) 8 in die Gleichung mit ein, so erhält man 


BE wa. (4.1.1) 


Der Proportionalitätsfaktor, die Konstante c,, wird als Widerstandsbeiwert bezeichnet. Er richtet 
sich nach der Form des Massenpunktes und wird experimentell bestimmt. So ganz von Dimensions- 
losigkeit des Massenpunktes können wir hier also nicht sprechen. Einige der häufigsten Widerstands- 
beiwerte gibt Tabelle 4.1 wieder. 


Tabelle 4.1 Widerstandsbeiwerte 
Form Cv 


Zu jedem Zeitpunkt seiner Bahn, ist der Massenpunkt damit den in Bild 4.1 dargestellten Kräften 
ausgesetzt. Da der Luftwiderstand stets der Bewegung entgegengesetzt gerichtet ist, ergibt sich für 
die einzelnen Komponenten 

1 


Eiae zw 5 Av? cosa (4.1.2) 


Be 
2 


wöAv,v (4.1.3) 
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Bild 4.1 


Ballistische Kurve 


Fus= -Icw5Ay, yı“ +W (4.1.4) 
und 
Fu -iwöAn, v+v. (4.1.5) 


Für den Bewegungsansatz wird damit 


1 
mau wdAu Via tv (4.1.6) 
dv 
mai Io dAn Var (4.1.7) 
5A 
um v. v?+w dt (4.1.8) 


1 

may==zcwdAn, Vi tw omg (4.1.9) 
wöA 

dv, --( a War ro)a. (4.1.10) 


Wir erhalten zwei gekoppelte Differentialgleichungen, (4.1.6) und (4.1.9), deren analytische Lösun- 
gen nur schwer zu integrieren sind. 
Nach der Euler-Cauchy-Methode gilt näherungsweise für die allgemeine Differentialgleichung 


Inte (4.1.11) 
die numerische Darstellung 


Av 
un. (4.1.12) 
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In unserem Fall 


5A \ 
Av, = -(® vv“ +2 ) Ar 


exakte 


(4.1.13) Lösung 
und numerische 
Cy 5A Vo 
Ay = -( ET war? +g) At. hi 
> ) 
3% 
(4.1.14) RR 
Bild 4.2 = 
v-t-Diagramm 
Y. 


Unter Betrachtung der funktionalen Verhältnisse nach Bild 4.2 folgt, daß wir uns mit fortschrei- 
tender Berechnung immer mehr von der exakten Berechnung in zwei Richtungen entfernen. 

Die Gleichungen (4.1.13) und (4.1.14) sind die Grundgleichungen des in Bild 4.3 folgenden Berech- 
nungsprozesses. Es folgen die Tabellen 4.2 und 4.3 mit Speicherbelegung und Rechnerprogramm. 


Um nicht noch einen Wert im Programm eingeben zu müssen, wird das Programm einfach durch 

R/S abgebrochen, wenn die benötigten Daten errechnet sind. 

Außer in der Hinsicht, daß Luftströmungen, Erdrotation und Luftdichteschwankungen unberück- 
sichtigt bleiben, ist unser Berechnungsprozeß unter Auslassung eines weiteren Parameters idealisiert. 
Der Auftriebskraft durch das, den Massenpunkt umgebende, Medium, Solange sich die Bewegung 

in Luft vollzieht, ist die Auftriebskraft vernachlässigbar klein. Im Wasser z.B., muß sie jedoch Be- 
achtung finden. Da sie der Gewichtskraft der vom Massenpunkt verdrängten Mediummenge ent- 
spricht und stets der Gewichtskraft des Massenpunktes entgegen wirkt, ergibt sich unter Betrachtung 
von Bild 4.4 nur eine Beeinflussung der vertikalen Komponente 


may == 5 0 im Auy v2 +v2 -mg-F.. (4.1.15) 


Die darin enthaltene Auftriebskraft F, ist 
Fa = VK 5m 9. (4.1.16) 


Der Index k steht für Körper, m für Medium. Das Körpervolumen des Massenpunktes ergibt sich 
aus 


m 
Ya (4.1.17) 


Tabelle 4.2 Speicherplatzbelegung 


00 8Ac,,/2m 04 c„ 08 vu 12 t; 16 Atyorn. 
O1 m 05 At 09 vi 13 Av, 179g 

02 A 06 v; 10 x; 14 Ay, 18 vu; +v2 
03 5 07 a; 11 y; 15 Aton 19 
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LBLA/1. Start 


LBLB /2. Start 


Eingabe: 
%o. Yo Aton dient zur Angabe der 
Ag, Vo Schrittweite At für den Ausdruck 
to, At der errechneten Werte 
Alor 


Ya = Ya. + Avyu 


Wi=Wi-ı + Aw 


Bi 4.3 


Flußdiagramm zur Ballistischen Kurve 


Der Wert t; wird für jeden 
Berechnungsumlauf kurz angezeigt, 
und zwar unabhängig vom Ausdruck 
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Tabelle 4.3 
Programm Ballistische Kurve 


1. Start Ay,Vj 
on 7& LEL 38 94 +7- 145 42 RCL 
oo 11 A 0939 95 = 146 09 09 
002 47 CMS 100 42 570 147 65 x 
0c2 04 4 101 13 13 148 43 RCL 
Dü4 42 57T 102 44 SUM 149 05 05 
105 cıÜ nie] 103 us 08 150 95 = 
DoE 91 R-S 151 44 SUM 
00799 PRT Avyyi. Yyi 152 11 11 
os FE ST# 104 43 RCL 
u n 7 N Lu 105 00 00 ti 
ll 97 DEE 106 65 ® pr 
ou 0 107 42 RCL IS ErRe 
ac a 102 092 09 155 44 SuM 
Een 103 85% 156 12 12 
a 110 42 REL 15? 44 SUM 
= RC 111 12 18 SE 
a4 14 112 85 + 158 16 16 
ss & 112 43 R. EL Abfrage auf Ausdruck 
ee 159 43 RCL 
& lie. 08 x te0 15 15 
43 RC a ENER 1el 32 KT 
Er - 162 43 RL 
DE oc 118 05 05 163 16 16 
02 2 en 164 66 PAU 
En 121 42 sT0 Ice 7 
D1 Ic 123 44 sum 167° 43 RCL 
= >TO TE 124 09 09 Ausdruck 
Een = EN 
03% org 
3. cr73 Leon 293 170 10 10 
02 8 ca en 171 99 PRT 
oil 1 ESS 1 138 172 43 RCL 
42 570 NF-F 123 cs 173 11 11 
ıF ır ) 130 a ; 174 99 PRT 
us 131 37 Pr 175 43 RCL 
2. Start ü 132 42 sı 176 07 07 
= 5 n cz 177.99 PRT 
na37 76 re irE 43 RCL 
028 12 134 5 179 06 06 
033 81 au 1a ‚ee 180 99 PRT 
040 93 43 RÜL ie 181 43 RCL 
041 42 aa mi Be 182 12 12 
042 05 Er ER 183 99 PRT 
[NE = 51 4: 127 42 REL ie4 00 0 
044 99 08 128 0808 155 42 STO 
453 42 65 132 69% 186 16 16 
üdE 0. 43 140 43 RCL 18? 61 6T0 
4? 32 18 idi 05 05 158 42 RCL 
048 91 &5 142 95 = 
043 92 43 142 44 5UM 
casa 4z (iS 144 10 10 


70 


Damit ergibt sich in (4.1.15) eingesetzt 


ö 
May = 0 im Any N Fmaomyn 


(4.1.18) 


dm 
„ale 5-RN v2+rw +mg (1-5®). 


2 K 


(4.1.19) Bild 4.4 Sinkender Massenpunkt 


Unter Definition einer reduzierten Erdbeschleunigung 


6 
ar ( 2) (4.1.20) 
% 
erhalten wir 


wiederum die alte Gleichung (4.1.9). Damit kann das zuvor aufgestellte Programm auch für diese 
Fälle Anwendung finden. 
Nach Eingabe der Grundwerte c,,,ö, A und m durch Programmteil A, erscheint in der Anzeige 9.81 
und weist damit darauf hin, das die Erdbeschleunigung im Speicher 17 abgelegt wurde. Sie kann nun 
mit 

9 STO17 
korrigiert werden. Sie können sich auch noch ein zusätzliches Hilfsprogramm schreiben, das g' be- 
stimmt und in Speicher 17 ablegt. 
4.1.2 Das mathematische Pendel als erzwungene Bewegung 


Unter einem mathematischen Pendel versteht man die idealisierte Aufhängung eines Massenpunktes 
an einem gewichtslosen und unelastischen Faden (Fadenpendel). Durch Auslenkung aus seiner sta- 
bilen Gleichgewichtslage voltführt das Pendel unter Vernachlässigung von Luftwiderstand, etc. eine 
schwingende Bewegung um diese Ruhelage. Nach den Ausführungen unter 3.2 der Kinematik hatten 
wir für die Bewegung eines Massenpunktes auf gekrümmter Bahn eine Tangential- und Normalbe- 
schleunigung gefunden. In unserem Fall ergeben sie sich aus dem Ansatz 


ma =-mgsiny (4.1.22) 
und 

man = $S- mg cosy. (4.1.23) 
Für eine Kreisbewegung ergibt sich die Tangential- und Normalbeschleunigung aus 

,=1@ (4.1.24) 
und 


a = 1. (4.1.25) 
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Bild 4.5 
Mathematisches Pendel 


Eingabe 


l,m, Yo, wo 
to, At, Ator 


g 
Aw; == Tsinlpi.,) At 


= W.1+Aw;, 


pi =pi-ı + äp 


Bild 4.6 
Flußdiagramm zum math. Pendel 


72 


Mithin erhalten wir eingesetzt 


mi@=-mgsinyp (4.1.26) 
=" sing (4.1.27) 
dw g 

dw__I_ (4.1.28) 

Fr 

und 

mid? = S-mgcosyp (4.1.29) 
S= m(l9? +g cosy) (4.1.30) 


liefert die Seilkraft. Nach der bewährten Euler-Cauchy-Methode folgt 


g 
Aw = 7 sin (p) At (4.1.31) 
( (&) ) 
S=miIl At +9cosp). (4.1.32) 


und damit der in Bild 4.6 dargestellte Berechnungsalgorithmus. 


Tabelle 4.4 Speicherplatzbelegung 


00 Zähler 03 t 06 | 09 Ay 
01 Aton 04 w 07 m 10 s 
02 At 05 y 08 g 11 Atom. 


4.1.3 Gravitation und Planetenbewegung, als die Bewegung eines Massenpunktes 
unter Einfluß einer Zentralkraft 


Aus den von Kepler gefundenen Gesetzen über die Planetenbewegungen leitete Newton ein Gesetz 
über die Kraftwirkung zwischen Sonne und Planet ab. Dies hat allgemeine Gültigkeit für die Wech- 
selbeziehung beliebiger Massenpunkte und ist unter dem Begriff Gravitationsgesetz bekannt. Danach 
schreiben wir jedem Massenpunkt ein Gravitationsfeld zu, in dem dieser andere Massen anzieht. 


vı 


Bild 4.7 
Planetenbewegung 
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Tabelle 4.5 
Programm Mathematisches Pendel 


Start/Eingabe 


nano 65 05% 05 05 08? 43 RCL 
cc1 42 060 39 CDS 088 11 11 
onz [| 05 061 95 = 089 66 PRU 
23 [1 235 ÜbE 65 x 090 22 INV 
aud4 EI KM 063 42 RÜL 091 77 GE 
cos 42 RÜL 64 07 r 092 43 RCL 
u0£ DE 08 065 95 = 
NOT a5 = 066 42 STD Ausdruck 
INES} 44 SuM 067 10 [2] 
OC9 c4 04 Ü6S 32 MIT 02 2 nal 
210 SE R (69 65 x 055 05 05 
o11 oo Pi Si 07 01 1 096 99 PRT 
012 00 039 43 RCL "1 08 8 097 43 RCL 
212 oe 40 04 04 ürz 000 098 04 04 
o14 DW O4l 65 % 073 55 + 099 99 PRT 
c1S 3 ü43 453 RÜL ür4 83 0 100 65 x 
Jule R 43 02 02 07 95 = 101 43 RCL 
c1r 5 044 95 = 76 44 SUM 102 06 06 
118 1 045 59 - ir? 05 05 103 95 = 
219 STO 046 32 MIT R 104 99 PRT 
DE 03 04? 43 RCL i 105 43 RCL 
48 02 032 Ort 43 REL 106 10 10 
Start/Berechnung 49 35 = 079 032 03 107 939 PRT 
Aw, wi 050 33 X vsa 44 SUM 108 43 RCL 
us? 43 RLL üsz 44 SUM 110 99 PRT 
uüS3 08 0% 083 11 11 111 000 
eh .s x 112 42 STO 
[= = RCL Abfrage/Ausdruck 13 ı 11 
se 0208 084 43 REL 114 61 GTO 
057 689 8 085 01 01 115 43 REL 
SE 43 RÜL 086 32 KIT 


Wir wollen für unser Programm den umgekehrten Weg gehen und über ein Zweikörperproblem zu 
den Planetenbewegungen kommen. Nach dem Gravitationsgesetz herrscht zwischen den Massen m; 
und m» die Gravitationskraft 


F=fm m; et (4.1.33) 


Die uarin enthaltene Gravitationskonstante beträgt f = 6.67E-11 km? kg! s?2. Nach dem d’Alem- 
bertschen Prinzip erhalten wir für den Massenpunkt ma 


F=fmım; a = maa. (4.1.34) 
r 
Und daraus die infinitesimale Geschwindigkeitsänderung 


dv = fmı Sat. (4.1.35) 
r 


Nach altbewährtem Muster folgt 


Av = tm, Lat. (4.1.36) 
r 


4 


Bild 4.8 
Geschwindigkeitsplan 


Für die Geschwindigkeiten ergeben sich unter Betrachtung von Bild 4.8 und der Tatsache, daß es 
sich hier um ebene Bahnen handelt, die Beziehungen nach dem Cosinussatz 


v3 = Vya + Av? - 29 Av cos(a- pl. (4.1.37) 


Nach dem Sinussatz folgt weiterhin 


A 
Aa = arcsin (2 sin(a-y). (4.1.38) 
Durch erneute Anwendung des Cosinussatzes ergibt sich damit, nach Bild 4.7, auch die neue Ent- 
fernung 


rı = vra + Ar? - 2rg Ar cos (180 + p- a). (4.1.39) 


Und letztlich ergibt eine nochmalige Anwendung des Sinussatzes den in der Zeitdifferenz At über- 
strichenen Winkel 


Ay = arc sin (* sin 180 +9-0)). (4.1.40) 


Damit liegt unser Berechnungsalgorithmus in Bild 4.9 fest. Wie immer folgen Speicherplatzbelegung 
und Rechnerprogramm in Tabelle 4.6 und 4.7. 

Eine Diskussion der Keplergesetze erfolgt in Anwendungsbeispiel 4.1.5-5-. 

Das Programm wird auch hier wieder aus Einfachheitsgründen, nach Beendigung der gewünschten 
Ausgaben, durch R/S gestoppt. 

Dieses Prograrnm läßt sich auch leicht zu einem Mehrkörperproblem umfunktionieren. Es müssen 
lediglich die Gravitationskräfte der Massen beachtet werden. Deren Komponenten bezüglich eines 
gewählten Koordinatensystems sind dann zu bestimmen, ähnlich dem Programm Raketenbewegung. 


Tabelle 4.6 Speicherplatzbelegung 


00 Zähler 04 v 08 m; 12 @-w/180 +y-a 
01 Aton 05 « 09 m, 13 Ar 

02 At 06 y 10 f 14 Aa 

03 t 07 r 11 Av/F 15 Atyorh. 


Wird eine Masse tangential zur Erdkrümmung von der Erde abgeschossen, so ist die Startgeschwin- 
digkeit für ihre Flugbahn grundlegend. Aus dem Energiesatz, unter Vernachlässigung des Luftwider- 
standes, ergeben sich die in Bild 4.10 dargestellten und nur für die Erde zutreffenden Grenzgeschwin- 
digkeiten. Danach ist eine Mindestgeschwindigkeit von 7.9km/s notwendig, um die Erde umrunden 
zu können. Man spricht von der ersten kosmischen Geschwindigkeit. Die zweite kosmische Ge- 
schwindigkeit von 11.2 km/s erlaubt das Verlassen des Gravitationsfeldes der Erde. 
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LBLA/Start 


Eingabe 


f, m}, Ma, To, Yo 
ao At, to, Vo 


Aa; =arcsin(..... ) 


vo > 11,2 km/s 


p=ypi.ı+Yp Hyperbelbahn 


y 


Vo = 11,2 km/s 
Parabelbahn 


Men 
SE WW 
ri Vo 7 N N 


t=ti., +&%t 


Ausgabe 
9,7, Vi, 
Ft; 


Bein.) 
Bild 4.9 vo = 7,9 km/s 


Kreisbahn 


Fiußdiagramm zur Planetenbewegung 7,9 kmis <vo < 11,2 km/s 
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Ellipsenbahn 
Bild 4.10 


Flugbahnen tangential abgeschossener 
Erdsatelliten 


Tabelle 4.7 
Programm Planetenbewegung 


Start/Eingabe 

000 76 LBL 43 43 REEL 0935 65 = 142 

001 11 A n4a 05 05 096 52 & 143 

002 47 CMS os0 75 - os? 01 1 144 

003 01 1 051 42 RÜL 033 08.58 145 

004 00 0 052 06 e 39 000 146 

005 42 STD 053 54. 100 85 + 

006 0060 00 054 42 STO 101 42 FÜL t 

007 91 R-S 055 12 12 102 06 06 14? 492 RCL 

008 9% PRT 056 39 CD5 103 75 - 145 12 FE 

009 72 ST 05795 = 104 42 RÜL 143 44 SuM 

010 00 00 058 34 TH 105 05 05 ec. 48 75 

oil 97 182 059 42 570 106 54 9 1Sı 44 sum 

012 00 00 060 04 0 107 42 STD 152 02 02 

013 00 00 18 = 12 Er Se = 

218 2 ans 1 ge E03 Abfrage/Ausdruck 
ill 34 f# 133 42 RCL 

Berechnung 112 42 5T0 154 01 01 

EM 113 07 07 155 

016 76 LBL 156 

017 43 RCEL Ay, vi 157 

018 43 RCL 114 e 158 

019 10 10 115 159 

020 65 # 116 160 eG 

021 43 RCL 117 1eEl 43 FEL 

022 09 07 113 5 

023 55 + 1132 43 FCL Ausdruck 

024 43 RÜL 120 12 2 EZ 22 

025 07 07 121 35 SIN 182 52 

026 33 #2 ar; 122 935 = 164 38 

027 65 x 074 43 FÜL 123 35 = 1E53 43 

028 43 RCL 075 04 04 124 22 INW lE5 06 

029 02 02 076 65 x 125 38 SIN IE? 93 PR 

030 95 = Or? 42 RÜL 126 44 SuM IES 43 RC 

031 42 STO are 02 02 12? 06 06 ie3 ur 0r 

u 073 95 = 170 3% PR 

080 42 STD a; 171.42 

vi osı 13 1: 128 43 RCL ir2 08 

033 33 #8 1292 14 14 173 39 

034 85 + ' 130 44 SUM 174 42 

035 42 RCL 082 yE 131 05 05 17304 

036 04 04 083 + 176 99 

03733 88 os4 = REL Fj 1 4a 

038 75 025 07 07 122 42 FÜL IrE 11 

039 02 2 USE EHE i33 1C 10 1793 aa 

040 65 % 087 75 - 134 89% 180 43 

041 43 RCL 88 02.2 125 43 RLL 181 03 

042 04 04 083 65% 136 03 09 152 99 

043 65 x 020 2 RÜL 137 65 = 152 00 

044 43 RCL 031 07 07 125 42 REL 184 42 

045 11 11 092 659. 1393 05 082 155 15 

046 65 * 093 43 RÜL 140 55 - 156 61 

047 53 X ca4 13 12 141 42 RÜL 127 43 
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4.1.4 Die Raketenbewegung, als die Bewegung eines Massenpunktes mit veränderlicher 
Masse 


Die Raketenbewegung stellt das einzige uns bekannte Prinzip dar, die Beschleunigung eines Fahr- 
zeuges, ohne zu Hilfenahme von Reibungskräften, durchzuführen. Daher funktioniert dieses Prin- 
zip auch im kräftefreien Raum. Die erreichbare Raketengeschwindigkeit ergibt sich nach dem Im- 
pulssatz aus der Geschwindigkeit und der Anzahl der Masseteilchen austretender Verbrennungs- 
gase. Die Rakete erhält unter Betrachtung von Bild 4.11 das Impulsdifferential 


(m} + ma) dv = = v,dm. 
d 
a1) N 
Die Integration liefert 

Ve mı mq + m2 

"dv. _dm  _ [ dm 

Im m#+m | mı+m' 

va mim m Bild 4.11 IE 
(4.1.42) Raketenprinzip I 


mit dem Index e für End- und a für Anfangszustand. Unter der Voraussetzung konstanter Aus- 
strömgeschwindigkeit v, folgt 


1v1® = Tin (m) 1" "2 (4.1.43) 
s Va mq 


1 mı +mz my 
Mr (ve-v,) = In gen In v + m: (4.1.44) 


Mit der Annahme, daß die Startgeschwindigkeit v, = O ist, z.B. bei einem Start von der Erde auf 
einer Rampe, folgt 


er e) (4.1.45) 


Diese Formel ist für die Auslegung einer Rakete von großer Bedeutung. Bei Beachtung der kon- 
struktiv vorliegenden Ausströmgeschwindigkeit und der Raketenendgeschwindigkeit, als Vorausset- 
zung des Bahnverlaufs (siehe 4.1.3), ergibt sich unter Vernachlässigung von Strömungswiderständen 
der Nutzlast/Brennmasse-Quotient. 

Die Bewegungsgleichung der Rakete erhalten wir nach Bild 4.12 unter Einfluß der Gravitationskraft 
zum Zentralkörper M und der Impulskraft des Antriebes aus 


mas GT u + Mmz. (4.1.46) 


In Komponentenschreibweise heißt dies 


mau, cosa + fMm 5 c0s(180 + 4) (4.1.47) 
und 
d 
may = GW sina+fMm-5 sin(180+ y). (4.1.48) 
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Bild 4.12 
Raketenbewegung 


Unter der Annahme, daß die Menge der in der Zeiteinheit dt erzeugten Verbrennungsgase konstant 
ist, also 


dm 
Freie; u= const. (4.1.49) 


und damit für die Raketenmasse 
m= mg - ut (4.1.50) 


gilt, folgt nach altern Prinzip 


UV, cos“ cos (180 + y) 
Av, = (Ferm =G—) (4.1.51) 
mo — ut r 
und 
A (m were) 
= | —+ — 
wlmo-ut 2 (4.1.52) 


Damit liegt unser Berechnungsprozeß wiederum fest und die Programmierung folgt wie in den vor- 
angegangenen Ausführungen. In Bild 4.13 folgt das Flußdiagramm und in den Tabellen 4.8 und 4.9 
Speicherplatzbelegung und Rechnerprogramm. 


Tabelle 4.8 Speicherplatzbelegung 


00 Atyom 04 v 08 u 12 f 16 x 

01 Aton 05 a 09 v, 13 v, 17 F 

02 At 06 r 10 m 14 v„ 18 Zwischen- 
03 t 07 y 11M 15 y 19 werte 
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LBLA /Start 


Eingabe 
f,M,m,v,u,r 
y0,v,t, At 


UV, cosa@ == 
RER t 


Av ‚ 
= mo ut 


uv,sina sin (...) 
EURER ) Ar 


Av ‚-( 
y: mo - ut 


Die Umrechnung in 
«, v erfolgt durch den 
Befehl INVP>R 


Die Umrechnung in 
y, r erfolgt durch den 
Befehl INVP>R 


Bild 4.13 
Fiußdiagramm zur Raketenbewegung 


80 


Tabelle 4.9 


Programm Raketenbewegung 


Start/Eingabe 
000 76 LEL 
001 11 A 
002 47 CMS 
003 01 1 
004 02 2 
005 42 STO 
006 00 00 
007 931 R-S5 
008 99 PRT 
009 TZ ST* 
010 00 00 
011 97 DSEZ 
012 00 0ü 
013 00 ou 
014 06 06 
015 48 AI“ 
Umrechnung 

016 43 RÜL 
017 068 05 
018 32 %1T 
019 43 RÜL 
020 07 07 
021 37 PeR 
022 42 STD 
023 15 15 
024 32 X1T 
025 42 STO 
026 16 16 
027 42 RCL 
028 04 04 
029 32 44T 
030 43 REL 
031 05 05 
032 37 P-R 
033 42 STD 
034 13 13 
035 32 81T 
036 42 ST 
037 14 14 
Berechnung 

Avyi, Vi 

038 Te LEL 
039% 43 RÜL 
040 43 RÜL 
041 08 05 
042 65 x 
043 42 RÜL 
044 093 0% 
045 55 + 
046 43 RÜL 
047 10 10 
048 65 %* 
049 42 STD 
050 19 19 
051 43 RCL 
052 05 05 
053 39 CDS 


054 85 


059 
56 
057 
(38 
059 
a7] 
061 
(62 
063 
Ded4 
063 
Ü6E 
le 
LEE 
067 
rn 
ar 
ur 
073 
org 
AFS 
[TE 
Fr 
TE 
Ira 
nT=Jn] 
081 
DS 
Avyi Yy 


12 


C 
je BEN Ey Ed Kr Dad  ZSTIEE 5 


So 


-Ooocc 
(Ku ErEr Er Er enerErer, 


Media 


pin 
# Linn 


u DEET ENTER 


Gate fg 


ENT 


a a EN 


nn 


(Kr HE, Dee 


of 
U ENTE En 


ar} 


TO 
l = 


a 
u 
= 


1 


au 
Zi + 
Anm 


Lrı 
Ba | Bu 2 


1ll 85 
ill 42 
112 08 
114 935 
115 44 
1ls 16 
Ayuvi 
117 42 
112 12 
1193 69 
120 42 
121 08 
12295 
122 44 
124 15 
43 
LE 
a4 
te} 
43 
WE 
35 
44 
133 1u 
Umrechnung 
42 


 Bisfse Lin 


z 


43 
1 
: 
1 h 


[nz 28 


Loc Faie 


m 
fer 


163 69 8 
164 42 REL 
185 10 10 
lieb 55 + 
le? 42 REL 
1E2 02 06 
1ER 22 KE 
irü 95 = 
171 42 STO 
ır2 17 17 
tj 
173 42 REL 
174 02 02 
175 44 5UM 
Iire 03 05 
177? 44 5UM 
Ir8 00 00 
Abfrage/Ausdruck 
179 42 REL 
120 01 01 
181 32 817 
1852 42 REL 
18: 00 od 
124 E66 PAU 
189 22 INW 
186 77 GE 
157 42 RÜL 
Ausdruck 
185 25 AIV 
+2 FCL 
or Dr 
94 FRT 
4:3 FRÜL 
ü& U 
23 FRT 
4: FCL 
0005 
3 FRT 
43 FÜL 
a4 04 
Re 33 FRT 
& 42 RÜL 
E 10 10 
E 32 FET 
2 42 RÜL 
K ir 17 
& 93 FRT 
E 43 REL 
E 0203 
E s“ FET 
210 00 0 
2li 42 STO 
zl2 00 00 
213 &1 GTO 
214 43 FEL 


4.1.5 Anwendungsbeispiele 
Pe, DEN 


Ein Ball wird unter den Anfangswinkeln &, = 30° und &g = 60° mit der gleichen Anfangsgeschwin- 
digkeit vo = 20 m/s fortgeworfen. Seine Ansichtsfläche ist A = 0.01 m? und seine Masse m = I kg. 
Sein Luftwiderstandskoeffizient wird mit 0.4 angenommen. Die Luftdichte beträgt 1.293 kg/m?. 


A Eingabe: 27, 6065853894 
3. 334398783 
-18. 28376324 
13. 11594929 

1.6 


B Eingabe: 
21. 04341409 
1. 2863415927 

-24, 01681959427 
15, 20935427 

1.8 


Ausgabe: 


19, 67179151 
2 


37, 306.28321 
-1, 231955082 
-34. 02256163 
20. 6812462 

2.2 


229. 46521529 
1 7. 197018771 
r 91. 65762897 
4 
2 


me 
nd fa 
B 


vn 
Fi 
ı 
innen 


EREeER 
Omoıa 


Treginin 


15. 53903846 
3 


31. 399256536 
4. 349393918 
33. F2772996 
7, 12519042 
32 


I 


i 
Gaza 
"a 


33. 32321974 
1. 517744686 
-59. 10554626 
18. 72393803 
2.6 3.4 


35. 24466655 
-1. 895912267 

= 61. 93613172 
14. 1406747 20. 36683276 
= 3.6 


Fa in ac lee 
"Mina 


yim) Bild 4.14 
Ballistische Kurven 


Bild 4.14 zeigt das Ergebnis der Berechnungen. Eine in der Ballistik als Steilschuß bezeichnete Kurve 
(&o = 60°) erreicht die gleiche Weite wie ein ‚Flachschuß‘ (ag = 30°). Der Einfluß des Luftwiderstan- 
des ist hier gering sichtbar. Erst bei höheren Geschwindigkeiten weicht die Ballistische Kurve stärker 
von der Wurfparabel ab, 


Pe, 
Ein Fallschirmspringer, c,, = 1.5, A = 30m?, m = 100kg, fällt mit vo=0 im lufterfüllten Raum von 
5 = 1.293 kg/m?, aus 10000 m Höhe bei geöffnetem Fallschirm. 


Eingabe: 1. Cu 
1.2928 
30% A 
Loc, m 
GÜOL  Ar 
un vo 
U. ao 
Ü, xo 
000. yo 
Ü, to 

IE Alpre 


<R<x 


Natürlich wird dies kein vernünftiger Fallschirmspringer ausführen, denn er möchte ja nicht ‚stun- 
denlang’ obenbleiben. Bild 4.15 zeigt nämlich, daß die Sinkgeschwindigkeit sich einem Grenzwert 
nähert. Dieser läßt sich auch analytisch aus der Bedingung für den gleichmäßigen Fall 


ma Im5Ar (4.1.53) 
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[8] Ü, ‚ 0. 


9999, 243186 


2991. 2re4ll 
2, 
5. VERSE4ÜLE 
1.5 
Ü. 
3330, 712504 
-90. 
5. T94ecE2rd 
"] Ü. 0. 
3989, 558872 IIFT, 246843 Fan. 232081 
-30 -30 -30. 
5. 80071209 5. 508888723 5. 2068595113 
22 4,2 6.8 


0. 

3965. 171682 
-90. 

3. 806895117 


D. 

9964, 010303 
-90. 

5. 806895118 


Ö. 

99682. 848924 
-, 

5. 8065895119 
6.8 


D. 

996&1. 687545 
-90,. 

5. 80689512 


0, 

960, 526166 
-9,. 

5. 20689512 
7.2 


0, 

3259. 364787 
-90,. 

5. 80689512 
7.4 


berechnen, wurch v{m/s) 


/2mg 2) 
7 =]/ ——., (4.1.54) 
rn cw5A 5 


In unserem Fall also 
4 
2 100kg 9.81 mm? 
"= 21.5 1.293Kg 30m? a 


vo = 5.80689512°. 


Bild 4.15 


Freier Fall im Lufterfüllten 
Raum 


Unsere Näherungsberechnung entspricht 
damit der Rechnergenauigkeit. 


N En 
Ein Taucher c,, = 0.8, A = 0.2 m?, m = 100kg, ö« = 1500 kg/m? springt mit vo= 1 m/s ins Wasser 
von 5, = 1000 kg/m?. 

Entsprechend der Ausarbeitung beträgt die zu verändernde Erdbeschleunigung 


"Lac m /,_ 1000 
en ( 1500 


Wie Bild 4.16 zeigt, gilt auch hier 
nach der Formel (4.1.54): 


__/ 2100 kg 3.27 m m? 
ve” \ 208 1000kg 0.2m2 


vo = 2.021756662 ”. 


m 
= 3.2775. 
) Eu 


vim/s) 


Eine konstante Sinkgeschwindigkeit 
wird hier schneller erreicht. 


Bild 4.16 


Sinkgeschwindigkeit im 
Wasser 
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Eingabe: [a PR: 


Ausgabe: [u age 


ü. 
..zdededrgln 


1, 235765406 


IP 


l. 
-2. 479921861 

-3U, 
2. MIFIFFFET 


1.€ 


me rn<x 


[j 
DW 


2, 021240715 


2.4 


Er 

u 
-6. 11496559 
-30, 
2. 021719492 
3.2 


Ü. 
"5. 519211563 


6. FEZEERDILE 
-30, 
2. ÜELT4ERBE 
3. £& 


D. 
-7, 325010014 
-3n. 


u 

-9, 243 7ES187 
30, 

2.021 756469 
4.8 


u, 

9, "asl1649E 
-al, 

2, DE17IBES6E 
Te 


0. 

-10. 96717046 
-90,. 

2. 021756648 
3.6 


D. 

-11. 37152179 
-9,. 

2. 021756655 
5.8 


D. 

-11. 77587312 
-90, 

2. 021756658 


D. 

-12,. 18022446 
-90. 

2. 02175666 
6.2 


0. 
-12. 584537579 


-12. 25892712 
-90, 
2. 021756661 


-13. 79762978 
-30. 
2. 021756662 


0. 

14. 20138112 
-D. 

2. 021756662 
rT.2 


0. 

-14. 60633245 
-90,. 

2. 021756662 
7.4 


4 


Ein Fadenpendel (math. Pendel) mit einem 2 m langen Faden wird um y = 60° ausgelenkt. Der 
Massenpunkt hat m = 10kg. 


Eingabe: 10. m 


2. | 


60. yo 
0. wo 
0, to 

0.01 At 

0.1 Atprt 


Ausgabe: 


58. 66436294 » 
-. 4230498093 w 


-. 8460996196 V 
54. 24144219 $ 
0.11 


54. 9337676 
-,8351873851 
-1. 67037477 
69. 63890194 
0.2 


48. 9096734 
-1. 223665369 
-2. 447330738 

93. 51364141 
0.3 


40. 77579918 
-1. 572711235 
-3. 145422469 

122. 739927 
0.4 


30. 81587938 
-1. 863632016 
-3. 727264033 

152. 7613874 

0.5 


19. 42538461 
-2. 076683422 
-4. 153366843 

178. 0705392 

0.6 


7. 107079577 
-2. 194586747 
-4. 389173494 

193, 3807 
0.7 


-5. 556180858 
-2. 206647066 
-4. 413294131 
195. 2107304 

0.8 


-17. 34950204 
-2. 111722529 
-4, 223445055 
183. 13033423 
0,9 


-29, 48316851 
-1. 315559004 
-3. 827118008 

159. 9239933 
1. 


-39, 64261924 
17 


"130. 56862 
1.1 


-48, DISSETIT 


-1. 3059282395 
-2, 611096517 
100.6526364 

1.2 


-1. "aaa7 13053 
75. 04907206 


1.2 
-52. 33240312 
-. SIE2SIRZEIT 
-1. 032596521 
Sr, 2617r1lol 
1.4 


1 
4 
” 


-50, 44294069 
1. 139354525 
2. E74r09195 
87. 59302671 

1.5 


138. 2744336 


. 4701448906 
sl. E0872 


3.1 


ee 


. 141247724 


46. 

1: 

10”. 1022081 
2,7 


(grd) 


Bild 4.17 
Schwingungsverlauf des Fadenpendels 


Auch hier wollen wir eine Grenzwertbetrachtung anstellen. Aus der Energieumsetzung läßt sich 
analytisch die Geschwindigkeit des Massenpunktes im tiefsten Punkt seiner Bahn bestimmen. Nach 
Bild 4.18 wird das durch die Höhendifferenz h umgesetzte Energiepotential, durch die Gleichung 


v2 
mah=m, (4.1.55) 


beschrieben. Daraus folgt durch Umstellung 


max = V2gh (4.1.56) 
_1 2:9.81m-1m n En 
u ) ee 
2 VE Ymas | 
Ymax = 4.429446918 
max . s' Bild 4.18 Energiebetrachtung 
En 
Ein Satellit von der Masse m, = 1000 kg, soll auf einer elliptischen Bahn um die Erde (m, = 5.973E24 kg, 


f = 6.67E-20 km?/kgs?) gebracht werden. Dazu wird eine Anfangsgeschwindigkeit von vo = 10km/s 
gewählt. Der Satellit startet tangential (Yo = 0°, &o = 90°) zum mittleren Erdradius von rg = 6378 km. 


Die graphische Darstellung der Satellitenbahn zeigt Bild 4.19. Daran lassen sich nun anschaulich die 
drei Keplerschen Gesetze erklären. 
T. Gesetz 
Die Planeten bewegen sich auf elliptischen Bahnen, in deren einem Brennpunkt die Sonne steht. 
Dies heißt nichts anderes, als das die Strecken 

in, mit n=3,4,5,6,..., 


in Bild 4.19 gleich lang sind. 
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Eingabe: 6. 67- 


5.973 


Ausgabe: 


72. 20836605 
8654. 871392 
136. 3367104 
8. 236058417 
5. 318597327 

1000. 


106. 4798458 
12387, 89231 
161. 4244459 
6. 336350213 
2. 596112488 

2000. 


125. 1863725 
15777. 33403 
177. 763171 
5. 131882031 
1. 600483157 
; 3000. 


137. 6497628 
18626. 80871 
190. 8273146 
4. 314158425 
1. 148263182 

4000. 


2. Gesetz 


Die Flächengeschwindigkeit eines Punktes ist konstant. 


au f 
24 m 


-n<rne 


147. (427105 


z03E7, 42 
4 


0, PÜEZ0rES3 
scood, 


1 
> 


= 
v“ 
E 
> 
2 
2 
be 


45 
14000, 


In unserem Fall heißt dies, die schraffierten Flächen A, und A, stimmen überein. Dies stimmt für 
jedes Flächenelement der Bahn, das in der gleichen Zeiteinheit überstrichen wird. 


3. Gesetz 


Das Verhältnis der dritten Potenzen der großen Halbachsen der Bahnen zu den Quadraten der Um- 
laufzeiten ist konstant. 
Das heißt allgemein 


0 % 
bi by 


Dies läßt sich jedoch nur an mehreren unterschiedlichen Bahnen zeigen. 
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200fache Schrittweite 
der Berechnung 
oe 20fache Schrittweite 


der Berechnung 


Bild 4.19 
Satellitenbahn 


(Ein Größenvergleich, die durchschnittliche Entfernung zum Mond beträgt 384 403 km.) 


-6 — 

Eine Rakete soll eine Nutzlast von 50 kg auf eine elliptische Umlaufbahn um die Erde bringen 

{M = 5.973E24, ro = 6378 km, f = 6.67E-20 km?/kgs?). Eine Endgeschwindigkeit der Nutzlast 
von 10 km/s ist anzustreben (siehe 4.1.3). Die Ausströmgeschwindigkeit der Verbrennungsgase be- 
trägt 5km/s = const. mit einer Massenausströmung von u = 10kg/s = const. 

Aus Gleichung (4.1.45) läßt sich die Brennstoffrmasse der Rakete abschätzen 


ve 
ma=m; (e* z ı) (4.1.57) 
,ı1 
ma = 50kg (e? - ı) 
m; = 320 kg 
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Eingabe: 


6.67-20 
3.972 24 M 
370. mo 
EP vs 
10. u 
0. v0 
6378, ro 
45. u) 
0. vo 
0. io 
0.1 At 
1, Atprt 
Ausgabe: 


. 0005140347 
6378. 042354 
54. 56285634 
«1307393369 

360, 
3. 525705097 


‚0020305467 
6378, 137008 
56. 91772225 
. 265900 rl 

350, 
3. 427646148 


‚0046414744 
6378. 339497 
53. 32043401 
» 4053902567 

340. 
SESRIDERTE 


‚DOEIEEIOE4 
6378. 39952 
39, 31896623 
5493652359 

330. 
3. 231342395 


. 0132806971 
6378. 307162 
60. 08970832 
. 6980481876 

320. 
3. 133110874 


an3z3<R08% 


. DSEFSLFITF 
5391. 544297 
top 2-5 8 Kuletorc 
1. I3226E30272 
Zrü, 
2641377577 
10. 


50, 
2. 444512461 


ca 


210, >0. 
2, 0SOF1S753 1. .Errl20sd44E 
30, 


‚I2201021 36 ZSSE 
Babe, S93E15 3357 
Ed OHALSET4 Er 
2, 329104956 9. 1OREFEFFE Far 
zu, 1:20. Eu, 

1. 352285176 1. 264190308 SFUrFERdER 
17 24 21. 


55, 
5.50 


1.06511638 


FIR 
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v{km) 
150 


Erdoberfläche 


4 v..9,82 


- Beginn der Satellitenbewegung 


) 

| 
100 | h 

| 325 

I 

l 


50 


Bild 4.20 
Raketenbewegung 


6350 A 6400 6450 


Bild 4.20 zeigt graphisch die Flugphase. Nach 32 Sekunden ist der Treibstoff verbrannt und es be- 
ginnt eine Satellitenbewegung. Sie läßt sich mit dem Programm aus 4.1.3 weiterverfolgen. 


4.2 Kinetik starrer Körper 


Unter einem starren Körper versteht man idealisiert einen solchen, der bei Krafteinwirkung keine 
Formänderung aufweist. Diese Annahme hilft, viele Probleme ausreichend genau zu lösen. 
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4.2.1 Massenträgheitsmoment 


Wir haben bisher Probleme der fortschreitenden Bewegung (Translation) betrachtet. Dabei hatten 
alle Masseteilchen eines starren Körpers die gleiche Bewegung und wir hatten ihn deshalb als Massen- 
punkt eingeführt. Wir kommen nun zur Drehung des starren Körpers um eine feste Achse (Rotation). 
Die beschleunigte Drehung eines starren Körpers in der Ebene um eine feste Achse wird durch die 
Einwirkung eines Drehmoments M hervorgerufen. Dabei vollführt jedes Massenteilchen dm, nach 
Bild 4.21, eine beschleunigte Bewegung. Aus der Kinematik wissen wir, daß ein Massenteil auf ge- 
krümmter Bahn, einer Normal- und Tangentialbeschleunigung unterliegt. Dies führt nach dem 
d’Alembertschen Prinzip zu dem Ansatz 


dF, = dma, (4.2.1) 
und 


dF„ = dma.. (4.2.2) 


Bild 4.21 Beschleunigte Drehung 


Während das Normalkraftdifferential dF„ kein Drehmoment hervorruft, seine Wirkungslinie geht 
durch den Drehpunkt, ruft das Tangentialkraftdifferential dF, einen Drehmomentanteil von 


dM = rdF, (4.2.3) 
hervor. Für die Gesamtheit aller Anteile gilt damit 

M -[r ar. [tr dm a. (4.2.4) 
Mit 

a.=re. (4.2.5) 


Darin ist e die Winkelbeschleunigung, die für alle Masseteile gleich ist, also 

M=e fr dm. (4.2.6) 
Analog zur Massenträgheit bei der Translation, F = ma, bezeichnet man 

I - [r?am (4.2.7) 


als Massenträgheitsmoment eines starren Körpers. Genauer, da es sich auf eine Achse bezieht, als 
axiales Massenträgheitsmoment. Zu beachten ist, daß das Quadrat des Abstandes in die Gleichung 
eingeht. 


Tabelle 4.10 zeigt die Zusammenstellung der Massenträgheitsmomente einfacher Grundkörper. Auch 
komplizierter gestaltete Körper lassen sich mit den Gleichungen berechnen, da die Summe der Mas- 
senträgheitsmmomente einzelner Grundkörper gleich dem Massenträgheitsmoment des aus diesen be- 
stehenden starren Körpers ist. Dies liegt an der Eigenschaft des Integrals in Gleichung (4.27). Nicht 
immer fällt nun die Drehachse des Grundkörpers mit der des starren Körpers zusammen. Dazu be- 
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Tabelle 4.10 Massenträgheitsmomente (axiale) 


Massenträgheitsmoment 


Quader 


Hohlzylinder 


Hohlkugel 


Gerader Kegelstumpf ! 


A) 8 x 


Kreisring 
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(x, Yı) 


Bild 4.22 


N (x2, v2) 


trachten wir nach Bild 4.22 das Massenträgheitsmoment bezüglich einer zweiten Achse gegenüber 
der Schwerpunktsachse. Es gilt für den Radius r, die geometrische Beziehung 


12 = (x1 = X2)? + (yı= v2)? 14.2.8) 
=? -2xıXatxd +ty2 -2yı yatyZ (4.2.9) 
=r2 +2? 2(x1x2 + Yı Ya), (4.2.10) 


die eingesetzt in (4.2.7) 
Id [7 am +32 [am-2 [ixıxa+vı Ya) dm (4.2.11) 


ergibt. Darin ist f dm= m und - 2 (x1 %a + Yı Ya) dm das statische Moment des starren Körpers 
bezüglich des Schwerpunktes, also Null. Damit folgt die, als Satz von Steiner (Verschiebungssätz) 
bekannte Gesetzmäßigkeit 


Id, = Id, +ma2. (4.2.12) 


Nun können wir die in Tabelle 4.10 angegebenen Formeln programmieren und mit Hilfe des Steiner’ 
schen Satzes auf jede beliebige Achse umrechnen. Das heißt natürlich, daß die 2. Achse im Abstand a 
parallel zur 1. Achse verlaufen muß. 


Tabelle 4.13 Speicherplatzbelegung 


00 m 01 Zwischenspeicher 


Nachfolgend soll noch ein Programm zur näherungsweisen Berechnung rotationssymmetrischer Kör- 
per jeglicher Querschnittsform aufgestellt werden. Dazu betrachten wir Masseteilchen nach der in 
Bild 4.23 dargestellten Form und erhalten 


dm = 5 2nrhlr) dr. (4.2.13) 


In (4.2.7) eingesetzt 
r2 
Id=52r ron dr. (4.2.14) 


n 


95 


il Bei den LBL’'s treffen wir folgende Zuordnung: 


LBL__ /Start 


A: Quader 
B: Hohlzylinder 
Eingabe C: Hohlkugel 
benötigter D: Kegelstumpf 
Werte E: Ring 


R/S / Prgm. Ende 


Umrechnung nach Steiner für eine zweite Achse im 
Abstand a. 

Bedingung für den Start: 

Id} muß bei Programmaufruf in der Anzeige stehen. 


LBLE’/Start 


Eingabe 


a 


Id, = Idı + ma? 


"Ausgabe E: 
'd Bild 4.23 


Flußdiagramm zur Berechnung von 
R/S / Prgm. Ende 


Massenträgheitsmomenten 


Ersetzen wir hierin angenähert das Differential durch die Differenz, so erhalten wir 
n=rg-Ar 
ld=527 S rhtndar. (4.2.15) 
nen 
Diese Formel läßt sich noch genauer gestalten, wenn wir das arithmetische Mittel der beiden Höhen 
eines Streifenelements bilden 
n=n-Ar 


Id=52nr \ r> 


haln) + hafr, 
Bel DAIHE. (4.2.16) 


nen 
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Tabelle 4.12 
Programm Massenträgheitsmomente 


Quader: 
000 7&£ LEL O45 02 3 189 75 > 132 95 = 
001 11 A o4e 935 = ao 42 RÜL 134 98 AD“ 
002 91 RS O4? 32 84T ca9ı Cl 01 135 72% FRT 
003 #9 PRT [u 7 - Be u u = Per u 126 91 RS 
004 42 STD 049 42 FLL 032 032.23 
005 00 o00 cosco 00 00 ua4 95 = ER 
006 55 + 051 55 + 095 98 AMW Kreisring: 
007 01 1 0S2 02 2 Da6E 99 PET 137 T& LEL 
008 02 2 053 935 = a7 91 Res 138 15 E 
009 65 X 054 133 931 RS 
010 53 € 155 140 99 FRT 
011 91 RS 056 Kegelstumpf: 141 42 5T0O 
012 99 PRT 0S7 9% FR u = 142 00 nu 
013 33 x2 058 91 Res ass Te LEL 143 65 x 
014 85 + was 14 D 144 53 ( 
015 91 R#S en a Br 145 31 Rx 
016 93 PRT 93 FR + 
01? 33 x BER NU) 102 42 5STD r = 
018 95 = 059 FE LEL 102 00 50 193 85 
019 98 ADY el 12 C 104 €&5 1492 03 2 
020 993 PRT cell 91 RS 105 92. 150 ES x® 
021 91 RS 062 93 PRT i0e 02 3 re Ta 
063 42 STD 107 €5 * 152 33 PR 

: De4 00 Gl 108 53 % E23 2a WE 
Hohlzylinder: DES 65 i0o9 53 R nn ss 2 
022 76 LBL 066 92. 110 31 Re ıs5 04 4 
023 12 B Tr 111 933 PRT ISE 95 = 
024 9 RS NEE 85 ae: 1iz yS 8 157 aa Aw 
025 9° PRT 69 ad 113 32 81T 155 33 PRT 
026 42 STO nr0 52 & 114 05 5 153 91 RS 
Er g 00 ort 91 Res 115 75 - " 

= 72 939 FRET lie 91 Res j 

029 04 4 073 45 v® 11793 PRT ia 
030 65 x cr4 32 KIT 118 42 STO 1e0 FE LEL 
031 53 % ars 05 5 119 01 0 1E1l 10 E* 
032 91 RS Te 79 - 120 45 1E2 9393 FRT 
033 99 PRT off 91 Res 121 05.5 1E3 85 + 
034 33 x%8 O7 9% FRET 122 \ leE4 43 FÜL 
035 85 + 079 42 STO 123 jE5 00 an 
036 91 RS 080 0 01 124 16E& 939 PRT 
03799 PRT 081 4959 Y# 125 | U 
038 33 82 os2 03 5 126 1EE 1 RS 
039 75 - 083 54 127 ET ! 
040 32 KIT 084 55 + 128 1r0 
041 91 RS DES 53% 129 irı 55 
042 99 PRT ÜBE 32 ST 130 irz2 38 Al“ 
043 33 #2 DEr 121 173 933 FRT 
044 55 * nt=7=] 132 174 31 Fr 
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Angleichung nach (4.2.15) 


= 
I Angleichung nach (4.2.16) 


% LD 


ha (r,) 


Bild 4.24 


Bild 4.23 Masseelement 


LBLA /Start 


Eingabe 
ö (Dichte) 
ra, fa, hı1, Ar 


Anzeige 
ri 


Da h; {r,) =h, (r;4 4) ist, braucht immer 
| nur ha (rj) eingegeben werden. 


Eingabe 


Ausgabe 


h,+h 
amt 
r; Bu Ar 


Ausgabe 
Id 


R/S / Prgm. Ende 


Bild 4.25 

Flußdiagramm zur Berechnung des 
Massenträgheitsmomentes eines 
rotationssymmetrischen Körpers 
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Tabelle 4.13 Speicherplatzbelegung 


00 

016 

02 rıfr; 
03 r2- Ar 
04 Ar 

05 & 

06 h; (r,) 


Tabelle 4.14 Programm Massenträgheitsmoment eines rotationssymmetrischen Körpers 


Start/Eingabe: 
000 7r& LBL 024 De KE a4 “tt ja Tc} 
001 1i A ces 35 ADW NE Cr :TO dead 975 
002 47 CMS a4sS CE u DES 4: 
003 91 R-S Anzeige/Eingabe: Ode 54. DEE 04 
004 99 PRT us6 147 55 + 067 44 ZUM 
005 42 STO n>7 [1 EP 1 Por DES O 0E DE 
006 01 01 038 a4 35 = 163 E11 GTO 
007° 91 RS 123 050 93 FR [ro 43 REL 
008 9° PRT n®H cS1l 98 AIW 
009 4: STO 031 Endausgabe: 
010 02 02 032 Summe: o?1l Fe LEL 
oil 91 RS 33 ÜS2 BSH IT2 93 PRT 
012 39 PRT c34 052 42 RÜL Or2 42 RÜL 
013 75 - (35 939 PR oS4 04 04 74 05 05 
> 055 99 = 75 65 © 
015 993 PFRT aSE 44 5UM (TE 43 RÜL 
016 42 STO Berechnung: üsr 05 08 a 1 5 
017 04 04 USE OO SECHIT ME ES HM 
018 95 = ar 9300 = Abfrage/Ende: [u Dr a Po 
019 42 5TO DE FE juT-T u Be Ep 
020 03 03 032 520% 058 BR ga 
021 91 RS 040 42 RÜL [59 ag = 
022 99 PRT dl 02 08 ja1=20] a2 ADW 
023 42 STO ade 85 + ost 3 FRT 
ÜBE [np Al RS 


4.2.2 Das physikalische Pendel 


Wird ein Körper außerhalb seines Schwerpunktes drehbar gelagert, siehe Bild 4.26, so führt er, 
unter Auslenkung aus seiner stabilen Lage, Schwingungen ähnlich denen des Fadenpendels durch. 
Es gilt analog 


Ide=-mglsinp (4.2.17) 
dw mgl 
dt Ta sine. (4.2.18) 


Nach der Euler-Cauchy-Methode folgt 


Aw = - —— sinp (4.2.19) 
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die Differenzengleichung der Bewegung. Wenn wir diese mit der des Fadenpendels vergleichen, Glei- 
chung (4.1.31), und eine reduzierte Pendellänge, der Art 


Id 
Ira = (4.2.20) 


einführen, ergibt sich 


eek 122 
Ired 


Bild 4.26 
Physikalisches Pendel 


eine Übereinstimmung. Es ist also kein neues Programm notwendig, sondern wir können unter Be- 
nutzung der reduzierten Pendellänge nach (4.2.20) das Programm Fadenpendel benutzen. 
4.2.3 Reduzierte Masse und Schwungmoment 


Betrachten wir noch einmal Gleichung (4.2.7). Darin ist der Radius r abhängig von der Lage des 
Massenelements dm. Stellt man sich nun die gesamten Masseteilchen dm auf einem konstanten 
Radius vor, etwa wie bei einem Zylinder mit sehr geringer Wandstärke, so läßt sich schreiben 


ld = r? (am. (4.2.22) 


Die Größe f dm ergibt sich damit für ein bekanntes ] und einen beliebig gewählten Radius r zu 
[em en (4.2.23) 


und wird als reduzierte Masse bezeichnet. Mit der Annahme, daß die reduzierte Masse der Masse 
des Drehkörpers entspricht, ergibt sich nach (4.2.23) ein bestimmter Radius 
Id 

i= R m: (4.2.24) 
der als Trägheitsradius bezeichnet wird. Unter Definition eines Trägheitsdurchmessers 

D;, = 2i, (4.2.25) 
folgt eingesetzt 

mD? = Ald. (4.2.26) 
Diese Größe wird allgemein als Schwungmoment bezeichnet. Danach können wir das unter 4.2.1 auf- 


gestellte Programm um die in Bild 4.27 dargestellten Programmteile ergänzen. 


Tabelle 4.15 
Speicherplatzbelegung 


00 m 


100 


yes | Berechnung der reduzierten Masse 


LBLA’/Start 


Bedingung für den Start: 


j 
| 
| 
| 
| 
I 
u — 


Id muß bei Programmaufruf in der Anzeige stehen 


A| Berechnung von Trägheitsradius und Schwungmoment 


Bedingung für den Start: 
Id muß bei Programmaufruf in der Anzeige stehen. 
m wird aus Speicher 00 gelesen 


Bild 4.27 
Fiußdiagramm zu reduzierte Masse/Schwungmoment 


Tabelle 4.16 Programm reduzierte Masse/Schwungmoment 


reduzierte Masse: 


175 76 LBL 124 90 PRT da RÜL 
176 i6 AR" 185 91 Res am Qu 
177.99 PRT “3 FRT 
‚ss ER RS Schwungmoment: ar 3 
180 993 PRT 1865 FE LEL = AD" 
181 33 3 rn ne a FRT 
182 95 = 185 99 FRT BSH 
183 98 AIY 129 59 = je 


4.2.4 Deviationsmomente 


Bei rotierenden, starren Körpern wird auch nach der auftretenden Lagerreaktion oder den inneren 
Spannungen gefragt. Dies führt zu folgender Überlegung. Wird ein, nach Bild 4.28, um die z-Achse 
rotierender Körper, plötzlich zu einer weiteren Achse, z.B. zur x-Achse, fest fixiert, so würde ein 
Masseteilchen des Körpers das Momentdifferential 


dM, = zdF, (4.2.27) 
bewirken. Der Kraftanteil ergibt sich aus 


dF,=dF„sinp+dF,cosp. (4.2.28) 


Bild 4.28 


Die zum Mittelpunkt der Bahn weisende Führungskraft Fn, hier Zentripetalkraft genannt, ergibt 
sich nach dem d’Alembertschen Prinzip aus der Normal- bzw. Zentripetalbeschleunigung. Diese 
wiederum nach den Ansätzen aus der Kinematik zu 


dF„. = a, dm = rw? dm. (4.2.29) 
Analog folgt für den Tangentialkraftanteil 

df, = a, dm = redm. (4.2.30) 
Für alle Masseteile damit 

M, -)z dF, = [2r@2am sinp+ [zream cos p. (4.2.31) 
Da w und e für alle Teile gleich, folgt 

M, = w? [zv dm+e [2x dm. (4.2.32) 
Darin bezeichnet man allaereih analog zum Massenträgheitsmoment, 


Id. = [xzem (4.2.33) 


als Deviations- oder Zentrifugalmoment. Der Sonderfall 


Idun = |xam (4.2.34) 


wird als polares Trägheitsmoment bezeichnet. Analog zum Miassenträgheitsmoment wollen wir auch 
hier ein entsprechendes Programm für die in Tabelle 4.17 dargestellten Deviationsmomente einfacher 
Grundkörper aufstellen. Diese wenigen sollen als Beispiel genügen. Analog zum Steinerschen Satz 
gilt auch hier, siehe Bild 4,29, 


Idog = Id,, + mab. (4.2.35) 


Da die Programme den gleichen Aufbau wie Bild 4.23 haben, erspare ich mir an dieser Stelle ein 
Flußdiagramm. 
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Tabelle 4.17 Deviationsmomente 


Deviationsmoment 


Kugeloktant 


Tabelle 4.18 
Speicherplatzbelegung 


00 m 


Bild 4.29 
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Tabelle 4.19 Programm Deviationsmomente 


Quader: Kugeloktant: Verschiebung: 
nan FE r6& 37 76 LEL 058 76 LBL 
IT u 12 038 12 C 059 10 E’ 
co2 51 si 0133 91 RS 060 993 PRT 
00393 Er 040 29 PRT 061 85 + 
oo4 48 gez N! a4l 42 5STD 062 43 RCL 
005 Go oo c4e 00 00 063 00 00 
IOE 34 a4 42 934 +0 064 65 x 
gar 55 65 44 85% 065 91 RS 
DS 04 “1 45 02 2 066 99 PRT 
c03 65 94 [u Tu EB 067 65 x 
io 81 2 ü4r 531 Res ü68 91 R/S 
coll 993 Sa 04S 939 FRT 69 99 PRT 
E65 [1 0149 33 #8 070 95 = 
13 91 [ne [SG 55 *+ 071 98 ADYV 
old 939 35 os 05 5 072 99 PRT 
01S 95 as USE 589 * 073 91 R/S 
oje 9 39 053 849 m 
al? 99 “1 cS4 95 = 
als 91 055 98 AD% 

0SE 93 FRT 

0157 91 Res 


4.2.5 Massenkräfte beim Kurbeltrieb und ihr Ausgleich 


Der Kurbeltrieb dient zur Umwandlung von Schub- in Drehbewegung und umgekehrt. Wir wollen 
an dieser Stelle speziell den Bewegungsablauf eines Kolbenmotors betrachten. Der algebraische 
Ausdruck für die reale Kolbenbewegung ergibt sich, unter Betrachtung von Bild 4.30, aus folgender 
Ableitung 


x=|+r-1lcosß-rsinp (4.2.36) 

r  sinß 
= Dr) (4.2.37) 
cosß = Y1 - X cos?y (4.2.38) 
x = r(1- sing) + III -V1-A? cos?y), (4.2.39) 


bzw. bei allgemeiner Phasenverschiebung um « 


x=r(1-sinip-a)) + 1(1-V1-%? co(p-a)). (4.2.40) 


Geschwindigkeit und Beschleunigung ergeben sich wiederum angenähert aus den Differenzenquo- 
tienten 


Ax 
vs At (4.2.41) 
und 
Av 
er (4.2.42) 
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Bild 4.31 Schubstange 


Bild 4.30 
Kolbenmotor 


Den Bewegungen der Triebwerksteile gemäß, werden oszillierende und rotierende Massen unter- 
schieden. Danach ergibt sich die oszillierende Masse aus 


fst 
mo = Ma) + me +tme: (4.2.43) 


Darin ist mg, der Massenanteil der Schubstange, der nach Bild 4.31 durch den Faktor rg,/l seinen 
oszillierenden Anteil hat, m« die Kolbenmasse und mg die Kolbenbolzenmasse. Die rotierenden 
Massenteile setzen sich aus 

Ir r 


— + mw = +mz+tmn (4.2.44) 


MAR 7 Mg 1 


zusammen. Darin ist mg, (1 rg,)/l der rotierende Massenanteil der Schubstange, my, die Kurbel- 
wangenmasse, die nach Bild 4.32 durch den Faktor ryy/r auf den Drehmittelpunkt reduziert wer- 
den muß, mz die Kurbelzapfenmasse und m,, die Nadellagermasse. Die oszillierende Massenkraft 
heißt damit 


Fo= mo& (4.2.45) 
und die rotierende Massenkraft 
Fa=mara®. (4.2.46) 


Die auf den Kolben einwirkende Kraft, z.B. durch Zündung eines Gasgemisches und der damit ver- 
bundenen Druckzunahme, sorgt für eine Entspannungsbewegung des Systems, d.h. eine Vergröße- 
rung des Zylinderraumes durch Kolbenbewegung. Die Kraft liegt in der Regel indirekt als Indikator- 
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Bild 4.32 
Kurbelwange 


Pi 


dwW=A, p(S)dS 


Bild 4.33 
Indikatordiagramm 


Bild 4.34 
M Drehmomentenverlauf 


diagramm (Bild 4.33) vor. Diese praktische Meßwertaufnahme zeigt den Zylinderdruck über dem 
Weg. Die obere Kurve stellt die Entspannungsphase und die untere die Kompressionsphase dar. Die 
schraffierte Fläche ist ein Maß für die geleistete Arbeit. Die Kolbenkraft ergibt sich über die Kolben- 
fläche und den indizierten Druck zu 
2 
kein. (4.2.47) 
4 

Die Kompressionsphase wird durch die in einem Schwungrad bei der Entspannungsphase gespei- 
cherte Energie übernommen. Das Schwungrad ist für die Laufruhe eines Motors von entscheidender 
Wichtigkeit. Durch die Triebwerksbewegung und durch die Veränderung des indizierten Drucks er- 
geben sich Drehmomentenverläufe, wie sie Bild 4.34 wiedergibt. Daraus resultiert ein mittleres 
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Drehmoment Mdn. Die Abweichungen von diesem kennzeichnen das Arbeitsvermögen W. Dieses 
wiederum bestimmt das Trägheitsmoment der Schwungscheibe. Aus der vorhandenen Winkelge- 
schwindigkeit und einem angenommenem Ungleichförmigkeitsgrad ergibt sich das Trägheitsmoment 
aus der Gleichung 


Id= —;. 4.2.48 
32 ( ) 
Der Ungleichförmigkeitsgrad ist das Verhältnis der Differenzen der größten und kleinsten Winkel- 
geschwindigkeit der Schwungmassen Ymax BZW. Ymin ZU ihrem Mittelwert. Er wird aus Erfahrung 
bestimmt. Der Durchmesser der Schwungscheibe nach Bild 4.35 ergibt sich aus der Ableitung 


„aM na_,s 
Id= g (D’-d)b (4.2.49) 
zu 
_+/32ld , 4 
D = m, +d. (4.2.50) 
75 b 


Bild 4.35 
Schwungscheibe 


Kommen wir zum konkreten Programm. Dazu wollen wir, mittels vorhandenen Drucks, die Bewe- 
gungsverhältnisse des Kurbeltriebs wiedergeben. Die Druckeingabe soll dabei durch ein Unterpro- 
gramm geschehen, um diesen nach einem vorliegenden Indikatordiagramm oder mittels funktiona- 
ler Verhältnisse anzugeben. Bild 4.36 zeigt das entsprechende Flußdiagramm. 


Tabelle 4.20 Speicherplatzbelegung 


00 At 05 r 10 F; 15 Wi 

01 Ay 06 I 11 8; 16 &Md Ay 
02 y 07 d. 12 Fo 17 Zw 

03 ma 08 x 13 90-9; 18 n 

04 mo 09 Ax 14 t; 19 0 


4.2.6 Realer Stoß fester Körper 


Das Zusammentreffen zweier bewegter Massen m, und ma, bezeichnet man als Stoß. Dieser Vor- 

gang unterteilt sich in zwei Phasen. Dies zeigt Bild 4.37. Die in der Realität teilelastischen und teil- 

plastischen Massen m, und m,, verlieren einen Teil ihrer kinetischen Energie durch Umformarbeit. 

Am Ende der ersten Phase bewegen sich beide Körper mit der Geschwindigkeit 
mi vı + m2v2 


m, tm, (4.2.51) 
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LBLA/Start 


Eingabe 


du. ur, Mo 
Ma. #0. Ay 


xerli-singra)) til... 
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(1) 


M=vm-ı+rAp 


Ausgabe 
%, Fi, @i, Veiv 
at, Mdj, pi 
Öt, 


LBLB /Start Auswertg. 


Eingabe 
ö,d,b 


ur 
— = p (Dichte) 
g 

ZMd Ay 
mn Eier Eee 


Ausgabe 
m. Mdm 
Id 

D 


R/S / Ende 


Bild 4.36 Flußdiagramm zur Schubkurbelbewegung 


Tabelle 4.21 Programm Schubkurbel 


Start/Eingabe: 
000 76 LBL 
001 11 A 
002 47 CMS 
003 07 7 
004 42 STD 
005 00 00 
006 91 R/S 
007_ 939 PRT 
008 72 ST# 
009 00 00 
010 97 1SZ 
01i 00 00 
012 00 00 
013 06 06 
014 98 ADY 
Ausgabe 9: 

015 43 RCL 
ie 02 02 
017 42 STO 
018 19 193 


Berechnung x: 


019 71 SBR 
020 16 A' 
Rücksprungmarke: 
021 76 LEL 
022 43 REL 
Berechnung x: 
023 71 SBR 
024 16 A' 
025 43 RÜL 
026 02 02 
027 _ 99 PRT 
028 43 RCL 
0239 08 085 
030 99 FRT 
Eingabe p;: 

031 71 SER 
032 15 E 
Fi: 

033 65 % 
034 43 RCL 
035 07 07 
036 33 K2 
037 695 * 
038 89 ı 
039 55 + 
040 04 4 
041 95 = 
042 99 PRT 
043 42 STO 


Ü 


a 
De u Be Ba 
[ui 


m 


a Bu 
-m 


Der] [X Y] 
aıDe— 
mio 


nn 
par) 


ERFER 
NE 


< 
z 

w 
z 


SEES ET ET EPG EPG GEEFGFEFGFEFGFEFGRENG 


Rinne mnmnmn 
NIT II Or ASTFN EL 


Cena LEN U E EN 


ram 


"DIT Jin ec 


Hz Mc 


N eoinegiainen 


"Je &opf 


DOoQDLm 


DUO ATO RAM 


LOL EIN 


5. 
Ur] 


L 
1% 


mm 
ZN 
am 


= 
Am 


= 
Nomen 


r 


nm 
a 
4 


Abschluß: 

14E 69 8 

147 43 FREL 
142 01 01 
149 935 = 

15044 SUM 
151 18 18 
152 43 FÜL 
153 14 14 
154 93 FET 
155 01 1 


156 44 SUM 
157 18 18 
158453 FRÜL 
159 01 01 
le0 44 ZUM 


lei us 02 
ie2 95 AV 
163 61 GTD 
164 43 RÜL 
Unterprogramm 
Berechnung x: 
169 TE LBL 
lee 18 A’ 
1er 52% 
185 42 RÜL 
169 05 05 
170 €&9 8 
il 53% 
172 0 A 
173 759 > 
174 42 REL 
175 0202 
irE 28 5IM 
Ir7 945 
Ira 89 + 
179 43 FEL 
120 06 06 
1E1 E08 
182 53% 
123 CL 1 
184 75 - 
122 22 & 
126 01 1 
ler 727 
185 520% 
159 42 REL 
120 05005 
131 35 = 
132 45 

133 

124 

1939 

136 

137 

120 

133 

zuc 


Auswertungsprogramm 


4 } vE LEL 13 19? 284 55 > 
24 T l2E 8 S4 2855 89 ı1 
Ss4 92 ADW 3 = 286 91 RS 
Eu 2 ADW 33 FET 287 99 PRT 
45 Exil =Ü AD 5 ADY 283 55 + 
02 08 42 FEÜL 43 FÜL 289 91 R/S 
De: a 1r Y 1& 1E 290 939 PRT 
43 RUL su # So I=1 291 85 + 
203 08 08 43 RÜL 35 > 292 91 RS 
si0 545 15 18 4: 293 99 PRT 
211 4 5 = co 294 45 Y% 
212 aa FRT 23 295 04 4 
213 42 STO 55 296 95 = 
co 00 “1 297 45 Y*% 
Unterprogramm 42 RCL aa 298 93 . 
Eingabe p;: 1E 18 ss = 299 02 2 
24, mE sa + “2 ADY 300 05 5 
ac, ER 99 PRT 301 95 = 
gen 43 RCL 33 ADY 302 98 ADY 
32 PRT us ne ae ee 
> ETH 2 CL = a 304 91 R/S 
1. Phase 2. Phase 
[a ie TEN, nn ee Ann N 
\vı ) Y2 ZN u N 5 & Bild 4.37 
MM, EL Stoßvorgang 
m 
m, 2 "m tm, “m, ma 


Danach wird der elastische Anteil der Umformarbeit wieder in kinetische Energie umgesetzt und 
die Geschwindigkeiten der Massen lauten zum Ende der zweiten Phase 


k{vi v2) ma 


mı tm, (4.2.52) 


cu 


und 
k(v} = va) mı 


mı tm; (4.2.53) 


S-u + 
Dabei sind die Vorzeichen der Geschwindigkeitsrichtungen zu beachten. Aus der Energiebilanz 
folgt der Energieverlust beim Stoß 
2 

(1-82). (4.2.54) 


mım2 (viva) 


AE = mı+m, 2 


Der in diesen Gleichungen enthaltene Stoßfaktor k berücksichtigt die plastischen und elastischen 
Eigenschaften der Massen. Er ergibt sich aus den Relativgeschwindigkeiten vor und nach dem Stoß 
zu 

27 Cı 


k = mg (4.2.55) 
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Mit Hilfe dieser Gleichung läßt er sich durch ein einfaches Experiment bestimmen. Nach Bild 4.38 
fällt eine Masse m, aus der Höhe h, auf eine Unterlage ma. Aus der Relation der Rücksprunghöhe 
h, zur Fallhöhe h; , über die Umrechnung der potentiellen in kinetische Energien vor und nach 
dem Stoß, bestimmt sich der Faktor aus 


h} m, 
k -y Ba (4.2.56) 97 
ha 


Bild 4.38 2 

, So } h, 
Experimentelle Bestimmung Ion! 
des Stoßfaktors BIEERSR. CRD: ZER DR 


Häufige Stoßzahlen zeigt Tabelle 4.22. Der Berechnungsalgorithmus folgt in Bild 4.39. 


Tabelle 4.22 
Stoßfaktoren 
Material k 
Elfenbein 8/9 - 
Stahl 5/8 Eingabe 
Kork 5/9 my, V1,M2, v2 
Glas 15/16 
Tabelle 4.23 
Speicherplatzbelegung 
00 Zähler/m, tm 03 ma 
01 k 04 v4 
02 v2 05 m, 
06 K(v} - v2) 


Ausgabe 


u, C4, Ca, AE 


Bild 4.39 


Flußdiagramm zum Stoßvorgang 
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Tabelle 4.24 Programm Stoßvorgang 


non Tem LEL cdze 99 = Ss4 078 43 RCL 
ol 1i A üÜ sa + 4: 5TO 079 03 03 
alz 47 CMS in SEN le 06 080 55 + 
o03 05 5 ir 43 RÜL BIN 081 43 RCL 
ca4 de STO Ü 0S 05 43 REL 082 00 00 
005 Qu Cl c Sa + 203 083 65 x 
codes 31 Ra if 43 ECL SE + 84 53 € 
ol? 33 PET Ü ne 3 43 FÜL 085 43 RCL 
aus TE 5ST# [ an 0 086 04 04 
03 00 Gü 5 = 087 75 - 
oil ar DISE 94 FRT 088 43 RCL 
oil 00 cd SE HIT 089 02 02 
ol2 00 00 0980 54 > 
112 06 06 091 33 2 
ol4 98 Aw 092 55 * 
(15 42 FÜL b = 0983 02 2 
Ois 05 05 d= Fü 094 65 x 
DlF EI KH ol 01 095 53 <& 
I1S 42 RÜL EI .H 0986 01 1 
013 04 04 SE ı£ 097 75 - 
v2ü 85 + 42 FÜL 098 43 RCL 
2 42 RÜL a4 0 099 01 01 
03 03 TE 100 33 x2 
EI .M 4 FÜL 101 95 = 
42 RÜL sm ca de 102 99 PRT 
[u Po u oil S4 5 103 91 RS 


Verläuft die senkrecht auf den Berührungsflächen stehende Stoßnormale durch die Schwerpunkte 
beider Massen, so spricht man vom zentralen Stoß, dessen Gesetzmäßigkeit wir vorher abgehandelt 
haben. Beim exentrischen StoR unterscheidet man nach der Geschwindigkeitsrichtung beider Massen 
den geraden und den schiefen Stoß. Der gerade exzentrische Stoß nach Bild 4.40 läßt sich auf einen 
zentralen Stoß zurückführen, wenn man die reduzierte Masse 


Idı 
Med 7 22 (4.2.57) 
Bild 4.40 
und die Geschwindigkeit Gerader exzentrischer Stoß Y2 
vi = awı (4.2.58) 


einsetzt. Beim schiefen Stoß sind nur die Geschwindigkeitskomponenten bezüglich der Stoßnormalen 
am Stoß beteiligt. So gelten als Beispiel für einen schiefen Stoß gegen eine Wand nach Bild 4.41 die 
Gesetzmäßigkeiten 


Var = Yır (4.2.59) 
und 


Von=—-Kvmn. (4.2.60) 


vr 2 


Bild 4.41 Schiefer Stoß gegen eine Wand 


Gesetze über den Drehstoß lauten analog und sind im gleichen Sinne programmierbar. 


112 


4.2.7 Anwendungsbeispiele 
Bon, 


Das Massenträgheitsrmoment des in Bild 4.42 dargestellten Körpers ist zu bestimmen (Dichte 


7.85- 10°? kg/cm?). 


Der Körper unterteilt sich in drei Grundkörper, in zwei Quader und einen Hohlzylinder. Ihre 


Massen betragen 


m; = 246.6kg 
ma = 117.8kg 
ma = 157 kg 


Damit folgt durch das Programm 4.12 


Eingabe: 246. & mı Zuvor Programmteil 


30. RB (Hohlzylinder) 
10. r aufgerufen, 


so.h 


Ausgabe: 12330. Id, 
10275. Id, 


Bild 4.42 


Rotationskörper 


Drehachse 


(Maße in cm) 


se gut 


Da Drehachse-Hohlzylinder und Drehachse-Körper zusammenfallen, braucht eine Umrechnung 
nach dem Steinerschen Satz nicht durchgeführt werden. Für die beiden Quader ist dies jedoch der 


Fall. Dazu wird Programmteil E’ aufgerufen. 


2. Teil 
Programmteil A 
Eingabe: 117.8 m 
SD. 2 
30. b 


Ausgabe: 33376. bE6nT Ida 


Programmteil E’ 


Eingabe: 33376. 66657 
117.8 


30. a 


Ausgabe: 139396. 6667 Id,, 


Damit ergibt sich 
Id = 129.44 kgm? 


3. Teil 
Programmteil A 


Programmteit E’ 


S3641, E6667 
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nn 
Das Massenträgheitsmoment des in Bild 4.43 dargestellten Körpers ist zu bestimmen (Dichte 
7.85: 10 ®kg/em?). 


(Maße in cm) 


Bezugslinie für die Höhen 


Bild 4.43 
Rotationskörper 


Manuelle Be- 26.5 Ayo 
rechnung der 32, hes 
Höhendiffe- 4 1 . 5 hgo 
renzen: 


[5 


me 
JE nn 
“ [1 vo 
iv 
E2 
3 
oO 


Bı 
Ei: 

o 

a 


AinNinianin“ 
[°] 
a 


nl 
2 

o 

a 


Anwendung von Programm 4.14: 
Durch den Höhensprung bei r = 80, ist eine Aufteilung der Berechnung notwendig. 
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Eingabe: D. DG7TSS 5 =, su. 
zu. n 7.5 10.5 
50. 12 ZI9E50, 1250000. 
5. Ar 
8. hi , 
Ausgabe: : 5 
Eingabe: 20. BEINE, © 5 
Ausgabe: 8 hia 
El, Tim Ar E End- 
en NE Ausgabe: 
25. 12000. BSUOES zur Id; 
2 
125000, 45 
5 
BEO1l2S, 
Eingabe: 0. ES 
So. rn Damit folgt 
0. nn 
5. Ar Id = Id, + Id, = 1730 kgm?, 
41.5 hi 
ie Fa 2 
=9, h 
206802000.  ri’hmAr 
Ausgabe: go 
Eingabe: 26 


EI-Pigr 
Ausgabe: 23153218. 75 
End- 


Ausgabe: 
10799572.42 Ida 


r?hmıAr (cmP) Bild 4.44 


Verlauf des Massenträgheits- 
moments über dem Radius 


r(cm) 
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a 
Welche Schwingbewegung vollführt ein nach Bild 4.45 aufgehängtes Pleuel mit Id, = 35 kg cm? bei 
30° Auslenkung. 


Das Massenträgheitsmoment bezogen auf den Drehpunkt 
Drehpunkt ergibt sich aus (4.2.12) 


Ido = Idg + ma? = 194 kgem?. 


Damit ergibt sich wiederum nach (4.2.20) (Maße in cm) 
Ido 
Ira = ir 17.2 cm. 25,5 


Bild 4.45 


Physikalisches 
Pendel 


Mit Hilfe der zuvor ermittelten reduzierten Länge, läßt sich nun Programm 4.5 anwenden. 


Eingabe: Em 3. 323910178 
2a 12. 44685021 
vo . 2140953253 
1. wo 3, 962432152 
to 0.2 
0. At 
02 Atpre -23, K020EBOT 15.62561513 
-3, 840180447 11. 099659633 
-, 0EEOS11URT „1909147772 
Ausabe " 7. 1FOS1ETEE 9. 417638999 
0.22 


12 
3. 01258562 
0.24 


30. 1690215 
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Yard) 


t(s) 


Bild 4.46 
Schwingungsverlauf 


-10 


-20 


-% 
Bild 4.46 zeigt den so ermittelten Schwingungsverlauf. 


za 
Bei einem Kolbenmotor mit den Daten 
Ik =0.1m 
I =0.3m 
r =0.05m 
mo = 80kg 
ma "7 60kg 


wurde ein Indikatordiagramm, Bild 4.47, erstellt. Gesucht ist eine Analyse vorhandener Bewegungs- 
und Belastungsgrößen bei einem Ausgangswinkel yo = — 90° und schrittweises Vorgehen um 10°. 


p; (bar = 10° N/m?) 


Bild 4.47 
Indikatordiagramm 


01 117 


Eingabe: 


ja Den} 


Toner 


efddan® 
oo 


ae] 


in“ 


agıin® 


. UOOER4FEE 
‚ocie3a4047 
7. 274598038 
a 4relzb4ld 


el, 


0. 
. 0541220108 
= 


zanch, 


=, [POS2E209 

. 372151599 
.o0444sl 
-, SoO3202143 
-, A304EI1477 
5. 354260507 
24, 30494442 


40. 


„azrl3j‘ 


1.£&0 33112 
.GN4STSEIEE 


-. Go0c1>5003 
241 


Ss0. 
‚0124243123 
20000, 

137. 0796327 
131. 2095128 
5. 613358017 
1. 512091131 
. 0043575651 
. 0002104713 
« 0168377008 
4. 400155953 
39. 16826522 


60. 

. D07r422112 

55000. 

431. 9689899 

332. 1576497 

11. 43325048 

. 5746419067 

. 1064964885 

-. 1022168197 
-. 1773455762 
9. 235050295 

40. 04290712 


70. 

0, 003503173 

100000, 

785. 3981634 

753. 3701526 

15. 84687301 

0. 631039322 

. (067175501 

-. 5003503135 
-. 0280250805 
11. 32415617 

40. 67394645 


30. 

0. 000585279 
130000. 
1432. 25651 
1477.08943 
22.13924988 
. 45065686821 
. 5058089102 
. 0020162038 
.1512963018 
10. 82839923 
41. 12461513 


90. 

0. 

500000. 
4712. 38898 
4712. 38898 
39. 63327298 
‚2523132322 
. (035086506 
« 0091166818 
‚7293345417 


D. 
41. 237692838 


100. 

0. 0003835274 
830000. 
4948. 005429 
4897. 714267 
40. 40509155 
. 2474935613 
-. 0035769781 
. 0236295475 
2. 290363803 
-35. 20635005 
41.62442194 


110. 

0. 003503173 
700000. 
3497. 787144 
5273. 591065 
41. 92638504 
. 2385104445 
-. 0109760137 
. 0310218517 
2. 481748136 
-79, 26909318 
41. 86293239 


120. 
‚DOrr7422112 
. 735000. 
Sr772, 676501 
5240. 652228 
41. 79574232 
. 2392588203 
-. 0177173751 
. 0281760205 
2. 254081638 
-123. 4138537 
42. 10219121 


130, 
.0134243123 
750000, 
5890. 456225 
4920. 356729 
40. 49835157 
. 2469234476 
-. 0230115899 
‚0214407129 
1. 715257031 
-165. 0058482 
42. 34911465 


140. 
.D203157661 
r20000, 
5654. 866776 
4192. 510535 
37. 33319825 
0. 267499852 
-. 0257624533 
0, 01028362€ 
„8226900777 
-193. 1985133 
42.61661451 


150. 
„0281414478 
ESllld, 


Be DE URS 


Dir Ba On 


pure 


una“ 


Kr 


„Lich 


nasseoreas 
0150924159 
-, (3484570 
a Be 
54. 1593995 
aa Pa11a9rs 


180, 
.Gaglsenlls 


Ds 
-, ODE 


52, 41734637 
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Nach Abschluß der Berechnung 


Auswertung durch Programmteil B vektorielle 


Darstellung 
Ausgabe: mist) 
Md,n {Nm} 

Eingabe: LE 5 Maßstab 
Ausgabe: F10, FERITE  Idikgm?) 1000 N 
Eingabe: r350, p 

ol b 

0,05 d 

Bild 4.48 

Ausgabe: D (m) Rotierende Kräfte 


Damit lassen sich die nachfolgenden Diagramme erstellen. (Ohne Einfluß eines Schwungrades.) 


Ylord] 


geringe Kräfte beim 
Anfahrvorgang 


Bild 4.49 
Oszillierende Kräfte 
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ol[s'] Bild 4.50 
Winkelgeschwindigkeit der Kurbel 


yl[grd] 


v[m/s] Bild 4.51 
Kolbengeschwindigkeit 


» [grd] 
270 
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Md [Nm] Bild 4.52 
Arbeitsvermögen 
200 


0 90 180 270 


- 200 


5 — 

Zwei Massen, m, = 40 kg und m, = 25 kg aus Stahl, treffen mit den Geschwindigkeiten v; = 5 m/s 
bzw. va = 2 m/s zentral aufeinander. Wie verhalten sich beide Massen nach dem Stoß und welche 
Umformarbeit wurde geleistet? 


Den Sachverhalt gibt Bild 4.53 graphisch wieder. 


Eingabe: my 
v1 
ma vi v2 
v2 
u 
Ausgabe: 
4 C2 


Bild 4.53 
Geschwindigkeiten sich stoßender Massen 
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4.3 Mechanische Schwingungen 


Unter einer mechanischen Schwingung versteht man die periodische Bewegung einer Masse um 
eine Mittellage. Den einfachsten Fall bildet ein Feder-Masse-Systemn. Bei der Bewegung findet ein 
ständiger Energieumtausch zwischen potentieller und kinetischer Energie statt. Die potentielle 
Energiedifferenz wird auch als Federenergie bezeichnet. Die bei der Bewegung umgesetzte Wärme- 
energie, durch innere Reibung in der Feder, soll unberücksichtigt bleiben. Wirken auf ein schwin- 
gendes System keine äußeren Kräfte, bezeichnet man den Bewegungsvorgang als freie Schwingung, 
Andernfalls als erzwungene Bewegung. Die bei der realen Schwingung stets auftretende Wider- 
standskraft, Bewegung im Medium und Reibungskraft (Stokes’sche Reibung, im Gegensatz zur 
Coulomb’schen oder Newton’schen Reibung) ‚etc., soll in erster Näherung als geschwindigkeitspro- 
portional angenommen werden. Dies entspricht auch unserer Annahme in 4.1.1. 

Nach der Schwingungsform unterscheidet man Längs-, Biegungs- und Torsionsschwingung. 


4.3.1 Freie Schwingung 


Wir betrachten den allgemeinsten Fall, eine freie gedämpfte Schwingung. Die zum Zeitpunkt t an 
der Masse angreifenden Kräfte zeigt Bild 4.54. Danach wirkt am Körper die Federkraft 


F=fs (4.3.1) 


mit der Federkonstanten f 
und die Dämpfungskraft 


Fu = 2mds, (4.3.2) 
Fu 
Bild 4.54 Fr 
Freie gedämpfte Schwingung “ 
m'a 


mit der Dämpfungskonstanten d, als Maß für die Dämpfungsintensität. Nach dem d’Alembertschen 
Prinzip folgt 


ms=-fs-2mds (4.3.3) 
dv f 
“m s-2dv (4.3.4) 


und damit nach dem Euler-Cauchy-Verfahren 
f 
Av=- ms+2dv At. (4.3.5) 
Wir erhalten den in Bild 4.55 dargestellten Berechnungsalgorithmus. Er ist in der Programmform als 
Rumpfprogramm aufgebaut, d.h. die entscheidende Berechnung der Geschwindigkeitsänderung ge- 


schied im Unterprogramm. Damit können auch nachfolgende Probleme in dieses Schema integriert 
werden. Setzt man d = 0, liegt der Fall einer freien ungedämpften Schwingung vor. 
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LBLA /Start 


Eingabe Zur Eingabe wird nach Aufruf des 
mtdsw uU-___-_-- Programms A der Wert 8 eingegeben, 
to, At, Ator da B Werte einzulesen sind. Siehe 

Tabelle 4.25 
f 
Au=-(Es +2dv..,) At 
Bild 4.55 
Flußdiagramm zur freien gedämpften 
Schwingung 
Tabelle 4.25 
Speicherplatzbelegung 
00 Atyorh 05 5/5; 
01 Aton 06 d 
02 At 07 f 
03 toft; 08 m 
04 vo; 09 c=0 für diesen Fall 


Die Annahme einer linearen Federkennlinie ist nicht immer ausreichend genau. Daher soll nachfol- 
gend noch ein Programmteil mit der Federkraft 


Fı= fs(1+cs?) (4.3.6) 


aufgestellt werden. Der linearen Federkennlinie ist also eine kubische Parabel überlagert. Die Ver- 
gleiche werden im nachfolgenden Anwendungsbeispiel durchgeführt. Wir erhalten in ähnlicher 
Weise 


Au=-(4s+ 88420) ar (4.3.7) 


und damit lediglich einen zusätzlichen Term. Unser Unterprogramm bekommt die in Tabelle 4.27 
wiedergegebene Form. 
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Tabelle 4.26 Programm freie gedämpfte Schwingung 


Start/Eingabe: Ass: Ausgabe: 
DOO 7&E LBEL 013 43 RÜL ü42 7&ö LEL ıF ur 
001 11 A 020 04 04 ü43 93% PRT Es + 
002 47 CMS 021 65 & a44 00 0 42 REL 
003 91 RS 022 43 REL 045 42 5TD oe 08 
004 42 STO c23 02 082 dds G0 00 65 
005 00 00 c24 95 = üdY 42 FRÜL 43 RÜL 
006 931 RS 025 44 ZUM a4 05 05 os 05 
007_ 99 PRT zer 05 05 049 93 FRT sa + 
008 72 ST# SO 42 RÜL DE 58 
009% 00 00 t;: csi 04 04 EB 
010 97 DSZ e PEAFTR üs2 29 FRT 43 RÜEL 
011 00 00 Da ne 053 43 RÜL 00% 
012 00 00 Ds bz 02 054 02 03 65x 
013 06 06 BER 44 SUM 055 9% PRT 42 RÜL 
014 98 ADY vhn us 08 056 35 AM 04 04 

31 44 SLM 057 &1 GTO E; 5 

32 0 nl 4 
Berechnung: - = Ss 055 42 RÜEL a4 +. 
015 76 LBL Abfrage/Ausgabe: SBR fü &5 2 

= ür za PL 

016 43 RCL (22 42 RCL Av,v;: er Ber 
Aufruf/Unterprogr.: 034 iu 1 ol re 
01? 71 SBR 34 sum 
018 38 SIN 37? 00 00 24 \ 2 u 

cas 77 GE 243 52 Wü En 

0393 93 PRT 244 43 RÜL 

ü4ü &1 GTO 

041 42 REL 
Tabelle 4.27 SBR für nichtlineare Federkennlinie 
240 78 LBL 250 42 RÜL BECHE G4 + - 
241 238 SIN 251 05 05 ss + EIKE 
242 53 £ E52 85 + u Por 43 REL 
243 53 & 253 24 CE EIOH je u FE 
244 43 RÜL 254 65 8 42 FRÜL s4 
245 07 07 255 42 RÜL le Dr 44 SuM 
246 55 + 2SE 00 008 EI =M a4 og 
247? 43 RLL 287 89 42 FÜL “2 RETH 
248 08 08 E58 42 REL 4 04 
249 69 x 2593 05 05 4 


Wird die Dämpfungs- oder Reibungskraft, als Newton’sche Reibung mit 
Fa = esgn (v) v2, (4.3.8) 


also dem Quadrat der Geschwindigkeit proportional, angesetzt, ergibt sich die Differenzengleichung 
zu 


Av=- (# s + csgn (v) 2) At (4.3.9) 


und damit ein Unterprogramm nach Tabelle 4.28. 
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Tabelle 4.28 SBR für Newtonsche Reibung 


54 5 
94 +7- 
639 X 

42 RÜL 
2 02 
Ss4+ 5 

44 SuM 
ü4 04 
92 RTN 


4.3.2 Erzwungene Schwingung durch eine rotierende Masse 


Zum Ansatz betrachten wir das idealisierte Schwingungssystem nach Bild 4.56. Die im Abstand r 
außerhalb des Drehpunktes rotierende Masse m, , hat in Schwingungsr ichtung den Fliehkraftanteil 


Fe = mırw* sin (wt). (4.3.10) 
Der Ansatz gestaltet sich wie zuvor 
m$=-fs-2mdv-mı rw? sin (wt), (4.3.11) 


wobei m die Masse m, beinhaltet. Es folgt wieder 
Av=-(t +2dv-rolsi A 4.3.1 
v=-I75 vomtw sin (at) t (4.3.12) 


als Bewegungsgleichung. Der nachfolgende Berechnungsalgorithmus in Bild 4.57, ist für konstante 
Winkelgeschwindigkeiten ausgelegt. Ansonsten muß der Algorithmus eine zusätzliche Gleichung 
w = f(t) bestimmen. 


FERIEN A Tabelle 4.29 


Zusätzliche Speicherplatzbelegung 


09 w 
10 r 
11 mı 


Bild 4.56 


Erzwungene Schwingung 
durch Rotor 


Auch hier läßt sich wiederum der Programmrumpf von Tabelle 4.26 verwenden. Mit einer zusätz- 
lichen Speicherbelegung zu Tabelle 4.25 folgt Tabelle 4.29. Bei Aufruf des Programms A, muß an- 
schließend 11 eingegeben werden, da 11 Werte einzulesen sind. Tabelle 4.30 zeigt das entsprechen- 
de Unterprogramm. 
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Tabelle 4.30 SBR Erzwungene Schwingung durch Rotor 


240 76 LBL 59 FE nB 42 RC 
241 38 5SIH 22 65 03 = 
242 53 %& x 42 RÜ 5 \ 
243 53 % = FÜL 10 1u == 51h 
244 43 RCL O6 DE ESH 3 . 
245 07 07 EIN 43 FCL be] Er 
246 55 > 43 REL 230% BI 
24? 43 RÜL od 04 43 RÜL 
eg 05 05 TS - [1 us 
249 65 * 43 RÜL 
250 43 RÜL ii 1l 
251 05 03 sa + 

42 RÜL 


Eingabe Zur Eingabe wird nach Aufruf des 
mı,r,w,m,f,d Programms A der Wert 11 eingegeben, 
So, Vo, to, At da 11 Werte einzulesen sind. Siehe 

Ator Tabelle 4.29 


Bild 4.57 


Fiußdiagramm zur erzwungenen 
Schwingung 


4.3.3 Erzwungene Schwingung durch rotierende und oszillierende Massen 


Zum Ansatz betrachten wir wieder ein idealisiertes Schwingungssystem nach Bild 4.58. Der vor- 
handene Schubkurbeltrieb beinhaltet rotierende und oszillierende Massen. Sie werden als Massen- 
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kräfte 1. und 2. Ordnung bezeichnet. Ihr Anteil 
in Schwingungsebene beträgt für Massenkraft 
1. Ordnung angenähert 


Fı= mı rw? coswt (4.3.13) 


und für Massenkraft 2. Ordnung 


2. 


Fı= mırw cos? wt. (4.3.14) 


|. 


Die Differentialgleichung der Bewegung lautet 


ms=- (1s+ 2mdv- m; rw? (eos@t -7 coat). 
Auch hier beinhaltet m wieder m;. (4.3.15) 
Bild 4.58 


Erzwungene Schwingung durch rotierende 
und oszillierende Masse 


Aus dieser leitet sich wieder die nachfolgende Geschwindigkeitsänderung ab 


m 
v=- (# s+2dv- rw? (eosur-! cos®wr))At. (4.3.16) 


Beachten Sie bitte, daß im vorliegenden Fall der Winkel wt anders als bei reiner Rotation gemessen 
wird. Mit zusätzlicher Speicherplatzbelegung nach Tabelle 4.31 folgt das Unterprogramm für diesen 
Fall nach Tabelle 4.32. 

Tabelle 4.31 Zusätzliche Speicherplatzbelegung 


09 w 11 mı 
10 r 12 | 


Tabelle 4.32 SBR erzwungene Schwingung durch rotierende und oszillierende Masse 


2932 65% 
293 43 RÜL 
2934 09 09 
235 695 


z29E 43 RÜL 
297 03 03 


1] e933 54 2 

h 299 39 C05 
300 54 5 
201 54 5 
Sulz 44 +7- 
203 65x 
304 43 RÜL 
305 02 02 
306 54 >» 
307 44 SUM 


4.3.4 Schwingung unter Berücksichtigung der Federmasse 


Die Voraussetzung für die Betrachtung dieses Ansatzes besteht in der Annahme, daß jeder Teil der 
Feder einen gleichen Anteil an der gesamten Durchbiegung hat. Unter Betrachtung von Bild 4.59 
ergeben sich für die Durchbiegung an der Stelle s folgende Verhältnisse 


y ix (4.3.17) 


und auch 
u I 
yaıı (4.3.18) 
x“ 
Bild 4.59 
Schwingende Federmasse 
Der Arbeitsanteil eines Massenelements dm; beträgt zu diesem Zeitpunkt 
dm, .2 
dW = ru v”. (4.3.19) 
Daraus folgt durch Integration 
W= 3 [am y”. (4.3.20) 
Mit 
m; 
dm; = Ar ds (4.3.21) 
worin m, die Federmasse ist, gilt 
2 I 
=2 4.3.22 
We ): ds ( ) 
0 
wa.) (4.3.23) 
Er 3/72: ‚3. 


Dies heißt aber, daß die Federmasse nominell zu einem Drittel am Schwingungsprozeß beteiligt ist. 
Gleichung (4.3.5) ergibt sich damit ergänzt 


Av=- ( Ä mıst 2av) At. (4.3.24) 
m+7" 


Programmtechnisch ändert sich nichts, da statt der Masse m die ergänzte Masse 
n mi n n 
m'=m+ 3 eingegeben wird. (4.3.25) 
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4.3.5 Biegeschwingungen 


Wir betrachten ein Blattfederpendel nach Bild 4.60. Durch die Masse m wird eine statische Durch- 
biegung x erreicht. Nach Auslenkung aus seiner Lage, schwingt das Systern nach der Bewegungs- 


gleichung 


ms+fs=0. (4.3.26) 


Bild 4.60 Blattfederpendel 


Dies ist die allgemeine Differentialgleichung für freie Schwingung ohne Dämpfung. Die Federkon- 
stante f ergibt sich aus der Überlegung, daß die Kraft, im Bereich des Hookeschen Gesetzes, am 
Stab proportional zur Durchbiegung ist. Mit Hilfe der statischen Auslenkung ergibt sich also 


mg 


en 
Im Fall des Blattfederpendels 


3El 
re 


(4.3.27) 


(4.3.28) 


mit E = E-Modul des Trägerwerkstoffs und I = axiales Flächenträgheitsmoment des Querschnitts. 


Damit läßt sich dieses System als ein vorangegangenes allgemeines Schwingungssystem betrachten. 


So ergibt sich für einen Träger nach Bild 4.61 


48El 
f= a (4.3.29) 


Bild 4.61 
Schwingender Träger auf zwei Stützen 


5) 
37 


Berücksichtigt man die Federarbeit des Trägers, so ergibt sie sich aus der Entwicklung von x in y 


eingesetzt 
s 45? 
2 \ ah 
“3 (212 _ae21% 
ven (31° — 45°) x 
dm ax? 
dw = — 2? = — (9145? - 2412 5% + 1656) ds 
2 2gl 
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(4.3.30) 


(4.3.31) 


(4.3.32) 


(4.3.33) 


3 3 Be 
„nl Ss, 50/ Pr 
48El 384EI 


Bild 4.62 Be 
Schwingende Trägermasse yangls ( = 45) 5ai? ( _48 ) ( -35) 
16E1/ 31?/ 384EI 2 51? 


Es müssen also 2.stel der Trägermasse für deren Federarbeit berücksichtigt werden. Weiterhin wollen 


wir die bisher vernachlässigte Trägermasse beachten. Nach Bild 4.62 folgt auf gleichem Wege 


lee 
y- (1-2&)(-45)3 (4.3.35) 
dW = = y? (4.3.36) 
v2 5 
W= 2 [m = (2 m) en (4.3.37) 
o 


In diesem Fall müssen also 2 .stel der Trägermasse eingesetzt werden. Ein Geringeres als zuvor. 


Ein besonders wichtiger Hinweis zur kritischen Winkelgeschwindigkeit einer rotierenden elastischen 
Welle. Sie bestimmt sich, als statisches Träger-Masse-System aufgefaßt, aus der statischen Durchbie- 
gung durch 


&x=-Vx (4.3.38) 
Ein nachfolgendes Anwendungsbeispiel zeigt dies anschaulich. 


4.3.6 Drehschwingungen 


Ein Torsionspendel nach Bild 4.63 erfährt bei 
Auslenkung um y das rückstellende Moment 


GI 
Mg (4.3.30) 


Bild 4.63 
Torsionspendel 
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G ist der Gleitmodul des Fadens, Ip sein polares Flächenträgheitsmoment. Daraus folgt als Bewe- 
gungsgleichung für freie Drehschwingungen 


„__6lp 
Id-$=- T» (4.3.40) 
u GIp 
9-7 Tja (4.3.41) 
GIp 
dw = Id ydt (4.3.42) 


Aw =- — pt. (4.3.43) 


LBLA /Start 


Eingabe: 
G, Ip, I,Id 
Yo, Wo, to, At 


Das Flußdiagramm folgt analog zu Bild 4.64. 

Mit der Speicherplatzbelegung nach Tabelle 4.33 
läßt sich das Rechnerprogramm in 4.3.1 als Rumpf- 
programm auch hier einsetzen und wir erhalten 

ein Unterprogramm nach Tabelle 4.34. 


Tabelle 4.33 
Speicherplatzbelegung 


00 Atyomn 06 Id 
01 Aton 07 ı 
02 At 08 Ip 
03 toft; 09 G 
04 wolw; 10 
Tabelle 4.34 


Bild 4.64 
Flußdiagramm zur Drehschwingung 
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4.3.7 Anwendungsbeispiele 


Alle Programme in diesem Buch sind für die Anwendung mit dem Drucker geschrieben. Wollen Sie 
diese ohne ihn benutzen, dann müssen Sie in den Programmen die Befehle Prt durch R/S ersetzen. 
Das Programm hält an und Sie können das Ergebnis notieren. Danach starten Sie wieder mit R/S. 


Re 
Eine freie gedämpfte Schwingung hat folgende Daten: 
m = 50kg 
f = 80kg/s? 
d= 2, 0.4 und 0.045! 
sp =-5m 
Vo = 0 
to =0 
At= 01 Berechnungsdaten 
Aton = 1 
Eingabe: 50. m 1.451 
(d= 2) 80. ff SET set 
2. d re 
-5, so 
0. vo 
0. to 
0,1 At 
1. Ator 
Ausgabe: 


-3. 484475457 5; 
KESBOSEIS: x 
. j 


-2. 250307548 
1. 008055577 


2. 
Eingabe: 50. m 
(d= 0.4) 80 f 
0.4 d 
-5. so 
Ö vo 
0. tb = 
0.1 At 
1.  Atprt 
Ausgabe: SER DRAS 
VNZOSE 
-1. 989818698 a er 
4. 235892561 vi BaRuie 
1. 6 


1. 287383235 
1. 902153264 
2. 


Achtung! 

Nach Aufruf von A muß 
zuerst 8, für die Anzahl 

der Eingabewerte einge- 
geben werden. 


Zuvor Eingabe von 


B 


133 


Zuvor Eingabe von 
8 


Eingabe: 


Bild 4.65 
Schwingungsvorgang 


Aus der allgemeinen Lösung lassen sich Gesetzmäßigkeiten für diesen Schwingungsvorgang gewinnen. 
So liegt für den Fall 
Zi 


(4.3.44) 


m 


d> Y 


ein aperiodisches Schwingungsverhalten vor. In unserem Fall 


gilt dies für d= 25°". 
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Weiterhin liegen zwei Nulldurchgänge in derselben Richtung zeitlich um 


2 
Ta ——E- (4.3.45) 


auseinander. In unserem Fall mit d = 0.4 


T= ————— = 4,852225 


Weiterhin beträgt das Verhältnis der beiden Höchstausschläge zu einer Seite 


A; 


Td 
=e”, 4.3.46) 
Airı 


In unserem Fall mit d = 0.4 
eTd = 6.96493. 


Umgekehrt bezeichnet man die Größe 


A; 
In = Td (4.3.47) 
A; +1 


als Jogarithmisches Dekrement der Schwingung. Mißt man praktisch eine Reihe von Höchstausschlä- 
gen und die Schwingdauer, so läßt sich daraus die Dämpfungskonstante d bestimmen. Mit Hilfe der 
Gleichung (4.3.45) für T folgt bei gegebener Masse auch die Federkonstante. 


En 
Das unter — 1 — angegebene Schwingungssystem (d = 0.4 5°") wird auf unterschiedlichste Art er- 


regt. Der Schwingungsverlauf ist aufzuzeichnen. 
Die Art der Schwingungserregung wird durch die Größen so und vo wiedergegeben. 


1. Fall: Die Masse wird ausgelenkt und losgelassen (wie unter — 1 —) 
so=-5m; w=0 

2. Fall: Die Masse bekommt in der Ruhelage einen Stoß nach unten 
so=0; vo=-5mis 

3. Fall: Die Masse wird ausgelenkt und bekommt einen Stoß nach oben 
so=-5m; vo = 5m/s 

Diese drei Fälle mögen genügen. Den Schwingungsverlauf zeigt nachher Bild 4.66. 


Deutlich sieht man in Bild 4.66 die Einflüsse auf den Schwingungsverlauf des Systerns. 
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s{m) 


Bild 4.66 
Schwingungsvorgang 


-5 

Achtung! 

Zuvor 8 (Eingabewerte) 

angeben 

Eingabe: Sil, m -1. DAITSELET . 1443760255 

20. f 2, 424458 . 1434035718 

ü.4 od 2 : 
[so 

Pi . “ 
Es He . VILZIOITZE 

Ber & -. 1335993166 
L. Atprı 9 


Ausgabe: 


. 019936.85031 
-,. FA9ET1637 


10. 


«1343760958 
-. 2352749205 
8 


Eingabe: sn, 


-. 0616696129 
-. 1269023765 
9 


-. 0360645342 
« 0452593323 
10. 


en 
In diesem Beispiel vergleichen wir das Schwingungsverhalten des Systems aus — 1 —, d=0.45° I 
bei unterschiedlicher Federrückholung, d.h. bei linearer und kubisch überlagerter Kennlinie. (c = 0.1). 


Unter Benutzung des Unterprogramms für die kubische Federkennlinie ergab sich 


Achtung! 
Zuvor 9 (Eingabewerte) 
eingeben 


Eingabe: O.1e 


Ausgabe: 

1. 167327916 si 

5. 4802946322 u, 
1. t; 


kubische 
Federkennlinie 


lineare 
Federkennlinie 


ts) 


10 


Bild 4.67 
Schwingungsvorgang 
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-4- 
In diesem Beispiel wollen wir abschließend das System aus — 1 —, d= 0.45", bei unterschiedlicher 
Dämpfung betrachten. Dazu benutzen wir das Unterprogramm Newtonsche Reibung mit c = 0.1 


Achtung! 
Zuvor 9 (Eingabewerte) 
eingeben 
Eingabe: DL.le 1. 25075166 
SÜo. m -1. 20292091 
zo. f 5 
dd 
u) -. 4337673469 
J. vo -1. 714245239 
to ER 
O1 At 
1. A 
Rz -1. 329294807 
Ausgabe: + 0433376943 
-1. 553231172 5, 10. 


s(m) 


Bild 4.68 
Schwingungsvorgang 
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-5— 


Das Schwingungssystem von — 1 — wird durch eine außermittig rotierende Masse erregt. Die Daten 
sind 


m; = 20kg 
r = 0.05 m 
w= 605°! 
m = 70kg 
f = 500 N/m (also etwas steifer) 
d=04s' 
so =-5m 
Vo ” 0 
to =0 


Der Bewegungsablauf ist gesucht. 


Achtung! 
Zuvor 11 (Eingabewerte) 
eingeben 
Eingabe: 20. m, 2 
D.0S r 6 
El. w Er . 
?Ü. m 
soo. f RE EIER 
0.4 d ae 105 ne 
-5. 5 N un 
0. vo 
c. to 
0.1 At 
0.5 Atprt 


Ausgabe: 
-.0ois4182142 Si 
15. 401259824 Vi 
0,56 


8. 03678397 
13. 89884857 
1. 


10. 88307258 
1859381216 


S4a10llE 
1.5 


So0ssa4ad4 
4.9 


Aus Platzgründen müssen wir hier theoretisch einige Besonderheiten dieses Schwingungsvorganges 
betrachten. Man erkennt, daß die gedämpfte Eigenschwingung des Systems langsam abklingt und 
die erzwungene Schwingung immer mehr angenommen wird. Um dies deutlicher zu demonstrieren, 
wurde die Masse zum Start ausgelenkt auf — 5 m. Nach einer gewissen Zeit läuft die Bewegung nach 
dem Gesetz 


s=Asin(wt-y) (4.3.48) 


ab und damit mit derselben Frequenz wie die Erregerschwingung, nur um den Winkel % verschoben. 
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s{m) 


10 
erregte Schwingung 


=“ 


Erregerschwingung 
(stark vergrößert) 


freie 
Schwingung 


t(s) 
BE.) 


Bild 4.69 
Schwingungsvorgang 


Daher bezeichnet man % als Phasenwinkel. Er ergibt sich aus 


2 md 
tan y =. (4.3.49) 
f-mw 
Der Höchstausschlag A bestimmt sich durch 
mırw? 
A (4.3.50) 


Yt-mo2)2 +4m2 dw? 
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Bei Annäherung der Erregerfrequenz « an die Eigenfrequenz des Schwingungssysterns 


f 
x = jaze 


m (4.3.51) 


strebt A gegen ®. Dies veranschaulicht Bild 4.70. Daher bezeichnet man w, auch als kritische 
Winkelgeschwindigkeit. Der Vorgang selbst heißt Resonanz. Maschinenteile, die als elastische Bau- 
teile gelten können, müssen auf diesen Resonanzfall hin untersucht werden. Er kann zur Zerstörung 
des Bauteils führen. Zur Vermeidung gestaltet man > w,. Umlaufende Wellen sind solche Teile. 
Hier stimmen kritische Winkelgeschwindigkeit und Eigenschwingungszahl überein. Mit der statischen 
Durchbiegung 


mg A 
Vo. (4.3.52) | 
folgt | 
f g 
Rn -VE-y8 (4.3.53) | 
Yo | 


“x w 
Bild 4.70 Schwingungsausschläge 
-6- 
Dieses Beispiel soll einen Vergleich zwischen einer Schwingung, angeregt durch einen rotierenden 
Erreger und einer Schwingung, angeregt durch einen rotierend-oszillierenden Erreger darstellen. Da- 


zu benutzen wir das Schwingungssystem aus — 5 —. Zusätzlich wird die Schubstangenlänge | = 0.3m 
angenommen. 


Achtung! 
Zuvor 12 (Eingabewerte) 
eingeben 
Einbe: 0.3 | 8. 148437918 6. 98789366 9154477334 
20. mı 4. 018418821 3. 8432576 2, 257060809 
0.05 r 2. 5. 8 
60. ww 
sin, ? een 
0.4 a an 
-I.9 Er 
0. vo 
0. to 2. 311272375 
0.1 At 9. galt 
0.5 Atprı Ss 


Ausgabe: 


-, 3142700809 
14. 03739354 
0.5 


6. 673286165 
11. 7667563 
1. 


9. 395443469 
1. 232700043 
1.5 


E25 

36 

3.5 

"2.67 1D3T228 
-5. 741561797 
4, 


-10, 15920698 
1. 2244360453 
4,5 
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Bild 4.71 
Schwingungsvorgang 


Um beide Schwingungen überhaupt vergleichen zu können, mußte in 
Programm 4.32 die Funktion cos in sin abgeändert werdne 
Um beide Schwingungen überhaupt vergleichen zu können, mußte in 
Programm 4.32 die Funktion cos in sin abgeändert werden 


Bi, DER 


Auf einem Träger auf zwei Stützen, Stützweite 3m, wird mit einer neuen starren Stützvorrichtung, 
nach Bild 4.72, die Last von 1500 kg abgesetzt. Als Hauptträger werden 2 Doppel-T-Profile, für 
Opzul = 8000 N/cm?, verwendet. Gesucht ist die Eigenfrequenz des Systems und der Schwingungs- 
verlauf. 

Zunächst ergibt sich aus der zulässigen Biegespannung das notwendige Widerstandsmoment 


Mb 69 cm? 


W= = 
Opzul 


Unter Berücksichtigung, daß hier zwei Träger nebeneinander liegen, ergibt sich aus einem der übli- 
chen technischen Tabellenbücher der Träger 


1120, W,„ = 54.7cm?, I, = 328cm*, m = 33.3kg 
Damit ergibt sich für den in Bild 4.72 dargestellten Belastungsfall 


= nl 2). 
Ym BEI] 1 3R = 0.413cm 


Nach Gleichung (4.3.27) folgt 


Damit ergibt sich die Eigenfrequenz nach (4.3.51) zu 


f 
“x -ys =48.7°°' 


Den Bewegungsablauf zeigt Bild 4.73. 


Bild 4.72 
Schwingender Träger 


y 
Achtung! “= 


Zuvor 8 (Eingabewerte) 
eingeben 


Eingabe: 1516.78 m 
39E0000, Ff 


-.oo1lca=1T 
-. 13947659 


‚1 
je) 


11 

0. SE 

-0, 00412 so 1 
0. 

0. oo . DO413FLTTE 

0.005 At . 04S4E47 741 

0.02 Ator D. CE 


Ausgabe: 
-. 0019141438 5, 5 
1664209916 u u 
0.02 ı 


0,005 


s{m) 


PER IE EN 


Beet 


Bild 4.73 
Schwingungsverlauf 


Gesucht ist das Schwingungsverhalten des in Bild 4.74 dargestellten Systems. 
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fr = 2ro 
r2"fo 

Id, =6kgm 

Id; =3 kgm 
d=0,04 m 
I=0,6m 
G =8- 10 N/m? 
Ip = 2: 10. m® 


Bild 4.74 
Schwingungssystem 


Zur Betrachtung werden die Massen der Drehwellen auf ein Ersatzsystern nach Bild 4.75 reduziert. 
Nach Gleichung (4.2.23) folgt 


To 3 1 Ido 
Idıio = Id, rn = ad 


yı) 2 
Id2o = Id, (e) = Id, 


Ido = Idiot Idzo =45 kg m d an aD 
Ersatzsystem 
Achtung! 
Zuvor 9 (Eingabewerte) 
eingeben 
Eingabe: 8. 10 6 4.07% 
2. -07 Ip 52. FrF4g4S 
0.6 ] 
4,5 Jd 
4. vo .S24sg33021 .coos4sz 
D. wo 313. 1215119 207, 9138 
0. to 0. 06 
0.005 At , 
0.02 Atprı 


Ausgabe: 


-. 69816969323 vi 
-313. 7130933 wi 
0.02 tj 


vlard) 


Bild 4.76 
Schwingungsvorgang 
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Eine Grenzwertbetrachtung liefert REED PER RLELEG 


Fe 
I Id 
T= 27) Sp 
4.3.5 
T = 0.0812s A) 
Bild 4.77 


Torsionspendel zur Feststellung eines Id 


Abschließend sei noch erwähnt, daß in der Praxis ein Torsionspendel mitunter zur Feststellung des 
Massenträgheitsmomentes eines beliebigen Körpers benutzt wird. Den Versuchsaufbau zeigt Bild 
4.77. Der Vorteil liegt darin, daß lediglich die Massen m}, der Radius r, und die Schwingzeiten des 
Systems T und T, bekannt sein müssen. Die Zeit T ohne aufgesetzte Massen, beträgt nach Grenz- 


wertbetrachtungen 


PER 
I Id 
T=2 5 Sp‘ (4.3.54) 


Mit aufgesetzten Massen beträgt sie, unter Berücksichtigung des Steinerschen Satzes 


Ald+2m, 7) 


T,=2 GIp (4.3.55) 
Daraus folgt durch Gleichsetzung und Umstellung 
2 mir, 
Id = ze (4.3.56) 
1 
Ei 
(7) 


Wenn Sie eines der üblichen Stopuhrprogramme für Taschenrechner benutzen und dieses mit Glei- 
chung (4.3.57) verbinden, können Sie Ihren Taschenrechner zur Bestimmung eines Id benutzen. 
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VIEWEG 


Taschenrechner-Literatur 


Peter Kahlig 
Mathematische Routinen der Physik, Chemie und Technik für 
AOS-Rechner 


Teit I. Mit 71 Abb., 129 Beispielen und 34 Tabellen. 1979. VI, 178 S. DIN C 5 (Anwendung 
programmierbarer Taschenrechner, Bd. 3/1). Kart. 

Inhalt: Gamma- und Beta-Funktion, Kombinationen (Binomialkoeffizienten), Variationen 
(permutations, factorial powers) und ihre Logarithmen — Digamma-Funktion und ihre ersten 
sechs Ableitungen (Polygamma-Funktionen), Beta-Funktion und ihre ersten sechs Ableitungen 
— Exponentialintegrale, Integrallogarithmus, Integralsinus und -cosinus, hyperbolischer Integral- 
sinus und -cosinus. 


Dieser Band enthätl 13 ausgefeilte AOS-Programme für 30 oft benötigte Funktionen aus den 
Bereichen der Physik, Chemie und Technik. Zur Auflockerung und zur zusätzlichen Informa- 
tion dienen 71 Abbildungen, die fast alle vom Rechner selbst gezeichnet wurden, ferner 129 
Beispiele und 34 Tabellen. Der Tuning Kit im Anhang enthält 8 komfortable, universell ver- 
wendbare Sonderprogramme zum Zeichnen und Drucken. 


Teil 11. Mit 137 Beispielen, 71 Abb., 16 Tab. und einem Anhang: Logarithmisches Plotten und 
Erzeugung von Fehlerkurven. 1980. VIII, 1808. DINC 5. Kart. 

Inhalt: zeta-, xi- und Xi-Funktion von Riemann — eta-, kappa- und rho-Funktion, L-Funktion 
von Dirichlet — Polylogarithmen, chi-Funktionen von Legendre — Arcustangens-Integrale — 
Clausen-Integrale und Glaisher-Funktionen — Anhang: Arithmetische Funktionen — Logarith- 
misches Plotten von Kurven — Plotten der Ordinatenachse mit logarithmischer Teilung und mit 
inverser logarithmischer Teilung — Erzeugung von Referenzwerten und Fehlerkurven — Fehler- 
kurven zu Funktionsroutinen dieses Bandes. 

14 ausgefeilte AOS-Programme für 21 oft benötigte spezielle Funktionen. Zur Veranschau- 
lichung und zusätzlichen Information dienen 71 Abbildungen, die fast alle vom Rechner selbst 
gezeichnet wurden, ferner 137 Beispiele und 16 Tabellen. 


Hans-Joachim Ludwig 


Programmoptimierung für Taschenrechner (AOS) 


2., durchges. Aufl. 1980. X, 102 S. 12 X 19,5 cm (Programmieren von Taschenrechnern, Bd. 5). 
Kart. 


Inhalt: Wozu dient Programmoptimierung? — Techniken der Programmoptimierung — Ratio- 
nalisierung der Programmherstellung — Steigerung der Effektivität — Erhöhung der Betriebs- 
sicherheit — Verbesserung des Bedienungskomforts — Förderung der Flexibilität — Verringe- 
rung des Programmspeicherbedarfs — Verringerung des Datenregisterbedarfs — Verkürzung der 
Rechenzeit — Programmbeispiele. 

Ziel des Buches ist es, dem Benutzer von Taschenrechnern das geschickte Ausnutzen des Gerä- 
tes bis an die Grenzen seiner Möglichkeiten zu zeigen und dadurch die Leistungsfähigkeit eines 
Programms wesentlich zu erhöhen. Diese Anleitung zur optimalen Programmierung setzt keine 
speziellen Kenntnisse voraus. 


Lehr- und Lernsystem 
Mechanik und Festigkeitslehre 
von Alfred Böge 


Mit dem Lehr- und Lernsystem Mechanik und Festigkeitslehre liegt ein nach 
modernen didaktischen und methodischen Erkenntnissen gestaltetes Lehr- 
werk vor, das allen Ansprüchen, die heute von Dozenten und Studierenden 
gestellt werden, in jeder Beziehung gerecht wird. 


Alfred Böge 

Mechanik und Festigkeitslehre 

Unter Mitarbeit von Walter Schlemmer und Wolfgang Weißbach. Mit 605 Abb., 
26 Arbeitsplänen, 20 Lehrbeispielen und 16 Tafeln. 17., überarb. Auflage 1979. 
XIV,381 S. DIN C 5 (Viewegs Fachbücher der Technik). Gbd. 


Alfred Böge und Walter Schlemmer 
Aufgabensammlung zur Mechanik und Festigkeitslehre 


Unter Mitarbeit von Wolfgang Weißbach. Mit 516 Abb. und 907 Aufgaben. 
7., überarb. Auflage 1979. XII, 211 S. DIN C 5 (Viewegs Fachbücher der 
Technik). Kart. 


Alfred Böge 
Formeln und Tabellen zur Mechanik und Festigkeitslehre 


Unter Mitarbeit von Walter Schlemmer. 9., überarb. und erg. Auflage 1980. 
V,49S.DINC5 (Viewegs Fachbücher der Technik). Kart. 


Alfred Böge und Walter Schlemmer 
Lösungen zur Aufgabensammilung Mechanik und Festigkeitslehre 


Unter Mitarbeit von Wolfgang Weißbach. Mit 738 Abb. 2., überarb. Auflage 
1979. 176 S. DIN C 5 (Viewegs Fachbücher der Technik). Kart. 


Anwendung programmierbarer Taschenrechner 


Diese Reihe bietet den Benutzern programmierbarer Taschenrechner eine reich- 
haltige Palette von Aufgabenstellungen aus den Anwendungsgebieten der Natur- 
und Wirtschaftswissenschaften an, für die Programme zur numerischen Lösung 
entwickelt werden. 

Jeder Band behandelt ein in sich abgeschlossenes Themengebiet: Nach einer kurzen 
Einführung in die Theorie der jeweiligen Problemstellung wird der Lösungsalgo- 
rithmus entwickelt, das Programm dargestellt und kommentiert. 


Neben der direkten Nutzung der hier veröffentlichten Programme unterstützt diese 
Reihe den Leser wirkungsvoll bei der Ausarbeitung eigener Programmvarianten. 


Band 4: Statik — Kinematik — Kinetik 
für AOS-Rechner 


von Harald Nahrstedt 


Band 4 der Reihe Anwendung programmierbarer Taschenrechner zeigt an 30 ausge- 
testeten Programmen für den TI-59 den Einsatz des programmierbaren Taschen- 
rechners für wichtige Gebiete der Technischen Mechanik. Der Verfasser gibt zu 
jedem einzelnen Problem eine kurze Einführung in die Theorie und erklärt dann die 
Aufbereitung des Programms. Es folgen Beispiele für die Entwicklung eigener 
Programme. 


Ing. (grad.) Harald Nahrstedt ist Dozent an der VHS Hamm 
Anwendung programmierbarer Taschenrechner, 
Mikrocomputer. 


ISBN 3-528-04169-2 


